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Introduction 



Les moyens de production ont une performance qui s’est considérable- 
ment accrue ces dernières années (en facilité de réglage et en faibles dis- 
persions). Parallèlement, du fait d’une offre abondante du marché, les 
exigences des clients pour une demande en technologie de pointe sont 
de plus en plus fortes. En conséquence, la conformité des produits ne 
peut pas encore être garantie directement sans mesure. Pour ce faire, les 
entreprises utilisent un nombre important de matériels de mesure 
destinés : 

* à prouver la conformité finale des produits vis-à-vis des exigences 
clients (contrôle de réception) ; 

* à vérifier la faisabilité d’un processus de production (en phase de 
mise au point et en suivi de production). 

Réaliser une mesure n’est pas un acte anodin, car elle permet de prendre 
une décision technique et économique. Pour que cette décision soit 
prise à bon escient, il faut avoir foi dans sa mesure. Contrairement à de 
mauvaises habitudes, cette foi ne doit pas être aveugle, il faut savoir 
garder un esprit critique. Cet ouvrage se propose d’éclairer le processus 
de mesure en abordant les thèmes suivants : 

comprendre les exigences de mon client. Cela suppose que nous 
utilisons un même langage. Les normes de définition du produit 
(GPS 1 ) ne seront pas abordées dans cet ouvrage. Par contre, le con- 
cept d’unités « étalon » sera exposé avec quelques rappels historiques 
(chapitre 1) ; 

être capable de choisir un appareil bien adapté à cette mesure. 

Sans faire de la technologie métrologique, cela nécessite de connaître 
les différents principes de mesure et le vocabulaire de description y 
afférents (chapitre 2) ; 


1. GPS : NF EN ISO 14253-1 mars 1999, Spécification géométrique du produit (GPS) - 
Vérification par la mesure des pièces et équipements de mesure, Partie 1 : règles de 
décision pour prouver la conformité ou la non-conformité à la spécification. 
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s’assurer que les appareils utilisés ont bien la qualité métrolo- 
gique requise. Pour cela ils doivent être identifiés et accrédités par 
un service interne ou un organisme externe qui prouvera leur ratta- 
chement aux étalons nationaux (chapitre 3) ; 

analyser les sources d’erreurs potentielles pouvant influencer la 
mesure. Dans l’esprit d’une analyse 5M 1 , nous chercherons à mettre 
en évidence les différents acteurs perturbateurs de la mesure 
(chapitre 4). Nous mettrons l’accent sur l’influence du posage que 
nous retrouvons dans beaucoup de mesurage et comment y remédier 
(chapitre 5) ; 

valider un processus de contrôle. Pour cela nous devons quantifier 
l’erreur de justesse et de fidélité (répétabilité et reproductibilité). 
Après quelques rappels de statistiques de base appliquées à la mesure 
(chapitre 6), nous présenterons comment estimer la fidélité en rap- 
pelant la méthode industrielle du MSA 2 . Ensuite, nous préciserons 
comment estimer l’erreur de justesse et nous présenterons une syn- 
thèse avec l’approche « incertitude inertielle », permettant de lier ces 
deux concepts (chapitre 7) ; 

estimer l’incertitude de mesure associée au résultat donné. Nous 
rappellerons la méthode du GUM 3 permettant d’estimer une incerti- 
tude de mesure par expertises et statistiques. Nous l’appliquerons aux 
mesures avec correction (chapitre 8) ; 

maintenir cette qualité métrologique au cours du temps. Pour 
cela, on présentera différentes méthodes de détection de dérive éven- 
tuelle. Cette surveillance du matériel sera abordée au chapitre 9. 


1 . Duret D., Pillet M., Qualité en production : de l’ISO 9000 à Six Sigma , Paris, Éditions 
d’Organisation, 2005. 

2. MSA : Measurement Systems Analysis, manuel développé sous la responsabilité de Daimler 
Chrysler Corporation, Ford Motor Company et General Motors Corporation 
{www. carwin. co. uk/qs) . 

3. GUM : Guide to the Expression ofUncertainty in Measurement (NF ENV ISO 14253-2 
décembre 2002). 
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Système d'unités 



Système d'unités 


Le besoin d'étalons 

Mesurer, c’est faire une division. Le numérateur représente une grandeur 
(longueur, masse. . .) que l’on veut quantifier. Le dénominateur représente 
une certaine quantité, de même nature physique, que l’on a pu (ou que 
l’on a su) bien définir. C’est la notion même d’étalon (qui correspond 
dans ce cas à la grandeur unité 1 ), que nous développerons par la suite. 
Illustrons ce concept avec un exemple didactique : mesurer un tas de 
sable en nombre de brouettes. 

Pour faire une mesure, il faut : 

► avoir une grandeur physique à mesurer (ici un volume de sable) ; 

► choisir un domaine de cette grandeur comme étalon de mesure 
appelé également grandeur unité (attention ce n’est pas la brouette 
qui est l’étalon, mais le volume contenu dans la brouette) ; 

* dénombrer la quantité à retenir en précisant le nombre (ou mesure) 
d’unités définies par la quantité étalon. 

Figure 1.1 : Exemple didactique de mesure 



Mesure = ■ 



© 


1 . Les termes spécifiques à la métrologie (comme grandeur ou unité) sont définis dans le 
VIM ( Vocabulaire international de la métrologie ) - PR NF ISO 99999 octobre 2006 
(VIM), 3 e éd. 
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La mesure, en tant que résultat de la division, s’exprime par une gran- 
deur scalaire (un réel qui est estimé par un entier et une partie fraction- 
naire). 


Exemple 

Pour mesurer le temps qui passe (chronométrie), le choix de l'étalon cor- 
respond à la période d'un mouvement oscillatoire (du pendule à la 
vibration d'un atome). Il faut ensuite prévoir un système de comptage. 

Nous verrons par la suite que le choix de l’étalon est quelque chose de 
délicat. Sa définition et son appropriation ne sont pas toujours évidentes. 
Cela a conduit à une multitude d’étalons concurrents pour une grandeur 
physique donnée. C’est pourquoi nous ferons toujours suivre le résultat 
d’une mesure par le type de l’étalon retenu (par exemple : 2,35 mm). 

Nous voyons apparaître deux sources d’erreurs potentielles : 

* une définition faible de l’étalon (grandeur unité) ; 

un dénombrement incorrect dans le résultat de la division. 


Un peu d'histoire métrologique 

On peut retenir quatre grands domaines d’application nécessitant une 
métrologie de qualité : 

► domaine commercial (sous-entendu, être sûr d’en avoir pour son 
argent) ; 

» domaine architectural et cadastral (pyramides, cathédrales, palais. . .) ; 

► domaine scientifique (échanger et comparer ses résultats pour valider 
un modèle) ; 

» domaine industriel (passer de l’unitaire à la série, avec la possibilité 
d’interchangeabilité) . 

Concernant le domaine scientifique, rappelons les propos de Lord Kelvin : 

« Si vous pouvez mesurer ce dont vous parlez, et l’exprimer par un nombre, 
alors vous connaissez quelque chose de votre sujet. Si vous ne le pouvez, 
votre connaissance est d’une bien pauvre espèce et bien incertaine. » 
[William Thomson (1824-1907), mieux connu en tant que Lord Kelvin.] 
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A contrario, ce n’est pas parce que l’on a retenu un système de mesure, 
que l’on connaît parfaitement la grandeur à mesurer (par exemple, 
les tests de mesure du quotient intellectuel proposés par Binet (1857- 
1911)). 

L'héritage de l'Antiquité 

L’idée de choisir un domaine d’une grandeur [par exemple la lon- 
gueur] comme étalon est naturel et très ancien (le plus souvent lié à la 
morphologie humaine). On peut citer par exemple, les systèmes 
d’unités et de leurs multiples et sous-multiples égyptiens, gréco- 
romains, arabes et chinois. Ces derniers ont très tôt utilisé un système 
décimal. 


Figure 1.2 : La coudée royale égyptienne 



La métrologie et le pouvoir 

Pour quantifier les échanges commerciaux, il fallait bien utiliser des éta- 
lons de mesures (par exemple, étalons de volume pour la mesure des 
céréales). Dans les transactions commerciales, il pouvait y avoir des 
achats dans une unité et des ventes dans une autre unité plus avanta- 
geuse. La lutte contre ces pratiques injustes relevait naturellement du 
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garant de la justice dans le pays, c’est-à-dire le roi lui-même. Dans la 
plupart des pays, se trouvait un système de « poids et mesures » placé 
sous la haute autorité du roi. 

En 1671, Colbert avait cherché à rendre uniforme les poids et mesures 
dans tous les ports et arsenaux de France 1 2 . 

Notons que cette volonté d’unicité des différents étalons intéressait éga- 
lement les services de collecte des taxes et des impôts. 

Le système anglais 

Ce système, fondé sur des multiples qui ne nous sont plus habituels, est 
en fait très ancien. C’est un des premiers exemples de normalisation 
étendue à tout un royaume. 

Par exemple, pour les longueurs, les multiples et sous-multiples sont 
construits à partir du « pied » . Effectivement, cela correspond approxi- 
mativement à la longueur d’un pied humain (de grande taille), mais cela 
peut être vu également au sens d’instrument de mesure (comme pour un 
pied à coulisse) (tableau 1.1). 


Tableau 1.1 : Système anglais pour les longueurs (exemples) 2 . 


Nom anglais 

Nom français 

Ratio/ pied 

Équivalence SI 

Mile 

Mille 

5 280/1 

1 609,344000 m 

Yard 

Verge 

3/1 

0,914000 m 

Foot 

Pied 

1 

0,304800 m 

Inch 

Pouce 

1/12 

25,400 mm 

Point 

Point 

1/144 

0,176388 mm 


I 

e 


1. Hocquet J.-C., La Métrologie historique, Paris, PUF, 1995. ô 

2. http://fr.wikipedia.org. 
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La genèse du système métrique décimal 

En France, on trouvait un système d’unités très proche du système 
anglais. Le problème essentiel venait de la non-uniformité (dans les pro- 
vinces) de la définition du pied et d’une base de division qui pouvait 
varier (par exemple, la toise et le pied n’étaient pas divisés de la même 
façon). 

Le besoin d'unicité de l'étalon 

Bien que l’idée d’unicité ait déjà été proposée, son application ne 
pouvait se faire que suite à une volonté forte, de manière à vaincre les 
habitudes et l’inertie des usagers. Cependant, le développement du 
commerce, de la région vers l’international, va être un vecteur de redéfi- 
nition d’un système plus cohérent. Pour éviter toutes querelles de natio- 
nalisme, l’idée d’une définition universelle des étalons, rattachés à des 
constantes physiques, a germé dans l’esprit de la communauté scienti- 
fique. 

La fin du XVII e siècle voit s’affronter deux écoles. L’une propose de ratta- 
cher l’étalon à l’arc terrestre (Gabriel Mouton, 1618-1694), l’autre 
définit l’étalon comme le tiers de la longueur d’un pendule battant la 
seconde, appelé « pied universel » (Jean Picard, 1620-1682). 

Les deux écoles coexisteront pendant près d’un siècle et donneront lieu à 
différentes expériences (cela permettra d’affiner la définition de l’étalon 
souhaité). Notons que la définition à partir du pendule créait un lien 
entre l’unité de longueur et l’unité de temps. Ce lien sera à nouveau 
d’actualité, quelque trois cents ans plus tard. 

La nécessité d'une volonté politique 

Promouvoir un nouveau système de mesure ne doit pas uniquement 
être une affaire de spécialistes. Cette promotion doit être soutenue par 
une volonté politique forte de changement et des moyens financiers 
associés. Dès 1774, Turgot charge Condorcet (nommé inspecteur des 
Monnaies) de réfléchir au choix d’un étalon et de son usage dans le 
royaume à partir de copies (on ne parle pas encore de chaîne d’étalon- 
nage, mais le concept est né). 
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La Révolution française, et sa soif de changement, sera un catalyseur 
pour promouvoir ce nouveau système. L’uniformisation des « poids et 
mesures » devient un enjeu national. 

L’Académie des sciences est chargée de définir le nouvel étalon de lon- 
gueur. La définition de l’étalon à l’aide d’un pendule ou à partir d’un arc 
du méridien terrestre sera l’objet de controverses importantes. C’est la 
seconde proposition qui l’emportera. En effet, en 1791, l’Assemblée 
nationale accepte les termes du rapport de la commission composée de 
Borda, Lagrange, Laplace, Monge et Condorcet. La longueur du pen- 
dule battant la seconde est rejetée au profit de la dix-millionième partie 
du quart du méridien terrestre. 

De la définition à la matérialisation du mètre 

Il y avait déjà eu différentes mesures d’arcs de méridien terrestre par 
triangulation. Il a été décidé de mesurer par triangulation l’arc de méri- 
dien entre Dunkerque et Barcelone, deux villes situées à 45 degrés de 
latitude environ. Le choix des villes donnait une dimension internatio- 
nale à cette opération. 

L’Assemblée nationale, en juin 1792, va charger deux astronomes, Jean- 
Baptiste Delambre pour la partie nord et Pierre Méchain pour la partie 
sud, de faire ces mesures. Dans le contexte historique difficile de 
l’époque, ce fut une véritable épopée 1 . 

Le 22 juin 1799, l’étalon « mètre » est consacré. Pour des raisons pure- 
ment pratiques, le mètre sera matérialisé ensuite par une règle de platine. 

La Conférence générale des poids et mesures 

Le système métrique décimal du 22 juin 1799 représentant le mètre et 
le kilogramme mettra quelques années à s’imposer. Présenté lors des 
expositions internationales, sa diffusion sera rapide au sein de la com- 
munauté scientifique mondiale (citons en particulier Gauss, Weber, 
Maxwell et Thomson). 


1 . Pour en savoir plus, lire La Méridienne et Le Mètre du monde de Denis Guedj. 
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De ces différentes rencontres et associations internationales de scientifi- 
ques, naquit le Bureau international des poids et mesures (1876) qui tra- 
vaille sous l’autorité de la Conférence générale des poids et mesures 1 
(CGPM), rassemblant des délégués des états membres de la Convention 
du mètre et des états associés à la CGPM, et qui se réunit à Paris tous les 
quatre ans. 

En 1889, la première CGPM adopte un nouvel étalon du mètre, copie 
de l’ancien, mais de forme et d’alliage plus robustes aux influences (de 
même pour le kilogramme). Il représente le mètre à la température de la 
glace fondante. Un nouvel étalon va être adopté pour le temps : la 
seconde. Ces trois unités de base vont former le système MKS. 

En 1901, Giorgi proposera d’ajouter une unité électrique à ce système. 
Cela donnera le système MKSA adopté en 1 946. 

En I960 (onzième CGPM), ce système sera complété par l’introduction 
du kelvin et de la candela et prendra le nom de Système international 
d’unités (SI). La CGPM propose également une nouvelle définition du 
mètre (en accord avec l’ancienne) : 

« Le mètre est la longueur égale à 1 650 763, 73 longueurs d’onde dans le 
vide de la radiation correspondant à la transition entre les niveaux 2pl0 et 
5d5 de l’atome de krypton 86. » 

La définition renoue avec la volonté de ses créateurs, c’est-à-dire essayer 
de le lier à des constantes physiques. 

En 1967, la treizième CGPM rattache la seconde au monde de l’atome : 

« La seconde est la durée de 9 192 631 770 périodes de la radiation corres- 
pondant à la transition entre les deux niveaux hyperfins de l’état fonda- 
mental de l’atome de césium 133. » 

En 1971, la quatorzième CGPM lui ajoutera la mole. 

En 1983, la dix-septième CGPM propose une autre définition du mètre : 
« Le mètre est la longueur du trajet parcouru dans le vide par la lumière 
pendant une durée de 1/299 792 458 de seconde. » 

Le mètre et la seconde sont liés, nous retrouvons l’esprit du pied uni- 
versel défini à partir du pendule battant la seconde. 


1 . http-.Hwww. bipm. org/fr! 'convention! 'cgpm/. 
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D’une manière générale, en fonction des connaissances physiques du 
moment, on cherchera à rattacher les unités fondamentales à des phéno- 
mènes physiques 1 supposés invariants. Seule l’unité de masse (pour 
l’instant) reste fidèle à sa première définition. 

Nous n’avons décrit que quelques résolutions des différentes confé- 
rences générales des poids et mesures. Pour plus d’informations, le lec- 
teur pourra se reporter au site du Bureau International des Poids et 
Mesures (BIPM) 2 . 


Le système international d'unités : SI 3 

Il est divisé en unités SI de base et en unités SI dérivées. Choisir les 
unités de base juste nécessaires et indépendantes fait appel à des connais- 
sances physiques de haut niveau et qui dépendent donc de l’époque où 
elles ont été définies, héritage de l’histoire, avec ses habitudes et son 
inertie face aux changements. Les sept unités de bases seront considérées 
comme indépendantes et uniques. Le choix a forcément une part d’arbi- 
traire, comme l’illustre la figure 1.3. 

Figure 1.3 : Choisir ses unités de base 



I 

2 

1. http://www.bipm.fr/fr/si/si_constants.html. w 

2. http://www.bipm.org/fr/publications/resolutions.html. §■ 

3. Pour plus de détails, le lecteur est invité à télécharger la notice « Le système interna- o 

tional d’unités » BIPM sur le site : http://www.bipm.org/fr/si/si_brochure/general.html. 
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Unités SI de base 

Elles sont au nombre de sept. Leur mise en œuvre pratique est définie 
sur le site Web du BIPM à l’adresse : http://www.bipm.org/jr/si/si_brochure/ 
appendix2/. 


Tableau 1.2 : Unités SI de base 


m Le mètre (unité de longueur) 

Le mètre est la longueur du trajet parcouru dans le vide par la lumière pendant une durée 
de 1 /299 792 458 de seconde. 

kg Le kilogramme (unité de masse) 

Le kilogramme est l'unité de masse ; il est égal à la masse du prototype international du kilo- 
gramme. 

s La seconde (unité de temps) 

La seconde est la durée de 9 1 92 63 1 770 périodes de la radiation correspondant 
à la transition entre les deux niveaux hyperfins de l'état fondamental de l'atome 
de césium 133. 

A L'ampère (unité de courant électrique) 

L'ampère est l'intensité d'un courant constant qui, maintenu dans deux conducteurs parallèles, 
rectilignes, de longueur infinie, de section circulaire négligeable et placés à une distance 
de 1 mètre l'un de l'autre dans le vide, produirait entre ces conducteurs une force égale 
à 2 x 1 0 - 7 newton par mètre de longueur. 

K Le kelvin (unité de température thermodynamique) 

Le kelvin est la fraction 1 /273, 1 6 de la température thermodynamique du point triple de 

mol La mole (unité de quantité de matière) 

1 . La mole est la quantité de matière d'un système contenant autant d'entités élémentaires 
qu'il y a d'atomes dans 0,012 kilogramme de carbone 1 2. 

2. Lorsqu'on emploie la mole, les entités élémentaires doivent être spécifiées et peuvent être 
des atomes, des molécules, des ions, des électrons, d'autres particules ou des groupements 
spécifiés de telles particules. 

cd La candela (unité d'intensité lumineuse) 

La candela est l'intensité lumineuse, dans une direction donnée, d'une source qui émet 
un rayonnement monochromatique de fréquence 540 x 1 01 2 hertz et dont l'intensité 
énergétique dans cette direction est 1/683 watt par stéradian. 
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Unités SI dérivées 1 

De nombre illimité, elles sont le plus souvent définies à partir des unités 
de base (par exemple, la vitesse exprimée en m/s). Vingt-deux d’entre 
elles ont reçu un nom propre, le plus souvent correspondant au nom 
d’un scientifique illustre (par exemple, la force exprimée en newtons). 
Notons que deux unités sans dimension (correspondant à un rapport de 
même grandeur, comme pour un pourcentage) ont un nom propre, à 
savoir le radian (angle plan) et le stéradian (angle solide). 


Préfixes des multiples et sous-multiples 

La onzième CGPM (I960) a défini des préfixes pour les multiples et 
sous-multiples des diverses unités. Au-dessus du kilo et au-dessous du 
milli, on retrouve une base en 10 3 (fig. 1.4). 


Figure 1 .4 : Choix des multiples et sous-multiples 


Préfixes SI 



Nom & symbole 


Le kilogramme en tant qu’unité de base devrait avoir un nom propre 
(comme par exemple, la proposition historique du « grave », qui a laissé 
la place au « gramme » pour revenir au « kilogramme »). En consé- 
quence, on écrira un mg et non un /ikg. Il faut faire attention de ne pas $ 

I 

Q_ 

ô 

1 . http:llu>ww. bipm. org/fr/si/si_brochure/chapter2/2-2/table3. html. 
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confondre les préfixes avec les unités. Ainsi un « mN » correspond à un 
milli newton et non pas à un couple, qui s’exprimera, pour éviter toute 
confusion, en « Nm ». 

Principales règles d'écriture 

* Les symboles d’unités sont en minuscules sauf si elles dérivent d’un 
nom propre (par exemple, 1 A). 

► Le préfixe est sans séparateur (par exemple, 1 jim). 

► Lorsque l’unité est écrite complètement, elle doit être en minuscule 
comme un nom ordinaire (par exemple, 2 ampères). 

Remarque : le litre peut s'écrire I ou L pour éviter toute confusion avec 
le chiffre 1 (CGPM, 1979). Il n'y a jamais eu de Monsieur Litre 1 . 

Unités hors du SI 

Par habitude, certaines anciennes unités résistent. Il est difficile de 
s’approprier une nouvelle unité. Par exemple le radian, dans le langage 
courant, a du mal à s’imposer devant le degré. En revanche, la minute 
d’angle commence à laisser le pas au l/10 e et au 1/1 00 e de degré (calcu- 
latrice oblige). La mesure du temps utilise toujours une base non déci- 
male (heure, minute, seconde). D’autres unités comme l’hectare, le litre, 
la tonne sont toujours utilisées. 

Les unités précitées ont un « usage » accepté par le SI. 

Il existe un grand nombre d’unités en dehors du système SI. Certaines 
sont tolérées pour des raisons de commodité d’usage, d’autres par contre 
ne sont pas recommandées et devraient donc disparaître progressive- 
ment (par exemple, le cheval remplacé par une expression en watts) 2 . 


1. http://www.student.math.uwaterloo.ca/ -*stat231/stat231_01_02/w02/section3/fil. 2.pdf. 

2. http://www. bipm. org/ff/si/si_brochure/chapter4/4-2. html. 
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Dans ce chapitre nous essayerons de classer les appareils de mesure sui- 
vant deux approches : 

la première sera basée sur la réalisation physique de l’appareil 
(méthode directe et méthode indirecte). Nous nous limiterons aux 
principes de base. Pour une approche plus technologique, le lec- 
teur est renvoyé aux ouvrages spécialisés portant sur l’instrumentalisa- 
tion ; 

* la seconde précisera si l’utilisation d’étalons annexes est nécessaire 
(méthode par comparaison ou à lecture autonome). 

Comme pour toutes disciplines, le métrologue doit utiliser un langage 
de communication si possible exempt de toute ambiguïté. Pour cela un 
minimum de normalisation est nécessaire. Nous détaillerons quelques 
termes utilisés dans ce chapitre en nous référant au Vocabulaire interna- 
tional des termes fondamentaux et généraux de la métrologie 1 . Ce docu- 
ment de synthèse est l’émanation de sept organisations mondiales 
(BIPM, CEI, FICC, ISO, OIML, UICPA, UIPPA). 

De la valeur d’entrée (grandeur réelle ou mesurande) à la valeur de sortie 
(mesure), transitent des données intermédiaires qui peuvent prendre 
différentes formes physiques. Nous recenserons les éléments qui peu- 
vent intervenir dans cette « chaîne de mesures ». 


Principe physique de mesure 

Méthode directe 

C’est la méthode la plus ancienne et la plus naturelle. Par exemple, si 
nous mesurons la longueur d’un objet avec un double décimètre, nous 
obtenons directement une mesure de l’objet. Il n’y a pas de transforma- 
tion physique de l’entrée (déplacement). 


1 . Les termes spécifiques à la métrologie sont définis dans le VIM ( Vocabulaire international 
de la métrologie ) - Projet de norme AFNOR X07-001PR (PR NF ISO 99999 
octobre 2006). 
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Exemples de mesures directes 

• Mesure à partir d'une règle (en métal gravé, en verre avec dépôt. . .). 

• Mesure de masse (balance type « Roberval »). 

• Mesure de volumes (pipette graduée). 

• etc. 

Méthode indirecte 

C’est la méthode la plus répandue dans les appareils de mesure actuels. 
La grandeur d’entrée (mesurande) est transformée en une succession de 
grandeurs intermédiaires de nature physique différente jusqu’à la gran- 
deur finale de sortie qui sera quantifiée (mesure). 

Plus particulièrement, cela permet de passer à des grandeurs de type 
électrique qui facilitent le transfert vers un traitement informatique. 

Exemples de mesures indirectes pour les longueurs 

Le déplacement (entrée) pourra être transformé par exemple : 

• en rotation (micromètre 1 ) : bien que l'utilisateur ait l'impression de 
lire directement un déplacement, il compte en réalité un nombre de 
tours et en déduit indirectement le déplacement du fait du pas de la 
vis micrométrique (fig. 2.1) ; 

Figure 2.1 : « Entrée-Sortie » du micromètre 
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• en ondes sonores (par exemple la conformité dimensionnelle d'une 
pale de turbine sera donnée suite à l'analyse du spectre du son pro- 
voqué par sa vibration) ; 

• en variations de pression (analyse des pressions des fuites d'air, fonc- 
tion du jeu entre l'appareil et la pièce à mesurer) ; 

• etc. 

Certaines mesures de grandeurs nécessitent une méthode indirecte 
comme par exemple : 

► la mesure des forces (dynamomètre de type « peson », jauge de con- 
traintes, anneau de déformation avec mesure de pression d’air, cap- 
teur piézoélectrique. . ; 

► la mesure des températures (dilatation de liquide, solide ou gaz, 
variation de résistance électrique, thermocouple, effet photovol- 
taïque, rayonnement. . .). 

Méthode directe globale avec lecture indirecte locale 

En métrologie dimensionnelle, la lecture directe par règle est très 
employée. Il suffit de compter directement un nombre de pas. Il reste le 
problème de mesurer une fraction de pas. Cette fraction peut être 
mesurée directement par des moyens de grossissement (type micros- 
cope) ou en utilisant une lecture par amplification géométrique de type 
vernier 1 (fig. 2.2) ou franges de Moiré. 

Figure 2.2 : Exemple d'une lecture (30,80 mm) avec vernier au 1/50 e 


© 



1. Inventé par le géomètre français Pierre Vernier (1580-1637). 
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Lors de l’utilisation de règles optiques, on peut compter le nombre de 
pas (réseau de traits distants de 0,25 à 10 flm environ, correspondant à 
des couches de métal sur un support en verre ou en céramique) et déter- 
miner la fraction de pas complémentaire par analyse du signal électrique 
correspondant. 


Étalons externes ou internes 

Mesure par comparaison (étalons externes) 

C’est le type même de l’appareil utilisé en grande série. En effet, dans ce 
cas particulier, l’appareil sera utilisé toujours pour la même tâche. On 
lui demandera de permettre l’estimation d’écarts par rapport à une pièce 
supposée connue (étalon de travail), sensiblement de même dimension 
(voisinage de mesure). L’étalon de travail peut être une pièce particulière 
de la production que l’on a retenue et promue au rang d’étalon après 
l’avoir mesurée finement, ou une pièce de forme géométrique plus 
simple (cales-étalons, cylindres étalons...). 

En réalité, l’appareil a de fait un étalon embarqué (interne), mais il est 
de faible étendue. Il doit simplement permettre de couvrir l’étendue des 
écarts par rapport à l’étalon de travail. Ces appareils sont souvent 
appelés « comparateurs ». Par exemple, dans un comparateur méca- 
nique d’étendue de mesure d’environ 10 mm, la crémaillère taillée dans 
la tige joue ce rôle d’étalon interne. 


I 

e 


o 
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Figure 2.3 : Exemple de mesure de la longueur d'une pièce 
avec un comparateur mécanique et une cale-étalon 



Dans l’exemple de la figure 2.3, l’incertitude sur le résultat de la mesure 
va dépendre entre autres de la qualité de la cale-étalon et du compara- 
teur. 

Mesure autonome (étalons internes) 

Dans ce cas de mesure, l’appareil devra fournir un résultat seul. Cela 
nécessite un système d’étalons propres, embarqué dans la structure de 
l’appareil. Ce système doit permettre de couvrir l’étendue de mesure 
d’utilisation de l’appareil. 

Exemples 

• Pied à coulisse (règle gravée). 

• Micromètre (vis micrométrique). 

• Machine à mesurer tridimensionnelle (règle optique). 
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Terminologie 

Nous préciserons simplement quelques termes utilisés régulièrement 
dans cet ouvrage. Pour la définition normalisée des différents termes, le 
lecteur est renvoyé au Vocabulaire international des termes fondamentaux 
et généraux de la métrologie 1 . 

Mesurande : valeur vraie de la grandeur à mesurer (par définition 
cette valeur ne sera jamais exactement connue). 

Mesure : valeur estimée du mesurande suite à un mesurage. 
Mesurage : activités permettant d’obtenir une mesure. 

Étalon : entité physique censée représenter en un lieu donné, le 
mieux possible, une fraction de l’unité de base. Au niveau de l’entre- 
prise, on distinguera : 

- les étalons de référence : matérialisant la plus haute qualité 
métrologique de l’entreprise (généralement situés en métrologie 
centrale), 

- les étalons de travail : utilisés pour étalonner ou vérifier les appa- 
reils utilisés sur les sites de production. Ils sont eux-mêmes éta- 
lonnés par rapport aux étalons de référence (ou aux étalons de 
transfert quand il y a plusieurs niveaux d’étalons). 

Intervalle de mesure : correspond à l’ensemble des valeurs que peut 
prendre la grandeur à mesurer. Par exemple, pour un thermomètre 
[— 20 °C, + 80 °C], cela donne une étendue de mesure de 100 °C. 
Répétabilité : si la mesure est répétée plusieurs fois par le même opé- 
rateur, dans les mêmes conditions (même appareil et même labora- 
toire), dans un laps de temps court, l’estimation de la variance à 
partir des résultats permet de quantifier ce concept de répétabilité. 
Justesse : la justesse d’une méthode de mesure est estimée par la dif- 
férence entre une valeur vraie de la grandeur mesurée et la moyenne 
des différents résultats de mesure, elle est exprimée en termes de 
biais. 


1. VTM {Vocabulaire international de la métrologie) — Projet de norme AFNORX07-001PR 
(PR NF ISO 99999 octobre 2006) . 


I 
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Sensibilité : la mesure S (sortie) est fonction du mesurande E 
(entrée). Soit S =/(ü) la fonction correspondante. La sensibilité sera 

égale à la dérivée de cette fonction : s — . Si s est une constante, 

l’appareil aura un comportement linéaire. 

Chaîne de mesure : suite des éléments de l’appareil conduisant le 
signal de mesure de l’entrée (capteur) jusqu’à la sortie (communica- 
tion avec l’utilisateur). 

Résolution : la plus petite différence d’indication d’un dispositif affi- 
cheur qui peut être perçue de manière significative ( VIM ). La résolu- 
tion est un indicateur de qualité de l’appareil, mais il ne faut pas se 
focaliser uniquement sur cet aspect et regarder le rapport résolution/ 
étendue (pouvoir de résolution). Nous obtenons par exemple : 

- pied à coulisse : yyj = 1,3-10 4 , 

- micromètre : = 4 • 10~ 4 , 

25 

- comparateur mécanique (course 10 mm) : — = 1-10 3 . 


Notions de chaîne de mesure 

L’appareil peut très bien être utilisé sans que l’on sache comment est 
traité le signal de mesure. On ne s’occupe alors que de la fonction de 
transfert entre l’entrée et la sortie (fig. 2.4). 

Figure 2.4 : L'appareil de mesure considéré 
comme une « boîte noire » 


1 

^ Grandeur d'entrée 

Fonction 

Grandeur de sortie 

Q_ ^ 

o Mesurande 

(J 

© 

de transfert 

Résultat de mesure 
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La mesure n’est qu’une interprétation subjective de la réalité (mesu- 
rande). Pour être conscient de cette subjectivité (et éventuellement être 
capable d’esprit critique), il est souhaitable d’avoir une vue de 
l’ensemble des composants de la chaîne de mesure. Cette dernière peut 
être illustrée par les composants suivants (fig. 2.5). 

Figure 2.5 : Exemples d'éléments d'une chaîne de mesure 



Capteur 

C’est l’élément clef de la chaîne de mesure qui donne le signal de base. 
Les variations du mesurande doivent être « captées » de manière signifi- 
cative afin d’être « interprétées » par la suite. Cette saisie d’informations 
ne peut se faire souvent qu’avec un apport d’énergie extérieure (énergie 
d’activation). Par exemple, avec un pied à coulisse, il faut exercer une 
pression au niveau des becs pour garantir une bonne saisie. Cet apport 
d’énergie est nécessaire, mais peut également être néfaste en perturbant 
le mesurande que l’on veut quantifier (chapitre 4). 
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La conversion physique de la nature du signal commence souvent au 
niveau du capteur, par exemple convertir une grandeur mécanique en 
une grandeur électrique (fig. 2.6). 


Figure 2.6 : Capteur LVDT 
( Linear Variable Differential Transformer ) 



Amplificateur d'instrumentation 

Le niveau énergétique dans l’environnement du capteur est souvent très 
faible et fortement influencé par les perturbations. Pour pouvoir com- 
muniquer avec l’utilisateur (enregistrement graphique, indicateur...), 
cela nécessite un seuil énergétique minimal. L’amplificateur doit fournir 
un nouveau signal le plus fidèle possible à l’original. Ce gain d’énergie 
est souvent obtenu à l’aide du réseau électrique ou pneumatique dispo- 
nible dans l’entreprise. 

Traitement du signal de mesure 

Nous nous limiterons dans ce chapitre à indiquer deux modes de 
traitement : le filtrage et la conversion analogique/numérique. Pour plus 
de détails, le lecteur est convié à se reporter à des ouvrages spécialisés. 
L’analyse et le traitement mathématique des données seront repris dans 
les chapitres suivants. 
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Filtrage 

Le filtrage est nécessaire lorsque le signal de mesure fourni est trop riche. 
Ce surplus d’information peut être vu comme une pollution vis-à-vis de 
l’objectif recherché. Il sera alors nécessaire de trier l’information fournie. 
Parfois, c’est l’inverse qui se produit, des informations utiles sont igno- 
rées par le capteur qui ne les perçoit pas. 

La mesure de la rugosité illustre bien cet aspect. La surface réelle palpée 
comporte des défauts de forme (premier ordre), des défauts d’ondula- 
tion (deuxième ordre), des défauts de rugosité (troisième ordre) et pour 
finir, des défauts d’ordre supérieur (arrachements par exemple). La tech- 
nologie employée va engendrer de fait un filtrage. Par exemple, un pal- 
page mécanique est tributaire de la forme du palpeur (rayon de pointe) 
et du type d’exploration de la surface (palpage avec ou sans patin). En 
fonction du critère que l’on veut mesurer (rugosité par exemple), cela 
nécessitera de choisir le filtre le plus adapté (par exemple : emploi d’un 
filtre gaussien surfacique ISO 16610-61) (fig. 2.7). 

Figure 2.7 : Exemple de représentation de rugosité 


Hm 



Convertisseur analogique numérique 

Le plus souvent le capteur traduit le signal d’entrée physique (déplace- 
ment, pression, température...) en un signal électrique de faible puis- 
sance sous forme analogique (signal continu). Pour des raisons 
d’affichage, de stockage d’information, de traitement par un calculateur 
il sera nécessaire de transformer cette grandeur analogique en une gran- 
deur numérique à l’aide d’un codage (par exemple, code binaire naturel, 
code de GRAY, code DCB. . .) (fig. 2.8). 
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Figure 2.8 : Exemple de conversion analogique/numérique 



0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 


Entrée (analogique) 


Affichage - mémorisation de la mesure 

L’appareil de mesure peut communiquer avec l’utilisateur sous diffé- 
rentes formes. 

► La mesure est du type grandeur scalaire : un affichage à aiguille ou 
numérique est suffisant. 

► La mesure est multi-composantes : 

- affichage à plusieurs cadrans (type multi-cotes), 

- affichage sous forme de graphe (diagramme de circularité sous 
forme de coordonnées polaires, rugosité. . 

- affichage sous forme de zones (relevé de température par thermo- 
graphie), 

- écriture sous forme de fichier (fichier de points sur machine à 
mesurer tridimensionnelle) avec traitement informatique associé. 

Les enregistrements (physiques ou informatiques) permettent de conserver 
une mémoire de cette information et de permettre un traitement ulté- 
| rieur. 

Q- 

o 
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Chapitre 3 


Gestion des moyens 
de mesure 



Les entreprises utilisent beaucoup d’appareils pour garantir la confor- 
mité de leurs produits et pour évaluer la performance de leurs processus 
(conception et production). Pour cela, il faut avoir confiance dans les 
mesures obtenues. Les référentiels de systèmes qualité, de type ISO 
900 1 , exigent que les appareils soient vérifiés régulièrement pour pou- 
voir assurer cette qualité de mesure. Cela demandera une gestion admi- 
nistrative et technique du parc d’appareils de contrôle. 

La mise en œuvre de la gestion et de la maîtrise des moyens de mesure, 
doit s’accompagner d’une réflexion plus globale de la « fonction 
métrologie », portant sur la stratégie de suivi (gestion externe ou 
interne) et, d’une manière plus générale, sur l’importance d’une culture 
métrologique dans l’entreprise. 


Situer sa « fonction métrologie » 

Pour cette réflexion, il est souhaitable de répondre aux questions sui- 
vantes. 

Est-ce que la qualité requise est bien spécifiée par mon client ? 

Des disparités très fortes peuvent se rencontrer. On trouve des secteurs 
industriels où le produit est défini avec une grande rigueur. Par 
exemple, dans le secteur automobile où l’on peut avoir plusieurs four- 
nisseurs (de langue et de culture différentes) pour un même produit, 
cela nécessite de préciser les variations de qualité acceptables en utilisant 
le mieux possible les normes internationales. Dans les entreprises de 
taille moyenne, où les études et la production cohabitent, il y a risque 
d’utiliser un langage de définition interne, pseudo-normalisé, engen- 
drant des ambiguïtés. Enfin à l’extrême, il reste des secteurs où la des- 
cription de la qualité relève du bon usage. Ce dernier cas est dangereux, 
car il engendre une subjectivité de décision de conformité, source d’une 
variabilité importante. 
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Est-ce que je comprends les exigences de mon client ? 

Cela sous-entend que l’on parle le même langage et que celui-ci est suffi- 
samment riche pour être le plus explicite possible. En toute rigueur, une 
définition de domaine de conformité ne devrait pas conduire à deux 
interprétations différentes. Par exemple, la lecture des spécifications 
géométriques du produit 1 nécessite de bien connaître la normalisation 
se rapportant à ce sujet. L’exigence à vérifier doit conduire au choix d’un 
appareillage de mesure approprié. Dans l’exemple suivant, la mesure 
d’une cale présentant deux défauts angulaires a jet (X 2 , ne donnera pas le 
même résultat suivant la technique de mesure utilisée. La décision de 
conformité exige de bien savoir ce que l’on veut réellement mesurer et 
ne pas se contenter de prendre l’appareil disponible (fig. 3.1). 

Figure 3.1 : Extraction de longueur avec machine à mesurer monodimensionnelle 
et colonne de mesure 



Quels moyens de mesure retenir et qui sera chargé de faire ces mesures (dispo- 
nibilités, compétences) ? 

La mesure n’apporte pas de valeur ajoutée physique au produit. À court 
terme, cela peut être perçu comme un investissement non prioritaire 


1. NF EN ISO 14253-1 mars 1999, Spécification géométrique du produit (GPS) - Véri- 
fication par la mesure des pièces et équipements de mesure - Partie 1 : Règles de déci- 
sion pour prouver la conformité ou la non-conformité à la spécification. 
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(souvent élevé) par rapport à des ressources de production. C’est un peu 
comme une assurance, l’investissement en ressources humaines et maté- 
rielles doit être mis en regard avec le risque de détection de produit non 
conforme chez le client. C’est un investissement rentable à long terme, 
cependant il faut éviter toute surenchère. Parfois, c’est le client, du fait 
de son poids économique, qui exigera tel ou tel moyen de contrôle. 
Dans le cas de transactions concernant les domaines du commerce, de la 
santé, de la sécurité et de l’environnement, les autorités publiques exi- 
gent une garantie de la qualité de la mesure, on parlera alors de métro- 
logie légale 1 . 

Bien répondre à cette question est une source d’économie importante. 

Comment s’assurer de la performance requise pour ces appareils, dans leur 
contexte de mesure ? 

Il s’agit de valider un processus de mesure. La performance du moyen 
sera tributaire des composantes humaine et technique et perturbée par 
les influences du milieu dans lequel la mesure se réalise. Le résultat de 
mesure n’est qu’une image de la réalité. Pour que les décisions découlant 
de cette mesure soient fondées, il faut avoir une perception de la robus- 
tesse de cette dernière. 

Comment vérifier qu ’ils ne se dégradent pas dans le temps ? 

Comme dans tout processus, des dérives en fonction du temps risquent 
d’apparaître. Il sera nécessaire de définir une fréquence de surveillance 
pour garantir que la performance de l’équipement ne s’est pas dégradée. 

Commet gérer cette information ? 

À ces questions techniques, il faut rajouter la partie administrative rela- 
tive à la maîtrise des enregistrements pour prouver que ces vérifications 
ont bien eu lieu et ont été réalisées par un personnel compétent et accré- 
dité. 

Tout cela dans un environnement économique souvent difficile. Il faut 
le juste nécessaire pour éviter toute surqualité. 


1 . Organisation internationale de la métrologie légale (OIML : www. oiml. org) . 
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La confiance dans la mesure 

La qualité de la mesure est définie par une approche de type 5M (cette 
analyse sera reprise en détail au chapitre suivant) (fig. 3.2). 


Figure 3.2 : Exemple de démarche 5M pour une mesure de rugosité 
Méthodes Milieu Main-d'œuvre 



Le MSA 1 retient une arborescence légèrement différente, mémorisée 
à l’aide de l’acronyme : SWIPE {Standard, Work Piece, Instrument, 
Person/Procedure, Environment). 


L'instrument 

Par instrument, il faut comprendre l’appareil, les équipements et les 
montages associés ainsi que les étalons utilisés. 

Toutes les entreprises n’ont pas le besoin (ni les moyens) d’être au 
niveau de la définition du mètre étalon. Une copie de bonne qualité 
d’un sous-multiple, sous une forme physique plus accessible (boîte de 
cales par exemple), est suffisante, sous réserve qu’elle soit garantie par un 
organisme mandaté. 

L’entreprise doit s’assurer que la maîtrise de la mesure est déclinée 
jusque sur les sites de production. Chaque site n’a pas forcément besoin 


1 . Measurement Systems Analysis : manuel développé sous la responsabilité de Daimler Chrysler 
Corporation, Ford Motor Company et General Motors Corporation {www.carwin.co.uk/qs). 
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des mêmes performances métrologiques. Il suffit qu’elles soient en adé- 
quation avec la qualité requise par la définition du produit (sauf si 
l’entreprise possède l’appareil et désire tout simplement le rentabiliser). 
Pour cela on utilise des appareils spécifiques ou standard qui devront, à 
une fréquence donnée, être comparés aux étalons de référence (c’est la 
plus haute qualité métrologique existante dans l’entreprise). Si la qualité 
métrologique exigée au niveau de l’appareil de mesure en production 
n’est pas très grande, on peut se contenter de le tester à l’aide d’étalons 
intermédiaires, sous réserve qu’ils soient vérifiés par rapport aux étalons 
de référence. Ces étalons sont appelés « étalons de travail ». On utilise 
parfois des « étalons de transfert » comme maillon secondaire dans la 
comparaison avec les étalons de référence. La gestion et la mise en place 
de ces tests et comparaisons constituent une chaîne de mesure interne 
(fig. 3.3). 

Exemple 

Imaginons une grande entreprise, implantée sur plusieurs sites. Elle pos- 
sède une métrologie centrale au service des autres sites, une métrologie 
générale et des métrologies d'atelier dans chaque site respectif 
(fig. 3.3). Le raccordement aux étalons nationaux sera développé au 
paragraphe suivant (chaîne nationale d'étalonnage). 

Figure 3.3 : Notion de chaîne d'étalonnage interne 
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Toutes les entreprises n’ont pas les moyens techniques et financiers pour 
mettre en place une chaîne de mesure aussi conséquente. Elles peuvent 
faire vérifier (ou étalonner) tout ou partie de leurs appareils de mesure 
par un fournisseur de prestations métrologiques apportant les garanties 
nécessaires de liaison aux étalons nationaux. Rappelons le paragraphe a 
du chapitre « 7.6 Maîtrise des dispositifs de surveillance de mesure » de 
la norme ISO 9001 1 : 

« [...] Lorsqu’il est nécessaire d’assurer des résultats valables, les équipe- 
ments de mesure doivent être : 

a) étalonnés ou vérifiés à intervalles spécifiés ou avant leur utilisation, par 
rapport à des étalons de mesure reliés à des étalons de mesure internatio- 
naux ou nationaux (lorsque ces étalons n’existent pas, la référence utilisée 
pour l’étalonnage doit faire l’objet d’un enregistrement) ; 

[•••]» 

Le milieu 

Avec l’instrument, c’est le second facteur qui vient spontanément à 
l’esprit, lorsque l’on rencontre des problèmes de qualité de mesure. En 
particulier, en cas de problèmes de métrologie dimensionnelle, 
l’influence de la température est souvent perçue comme la cause 
majeure, alors que la contribution des autres « M » est peut-être plus 
importante et nettement moins coûteuse à corriger. 

L’influence du temps peut engendrer une augmentation de la dispersion 
des mesures, mais également une dérive (fig. 3.4). Elle est également 
appelée constance ou stabilité. 


I 

e 


l . Reproduction faite avec l’autorisation d’AFNOR. 
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Figure 3.4 : Influence du temps (constance) 



La main-d'œuvre et la méthode 

Ces deux « M » nous semblent inséparables. En plus de l’aptitude phy- 
sique (dextérité, sens expérimental), une formation scientifique et tech- 
nique est nécessaire pour faire des mesures de qualité. Elle s’appuie sur 
des connaissances technologiques et normatives. En effet, si l’on est 
censé valider une conformité, cela va de soi que l’on comprend ce que 
demande le client (par exemple avec l’emploi des normes GPS). 
L’auto-maîtrise dans un processus de production impose de normaliser 
le sous-processus de mesure (appelé également protocole de mesure). 
Dans bien des secteurs, souvent avec une main-d’œuvre très compé- 
tente, ce point est négligé, car cela semble faire offense aux opérateurs. 
Pourtant suite à des études GRR (chapitre 7), on constate que des diffé- 
rences de comportement mineures, a priori sans influence, engendrent 
une variabilité de mesure forte. D’où l’importance de la communication 
entre les différents opérateurs sur les méthodes de mesure retenues (réu- 
nions, supports visuels, instructions. . .). 
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La matière (pièce) 

L’incidence de la pièce peut provenir de : 

» ses qualités physiques (par exemple de la souplesse du matériau) ; 

► de l’interaction « opérateur-matière ». Par exemple, un opérateur 
recherchera systématiquement le diamètre maximum pour une pièce 
ovalisée, et ceci de façon intuitive ; 

» de sa définition de conformité lorsque des caractéristiques sont liées. 

Il est dangereux de chercher à quantifier séparément l’influence de 
chaque « M ». Une démarche type « Analyse de la variance » peut être 
retenue pour voir s’il y a des interactions (chapitre 6). 


Chaîne nationale d'étalonnage 

La traçabilité des raccordements, fruit de deux organismes, permet : 

► de retrouver les enregistrements d’étalonnage à différents niveaux, de 
l’entreprise au niveau national, c’est la garantie du COFRAC ; 

» de garantir le lien physique avec les étalons nationaux, c’est le rôle du 
LNE. 


Le Comité français d'accréditation (COFRAC 1 ) 

Le rôle du COFRAC est d’accréditer les organismes délivrant des certifi- 
cations. L’accréditation d’un laboratoire est la reconnaissance d’une cer- 
taine compétence dans un domaine précis (par exemple, vérification des 
tampons filetés de diamètre 5 mm à 20 mm) et l’assurance de la péren- 
nité de cette compétence par un organisme faisant lui-même autorité en 
la matière (par exemple, système de management de la qualité suivant le 
référentiel NF EN ISO/CEI 17025 : septembre 2005 2 pour le labora- 
toire accrédité). 


1 . http:llwww. cofrac.fr. 

2. NF EN ISO/CEI 17025 septembre 2005, Exigences générales concernant la compé- 
tence des laboratoires d’étalonnages et d’essais. 
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En simplifiant, on peut dire que l’accréditation « COFRAC » apporte la 
garantie que les laboratoires retenus (possédant cette accréditation) sont 
fournisseurs de prestations de qualité car ils possèdent une technologie adé- 
quate et qu’ils mettent en œuvre un suivi documentaire prouvant le lien 
métrologique des étalons du client aux étalons nationaux. Cette accrédita- 
tion n’est pas imposée par l’ISO 9001, cela simplifie l’audit 1 pour la 
partie métrologie. 

Il existe une reconnaissance mutuelle entre divers organismes d’accrédi- 
tation 2 , c’est-à-dire que l’on estime que les écarts constatés entre ces 
laboratoires ne sont pas significatifs. 

Laboratoire national de métrologie et d'essais (LNE 3 ) 

Pour assurer la confiance dans l’appareil utilisé, celui-ci doit être éta- 
lonné. Cet étalonnage se fera souvent à l’intérieur de l’entreprise en réfé- 
rence aux étalons de la métrologie centrale (étalons de référence). 

Cette confiance ne sera totale que si l’on peut valider les étalons de 
l’entreprise. Le LNE a mis en place les structures et procédures permet- 
tant d’assurer la valeur technique et scientifique des étalonnages et de 
fournir une garantie officielle. 

Le Laboratoire national de métrologie et d’essais (LNE) remplace 
depuis le 27 janvier 2005 le Bureau national de métrologie (BNM), 
organisme interministériel, qui avait été créé en 1969. Fédérant les labo- 
ratoires de recherche détenant les étalons nationaux, le LNE a une mis- 
sion de soutien technique et scientifique à l’industrie et au monde de la 
recherche. Il contribue également au développement des programmes 
de métrologie européens et internationaux. 

Les quatre laboratoires nationaux sont : 

» LNE (Laboratoire national de métrologie et d’essais) ; 

» LNHB (rattaché au Commissariat à l’énergie atomique) ; 


w 1 . Pour plus de détails, le lecteur pourra consulter L Audit de la métrologie dans le cadre des 

§- certifications de systèmes de management. Collège français de métrologie, septembre 2005. 

ô 2. www.european-accreditation.org. 

® 3. www.lne.fi. 
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» INM (rattaché au Conservatoire national des arts et métiers) ; 

» SYRTE (rattaché à l’Observatoire de Paris). 

Les laboratoires accrédités 

Des organismes intermédiaires sont chargés de divulguer la métrologie à 
differents niveaux. Cela permet une approche plus efficace et plus éco- 
nomique. 

En 1977, on avait précisé les maillons intermédiaires de la chaîne d’éta- 
lonnage entre le BNM et les entreprises en créant : 

► des centres d’étalonnage agréés ; 

► des services de métrologie habilités. 

Dans le souci d’une harmonisation européenne, il n’y aura plus de label 
CEtA ou SMH, mais simplement des laboratoires accrédités par le 
COFRAC. Ces laboratoires, en situation concurrentielle, doivent faire la 
preuve que leur niveau d’incertitude est bien en conformité avec les 
prestations qu’ils fournissent (de petits laboratoires peuvent atteindre, 
dans un domaine particulier, un savoir-faire métrologique très pointu). 
Ils doivent pour cela être raccordés à un laboratoire accrédité d’incerti- 
tude plus faible. Si ce n’est pas possible, ils se raccordent alors directe- 
ment aux laboratoires nationaux ou aux laboratoires associés. 

En simplifiant, on peut illustrer la dissémination de la métrologie à 
l’aide de la figure 3.5. 


Niveau 1 : laboratoires nationaux 
(Etalon primaire) 


Niveau 2 : laboratoires associés 
(Etalon secondaire) 


Niveau n - 1 : laboratoires locaux 
(Étalon de transfert) 


Niveau n : entreprise 
(Étalon de référence 
de l'entreprise) 


Figure 3.5 : Évolution des incertitudes 
Incertitudes 
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L’objectif de cette chaîne de mesure est de quantifier la croissance de 
l’incertitude de l’étalon primaire à l’étalon de référence de l’entreprise. 


Choix d'une stratégie d'étalonnage 

Le rattachement aux étalons nationaux 1 nécessite à partir d’un certain 
seuil une externalisation des étalonnages. Reste le problème de savoir à 
quel niveau la faire. En simplifiant, on peut envisager deux solutions 
(fig. 3.6). 


Figure 3.6 : Exemples de chaînes d'étalonnage 


© 


Bureau international des poids et mesures (BIPM) 


NIST 


LNE 

Autres 

(USA) 


(étalons nationaux) 
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^ 

Laboratoire 


Laboratoire 

accrédité 


accrédité 

(incertitude très faible) 


(incertitude très faible) 


t — 

Laboratoire accrédité 

(possibilité de faire 
de la sous-traitance 

pour l'entreprise A) 


Entreprise B 



1 . http:llwww. metrologiefrancaise.fr. 
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En dehors de l’aspect stratégique de conserver ou non un savoir métro- 
logique d’étalonnage, ce choix peut être dicté par des considérations 
économiques : 

► raccordement externe (facturation de l’ordre de vérification, coût 
d’abonnement...) ; 

► raccordement interne (formations et salaires, coût des équipements, 
locaux adaptés, veilles normatives. . .). 

Les laboratoires (internes ou externes) peuvent être considérés comme 
des fournisseurs. Ils fournissent une prestation métrologique. Pour 
garantir la satisfaction du client (au niveau organisation et compétence 
technique), il doivent mettre en place une démarche qualité. Cette 
démarche s’appuiera sur un référentiel de système qualité reconnu ayant 
fait ses preuves. De nombreuses normes et guides concernent les labora- 
toires, on peut citer en particulier : 

N F EN ISO/CEI 17025 (septembre 2005) : Exigences générales 
concernant la compétence des laboratoires d'étalonnages 
et d'essais 

Bâtir son système qualité sur cette norme, permet de garantir : 

► l’efficacité du management du laboratoire (dans l’esprit de l’ISO 900 1 : 

2000 1 ) ; 

» la compétence scientifique et technique (qualité des ressources 
humaines et matérielles ainsi que des méthodes retenues). 

Avoir son système qualité construit en accord avec ce référentiel sera une 
aide précieuse pour se faire accréditer. 

N F EN ISO 10012 (septembre 2003) : 

Exigences pour les processus et les équipements de mesure 

Au sein de l’entreprise, la norme ISO 9001 : 2000 préconise un mana- 
gement de l’entreprise par une approche processus, avec la volonté de 
satisfaire le client après avoir bien compris ses exigences. L’idée générale 
est d’appliquer la même méthodologie aux processus de mesure qu’aux 
autres processus. On retrouve les mêmes chapitres que dans la norme 


1 . NF EN ISO 900 1 décembre 2000, Systèmes de management de la qualité - Exigences. 
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ISO 9001, seul le chapitre 7 est réintitulé « Confirmation métrologique 
et mise en œuvre des processus de mesure ». On retrouvera entre autres 
dans ce chapitre 7, des recommandations sur la périodicité d’étalon- 
nage, la maîtrise des équipements et logiciels, la conception des appa- 
reils, l’estimation des incertitudes, etc. 


La mesure, un état d'esprit 

Le concept mesure ne concerne pas uniquement l’opérateur métro- 
logue. Ce concept doit être intégré de la conception à l’élaboration du 
produit. En particulier : 

► le bureau d’études ne doit pas définir des paramètres ou des caracté- 
ristiques non quantifiables. Il doit aussi estimer le mieux possible les 
variations acceptables, compatibles avec un fonctionnement correct 
du mécanisme (fonction perte de Taguchi). Cette distribution des 
tolérances doit être faite avec l’aide d’experts en production. Citons 
dans cet esprit la démarche GPS (spécification géométrique des pro- 
duits) qui a comme objectif d’être un langage unique (cohérence des 
différentes normes) pour tous les services, des études au contrôle 
final ; 

► le service méthodes doit vérifier la compatibilité entre les valeurs de 
définition et les capacités des moyens de production et de contrôle 
(notion de capabilité). Suite à cette analyse, « les cotes de fabrication 
à risques » devront être bien répertoriées et clairement indiquées sur 
les plans de surveillance en fabrication. Les autres cotes pourront être 
contrôlées avec une fréquence beaucoup plus faible ; 

» les métrologues de terrain doivent garder un esprit critique (cons- 
tructif) pour l’utilisation des moyens de contrôle et savoir estimer 
leur incertitude de mesure. Tout constat de dérives ou de rebut élevé 
doit être répertorié et l’information retournée aux services de concep- 
tion et de production concernés ; 

J 
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» le service métrologie doit s’équiper (achat ou construction) d’un 
matériel compatible avec les exigences de qualité requises mais en 
évitant un degré de précision de mesure excessif (besoin réel). Gros- 
sièrement, le coût des appareils de mesure est une fonction puissance 
de l’inverse de la résolution. 


Mise en place d'une gestion 
des moyens de mesures 

Dans ce paragraphe, nous étudierons l’aspect purement gestionnaire. 
L’aspect maîtrise de l’aptitude à l’emploi (notion de conformité) des 
équipements de contrôles, de mesures et d’essais sera traité dans les cha- 
pitres suivants. 

Identification et inventaire des appareils 

Rappelons que la norme ISO 900 1 1 demande de : 

« déterminer [...] les dispositifs de surveillance et de mesure nécessaires 
pour apporter la preuve de la conformité du produit aux exigences 
déterminées ». 

A priori, tous les appareils ont une incidence, ils ont été achetés dans un 
but bien précis qui a toujours un lien plus ou moins direct avec le pro- 
duit. 

Pour cette identification, on peut faire des classes correspondant au 
risque encouru sur le produit en cas d’information erronée dans la 
mesure. Par exemple, un ohmmètre, destiné à la maintenance dans une 
industrie de sous-traitance mécanique, ne demande pas un suivi attentif. 
Les appareils non identifiés sont considérés comme des indicateurs (par 
exemple, un réglet pour un préréglage de machine). L’identification 
peut se résumer à un point de couleur précisant que ces appareils ne sont 
pas suivis. 


1 . Reproduction faite avec l’autorisation de l’AFNOR. 
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Le code d’identification contient plus ou moins d’informations (site 
d’affectation, type, service...)- Il permet d’identifier l’appareil (inscrip- 
tion indélébile sur l’appareil si possible) et d’avoir accès à la fiche de vie 
et au dossier le concernant (par exemple : PC 011 pour le 11 e pied à 
coulisse acheté). 

L’inventaire des appareils de mesure de l’entreprise représente une 
somme de travail importante. Dans certaines entreprises, l’appareil de 
mesure est perçu comme un bien précieux et personnel dans lequel on a 
une confiance parfois aveugle. En plus de l’aspect technique, il faut 
expliquer le bien-fondé de la démarche pour lever toute réticence 
d’identification. 

Cette tâche permet de bien cerner les besoins métrologiques de l’entre- 
prise et éventuellement d’éviter des doublons pour les acquisitions 
futures. En particulier, on pourra étudier l’évolution en analysant les 
dates d’acquisition. 

Cet inventaire donne également la possibilité de faire du tri. On con- 
serve parfois des appareils pour une hypothétique et peu probable utili- 
sation. Cela coûte cher en immobilisation et est contraire à la méthode 
des 5S. 

Il arrive que des appareils peu utilisés doivent être maintenus car on ne 
trouve pas une prestation identique sur le marché. De plus, certaines 
entreprises pensent qu’il est stratégique de garder une compétence 
métrologique importante. Comme nous l’avons déjà signalé, il peut être 
opportun alors de fournir des prestations métrologiques pointues à 
l’aide de ce matériel à des tiers. 


Indicateurs de qualité d'un appareil 

La performance métrologique est quantifiée par les notions de justesse et 
de fidélité. L’écart de justesse sera estimé par l’opération d’étalonnage ou 
borné par l’opération de vérification 1 . La fidélité, paramètre de disper- 
sion, devra être explicitée plus finement. On retrouvera par exemple les 


ô 1. NF EN ISO 10012 septembre 2003, Systèmes de management de la mesure - Exi- 

gences pour les processus et les équipements de mesure. 
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notions de répétabilité et de reproductibilité (on remarquera une ana- 
logie assez forte avec la maîtrise statistique du produit). Cela va provo- 
quer un flou dans l’exactitude de l’expression de la mesure, appelé 
« incertitudes ». Nous consacrerons les chapitres 7 et 8 à l’étude des 
incertitudes avec comme objectif principal de vérifier l’adéquation 
« tolérances-incertitudes ». 

Il existe beaucoup d’autres paramètres pour caractériser les appareils de 
mesure et de contrôle. Pour plus d’information et surtout d’exactitude 
dans le choix des termes retenus, le lecteur est renvoyé au Vocabulaire 
international de la métrologie (VIM). 

Il faudra vérifier que ces qualités ne se dégradent pas au cours du temps. 
Il sera nécessaire de mettre en œuvre un suivi des appareils de mesure 
(chapitre 9). 


Étalonnage ou vérification des appareils 

L’étalonnage a comme objectif d’essayer d’exprimer une valeur vraie 
(jamais exactement connue). Imaginons par exemple une cale sans 
aucune inscription. En la comparant à des étalons situés en amont dans 
la chaîne d’étalonnage, on pourra graver la longueur estimée sur la cale. 
Pour un appareil qui donne la mesure directement (un pied à coulisse 
par exemple), on comparera le résultat fourni par l’appareil avec une 
valeur supposée vraie de l’étalon. On en déduira une correction (positive 
ou négative) à appliquer aux résultats fournis par l’appareil. Une 
démarche à partir de statistiques affinera l’estimation de la correction. 
Un étalonnage nécessite de relativement bien maîtriser les grandeurs 
d’influences perturbatrices (fig. 3.7). 
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Figure 3.7 : Exemple d'étalonnage 


50.02 
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La vérification ne cherche pas à garder une connaissance intrinsèque de 
l’appareil. Elle a comme objet de maintenir ou de rejeter l’appareil. 
Cette décision est prise par le vérificateur à l’aide de la position relative 
de la courbe de correction obtenue en mesurant différents « mesurandes 
étalons » par rapport aux limites d’écarts maximum tolérées (fig. 3.8). 
La vérification entraîne une perte d’information, qui devient binaire 
(appartenance ou non à la zone acceptable). Comme par la suite, on ne 
sait plus comment on est positionné à l’intérieur des limites d’erreurs 
tolérées, on associera souvent à cette position une distribution uniforme. 


La notion d’écart maximum toléré (EMT) représente un indice de per- 
formance (par exemple le 1/10 ou le 1/100 de mm). L’EMT est défini 
par rapport à l’objectif de qualité du produit que l’on doit surveiller 
(voir la notion de capabilité traitée au chapitre 7). L’appareil sera déclaré 
conforme ou non, on n’enregistrera pas les écarts constatés. 


Figure 3.8 : Exemple simplifié de courbe de correction 


Écart 


Limite supérieure (+ EMT) 



Mesurande 


Limite inférieure (- EMT) 
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Dans le cas d’étalonnages internes, il faut définir des procédures d’éta- 
lonnage (certaines procédures sont normalisées, il est préférable de les 
utiliser). L’approche ISO « Spécification géométrique des produits 
(GPS) » s’intéresse non seulement à la définition du produit mais égale- 
ment à sa vérification. Par exemple, pour le comparateur mécanique à 
cadran, on pourra consulter la norme NF EN ISO 463 :juin 2006. 

Cas de la sous-traitance (vérification externe) 

Cette activité de suivi, d’étalonnage ou de vérification de tout ou partie 
des appareils n’est pas toujours possible en interne dans l’entreprise. Elle 
demande beaucoup de temps et des moyens humains et techniques qui 
peuvent faire défaut. Il est possible de contacter un laboratoire de 
métrologie pour sous-traiter tout ou partie de ces tâches. Par exemple, si 
l’entreprise possède de nombreux appareils identiques, l’étalonnage de 
ceux-ci pourra se faire en interne, à l’aide de cales-étalon vérifiées en 
externe. L’organisme doit, soit s’assurer par des audits que le laboratoire, 
en position de fournisseur, présente la garantie métrologique nécessaire, 
soit que le laboratoire possède l’accréditation COFRAC correspondant 
aux besoins. Dans ce second cas, la traçabilité du raccordement aux éta- 
lons nationaux est alors garantie. 

Suivi du matériel 1 

Une dérive des moyens de mesure risque d’entraîner une dérive de la 
production. Pour éviter ce danger potentiel coûteux, il faut étalonner ou 
vérifier régulièrement les appareils de mesure. Reste le problème de 
savoir à quelle fréquence faire ces vérifications (chapitre 9). 


1 . Références normatives : 

- FD X07-018 décembre 1997, Métrologie - Métrologie dans l’entreprise - Fiche de 
vie des équipements de mesure, de contrôle et d’essai ; 

- NF EN ISO 10012 septembre 2003, Systèmes de management de la mesure - Exi- 
gences pour les processus et les équipements de mesure. 
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Gestion des moyens de mesure 


Déclaration de conformité 

À l’issue de chaque vérification, l’appareil peut être : 
s> maintenu pour l’utilisation prévue (conformité) ; 

* ajusté ou réparé et remis en service ; 

► déclassé (utilisation nécessitant une précision plus faible) ; 

► mis au rebut. 

Les deux derniers cas relèvent de la maîtrise de produits non conformes, 
ils doivent être clairement identifiés. 

Documents de suivi 

La vérification engendre un constat de vérification 1 précisant le 
maintien ou non de l’appareil dans le service concerné. L’étalon- 
nage donne un document d’étalonnage 2 caractérisant l’appareil. 

Une étiquette (ou plaque) de confirmation permet immédiatement à 
l’utilisateur ou à l’auditeur de connaître pour chaque appareil son code, 
sa validité et sa date limite d’utilisation (fig. 3.9). 

Figure 3.9 : Exemple d'étiquette de confirmation 


Affectation : 
Référence : .. 
Étalonné le : 


Date de validité : 


Cette étiquette doit être si possible solidaire de l’appareil (auto-adhésive 
ou à œillet ou marquage indélébile). 


1. X07-011 décembre 1994, Métrologie- Essais - Métrologie dans l’entreprise - Constat 
de vérification des moyens de mesure. 

2. FD X07-012 novembre 1995, Métrologie - Métrologie dans l’entreprise - Certificat 
d’étalonnage des moyens de mesure. 
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On utilise parfois un système simplifié utilisant des points de couleurs 
(peinture indélébile). Chaque couleur fait référence à un échéancier qui 
associe une année de fin de validation de l’appareil à une couleur spéci- 
fique. 

Une fiche de vie 1 , référencée par l’étiquette, permet de connaître exacte- 
ment l’historique de l’appareil (entretien, maintenance, réparation, éta- 
lonnage, déclassement, concepteur, fournisseur, etc.). Une gestion 
documentaire gère cette fiche en assurant le suivi et planifie automati- 
quement les intervalles de vérification ou d’étalonnage. Elle permet éga- 
lement d’accéder au dossier technique de l’appareil. 

Il existe de nombreux logiciels facilitant cette tâche de gestion. Ils assu- 
rent par exemple les missions suivantes : 

- état du stock (en service, en vérification, réformés. . .) ; 

- édition d’une fiche de vie d’un appareil ; 

- prochain ordre de vérification, liste des retards. . . ; 

- aide à la détermination des fréquences de vérification ; 

- etc. 

Cette gestion des moyens de mesure, avec les suivis et ordres de vérifica- 
tion, peut être complètement externalisée. Beaucoup de fournisseurs de 
logiciels proposent cette option. 

Formation du personnel 

Le meilleur des matériels n’est rien si l’on ne maîtrise pas son emploi. 
Les mesures deviennent souvent très complexes du fait des caractéristi- 
ques à vérifier (tolérances géométriques par exemple), des techniques 
employées (métrologie tridimensionnelle) et de la culture mathéma- 
tique associée (statistiques, analyse du signal. . .). 

Même avec une technologie simple à mettre en œuvre, il est toujours 
bon de prendre quelques instants pour s’assurer que l’opérateur a bien 
compris sa mission et qu’il possède suffisamment d’informations pour 
garantir le respect de l’auto-maîtrise sur son poste de travail. 


1. FD X 07-018 décembre 1997, Métrologie dans l’entreprise - Fiche de vie des équipe- 
ments de mesure, de contrôle et d’essai. 
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Chapitre 4 


Erreurs de mesurage 



Lors d’un mesurage, le résultat de la mesure va être perturbé par diffé- 
rentes grandeurs d’influence engendrant des erreurs. Nous ne cherche- 
rons pas, dans ce chapitre, à quantifier les incertitudes associées à ces 
erreurs, mais simplement à donner des méthodes et règles d’analyse 
pour le métrologue, en recherchant les causes et comment les mini- 
miser, voire les supprimer. 


Analyse des sources d'erreurs 

Le mesurage (ou processus de mesure) est sensible à certains paramètres. 
Une partie de ces paramètres est immédiatement identifiée par le métro- 
logue. L’autre partie est plus sournoise, au sens où le métrologue n’en a 
pas forcément conscience. Cela va augmenter la part d’aléatoire sur les 
paramètres recensés, par ignorance de l’effet des paramètres non pris en 
compte (fig. 4.1). 


Figure 4.1 : Exemple de la non-prise en compte de l'influence 
de la température 






e> 

© 
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Recherche des influences par la méthode des 5 M 1 

Les paramètres influents sont classés en cinq types : 

► influence de l’appareil et des composants associés ( machine ou moyen 
de mesure) ; 

► influence de l’opérateur (; main-d'œuvre ) ; 

► influence de la matière ( matière ou mesurande ) ; 

► influence de l’environnement ( milieu ) ; 

► influence de la procédure de contrôle {méthode). 

Rappel : le MSA 2 3 retient un découpage légèrement différent en distin- 
guant les étalons de l'appareil et en rattachant la méthode et l'opéra- 
teur. Il est mémorisé à l'aide de l'acronyme : SWIPE (Standard, Work 
Piece, Instrument, Person/Procedure, Environment). 

Il est parfois difficile de pouvoir partager les influences aussi nettement. 

En particulier, la méthode de mesure peut être étroitement liée (interac- 
tions) à : 

* la main-d’œuvre : en dehors de l’expertise physique, il faut tenir 
compte de la compétence intellectuelle qui permettra de bien appli- 
quer les méthodes retenues ; 

► le moyen de mesure : de plus en plus de logiciels embarqués, qui appli- 
quent des méthodes de calcul et d’optimisation, font partie de l’ins- 
trument (par exemple, choix du type de filtrage pour calculer un 
critère de rugosité) ; 

* la matière ou mesurande : par exemple, pour la prise d’un diamètre 
sur un axe, l’opérateur recherchera systématiquement le diamètre 
maximum dans le cas d’ovalisation ; 


1 . Duret D., Pillet M., Qualité en production. De l’ISO 9000 à Six Sigma, Paris, Éditions 
d’Organisation, 3 e éd., 2005. 

2. Measurment Systems Analysis : manuel développé sous la responsabilité de Daimler Chrysler 
Corporation, Ford Motor Company et General Motors Corporation ( ioww.carwin.co.uk/qs ). 

3. Breissand G., « Une vision industrielle de la métrologie », Contrôles, Essais & Mesures, 
janvier 2004 propose la répartition suivante : mesurande (10 %), instruments (5 %), 
opérateur (15 %), milieu (20 %) et méthode (50 %). 
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Erreurs de mesurage 


* le milieu : une méthode de mesure est au moins systématiquement 
perturbée par une grandeur d’influence particulière. Si celle-ci est 
faible, elle sera négligée sous réserve que l’on ne recherche pas une 
finesse de mesure trop élevée (par exemple, l’influence de la poussée 
d’Archimède sur la mesure de masse). 


Comment lutter contre les grandeurs d'influence 

On peut recenser quatre méthodes pour contrer l’influence des gran- 
deurs perturbatrices. On peut procéder par : 

élimination : c’est la méthode la plus radicale mais souvent la plus 
onéreuse (par exemple, un laboratoire climatisé). De plus dans 
certaines situations, cette élimination est impossible (mesure sur 
site industriel). Par contre lorsque l’on donne un résultat, il faut 
garder à l’esprit que ce résultat correspond de fait à des niveaux 
bien définis des grandeurs d’influence de l’environnement. Par 
exemple, lorsque l’on caractérise la longueur d’une pièce, cela 
suppose 1 2 * que : 

- la température est à 20 °C, 

- la pression atmosphérique est de 101325 Pa, 

- l’humidité relative est d’environ 55 %, 

- l’effort de mesurage est nul, 

- la pièce n’est pas contrainte par l’ablocage, 

- etc. ; 

insensibilisation : l’appareil de mesure est insensible à la grandeur 
d’influence. Par exemple pour un instrument en invar 2 à très faible 
coefficient de dilatation, les variations de température ne sont pas 
perturbatrices ; 


1. Voir NF E 10-100-1 décembre 2002, Spécification géométrique des produits - 
Méthode de mesurage dimensionnel. 

2. Invar (TM) est une marque déposée d’Imphy Alloys depuis 1907. C’est un acier à forte 

teneur en nickel (36 %). Son coefficient de dilatation est de l’ordre de 1 . 1 0 6 K 1 . 
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compensation : on cherche à combiner deux signaux de mesure de 
manière à ce que le signal résultant (une somme algébrique par 
exemple) ne soit pas perturbé. Cette combinaison permet de s’astreindre 
à l’influence de certaines grandeurs perturbatrices comme la tempé- 
rature (par exemple, balancier compensé en horlogerie) ou le posage 
(développée au chapitre 5) ; 

correction : c’est la solution la plus ambitieuse. Elle nécessite d’avoir 
modélisé l’influence étudiée. Il peut s’agir d’une correction mathé- 
matique ou physique (par exemple, came linéaire mise en parallèle 
d’une vis pour corriger les erreurs de pas en agissant sur l’indexation 
du vernier). 

Sans être exhaustifs, nous donnerons à titre d’exemples quelques appli- 
cations des méthodes précédentes appliquées à l’influence de l’instru- 
ment de mesure et à l’influence du milieu. L’influence globale du posage 
est très importante en métrologie dimensionnelle, nous l’étudierons en 
détail dans le chapitre 5. 


Influence du moyen de mesure 

Par moyen de mesure, nous considérons bien sûr l’appareil, mais égale- 
ment les matériels annexes comme un marbre, un support, une potence, 
etc. 

Les appareils de mesure paraissent lourds et massifs. Cela vient du fait 
qu’ils sont calculés en déformations et non en contraintes. Les liaisons 
mécaniques sont réalisées avec un jeu très faible. Il est possible d’uti- 
liser, lorsque l’on n’a besoin que de faibles courses, la déformation de 
lames métalliques (par exemple, accéléromètres, montages de compara- 
teurs...) (fig. 4.2). 


I 

e 


o 
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Erreurs de mesurage 


Figure 4.2 : Liaison glissière réalisée par déformation de lames métalliques 



Le concepteur d’appareil de mesure cherche autant que possible à sup- 
primer les erreurs liées à l’appareil par élimination, c’est-à-dire qu’il vise 
un appareil parfait. Pour cela, il doit se préoccuper : 

► des erreurs de lecture (par exemple, défaut de parallaxe) ; 

* des erreurs de mobilité (incidence du frottement, phénomène de 
stick slip) ; 

► de la géométrie des palpeurs (filtrage ou non des défauts de rugosité, 
par exemple) ; 

► de la vitesse d’accostage lors de palpages mécaniques sur une machine 
à mesurer tridimensionnelle. On constate souvent des différences 
significatives entre la première mesure faite manuellement et celles 
faites en automatique à vitesse normale (bien souvent de meilleure 
qualité) ; 

» de la qualité géométrique des mécanismes (jeux et dispositions rela- 
tives). Ernst Abbe a établi un principe de conception qui est toujours 
d’actualité ; 

* des efforts retenus au niveau du contact (capteur) avec la pièce ; 

» etc. 

Les dilatations dues à la température modifient la structure. On ren- 
contre parfois l’utilisation de bilames (juxtaposition de deux matériaux 
de coefficient de dilatation différent) qui vont permettre en se défor- 
mant, de conserver une longueur ou un moment d’inertie. Nous verrons 
d’autres utilisations aux niveaux des posages utilisant ce principe de 
compensation (fig. 4.3). 
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Figure 4.3 : L'inertie de ce volant est pratiquement insensible 
à l'augmentation de température 



Principe d'Abbe 1 

« Le système de mesure doit être placé dans le prolongement (coaxiale- 
ment) de la pièce à mesurer (mesurande). » Pour illustrer ce concept, 
nous retiendrons un pied à coulisse, appareil qui ne respecte pas le prin- 
cipe d’Abbe (fig. 4.4). 

Figure 4.4 : Incidence du jeu J sur la mesure pour le pied à coulisse 



1. Physicien allemand, Ernst Abbe (1840-1905) était professeur et directeur de recherche 
dans la société Cari Zeiss. 
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Dans cet appareil, l’étalon embarqué (les graduations) est en parallèle de 
la pièce à mesurer (mesurande). Il est situé à une distance h qui peut être 
importante par rapport au mesurande x. Au lieu d’avoir y = x, nous 
obtenons : 


y — x— h ■ tan(ct) « x — h OC 


x-h- 


l 

L 


Note : pour ce calcul, les angles étant petits, nous n'avons conservé que 
les infiniment petits du premier ordre ( sin(a) ~ tan(a) ~ct et 
cos(a) ~ 1 ). Cette technique de calcul sera souvent retenue en métro- 
logie dimensionnelle (parfois en passant au deuxième ordre : 

cos (a)~ 1— R emarc l uons Q ue I e défaut angulaire pourrait égale- 
ment être dû au défaut de rectitude de la glissière (non pris en compte 
dans cet exemple). 


Figure 4.5 : Type d'appareil respectant le principe d'Abbe 



Le micromètre est un appareil qui respecte ce principe (fig. 4.5). Ce 
type de conception engendre des appareils encombrants, ce qui n’est pas 
toujours possible. Par exemple, sur une machine à mesurer tridimen- 
sionnelle type portique, seul l’axe Z respecte ce principe. Les axes Xet Y 
= retiennent l’architecture vue pour le pied à coulisse. 

Q- 

(5 

© 
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Influence des efforts de contact 

La mesure d’une pièce suppose que celle-ci ne soit soumise à aucune 
force. Tout effort de contact va engendrer des contraintes générant des 
déformations globales au sein de la structure de la pièce et du moyen de 
mesure, et locales par pression de contact. Nous avons vu que dans une 
chaîne de mesure, une énergie d’activation était nécessaire pour assurer 
la prise d’information. Cet effort de contact est généré par l’opérateur 
(cas du pied à coulisse) ou par le moyen de mesure. Parfois, un dispositif 
sur l’instrument permet simplement de limiter l’effort fourni par l’opé- 
rateur (cas du limiteur de couple sur le micromètre). 

Pour les déformations de contact entre les solides Sj et S2, elles vont 
dépendre : 

d du type de contact (ponctuel, linéique ou surfacique) ; 

► des rayons de courbure des solides en contact (Rj et R 2 ) ; 

► de la nature des matériaux (module de Young E et coefficient de 
Poisson u). 

Cet écrasement peut-être modélisé à l’aide des formules de Hertz 1 

(%. 4.6). 

Figure 4.6 : Application des formules de Hertz pour un contact ponctuel 
(sphère/sphère ou sphère/plan) 

Rayon R; : 
sphère 
de palpage 


Rayon R2 : 
sphère 
de référence 


I 

e 


a 

1. Heinrich Rudolf Hertz (1857-1894), ingénieur et physicien allemand. 


Écrasement total : 
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Dans l’exemple précédent avec des billes en céramique (Ri = 1 mm) et 
en acier (R 2 =15 mm), on trouve un écrasement total d’environ 0,3 /toi 
pour une force de contact d’1 N. 5 sera le plus souvent négligé ou com- 
pensé (notion de rayon dynamique). Cette formule est également 
valable pour le contact de type sphère plan, il suffit de faire R 2 = 

Maîtrise des influences dues à l'appareil de mesure 

Le métrologue devra choisir un appareil approprié à son problème 
(course utile, résolution, pièce déformable ou non, maîtrise de l’effort de 
contact, etc.). S’il reste une influence non maîtrisable, on pourra : 

► estimer le caractère systématique de l’influence et vérifier qu’elle est 
inférieure à un seuil qui est négligeable par rapport aux autres sources 
de variation. En particulier, on peut vérifier que l’on reste dans la classe 
d’exactitude spécifiée par un indice de classe. Dans le cas contraire, il 
faut envisager de faire des corrections ; 

► réduire les effets aléatoires (par exemple, l’influence d’un jeu) en tra- 
vaillant à partir de moyenne. L’incertitude est réduite en V«. 

Figure 4.7 : Diminution de la part aléatoire de l'incertitude en fonction 
du nombre de mesures 


1 ■ 

0 , 8 - 
0 , 6 - 
0,4- 
0 , 2 - 
0- 

0 5 10 15 20 25 30 35 40 45 50 


Nb mesures 


Influence du milieu 

Tout mesurage se fait dans un environnement qui va plus ou moins 
perturber les résultats de mesure. Recenser exhaustivement toutes les 
grandeurs d’influences dues au milieu est impossible. On demandera 
simplement au métrologue d’être conscient des principales sources 
d’influences. 
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Exemples de grandeurs d'influence dues au milieu 

La température : c’est celle qui est le plus souvent citée spontané- 
ment. Attention en particulier aux gradients de température qui peu- 
vent être engendrés par l’éclairage et provoquer des déformations 
non homogènes (marbre par exemple). En micromécanique, son 
influence reste faible puisque la dilatation est proportionnelle à la 
longueur de base. Nous l’étudierons plus en détail au chapitre sui- 
vant. 

Les vibrations du sol et des supports (ainsi que les chocs) : 

influence difficile à modéliser, on cherchera à l’éliminer par isolation 
(dalle flottante, support filtrant, suspension pneumatique, etc.). 
L’hygrométrie : le taux d’humidité, en dehors de la création néfaste 
d’oxydation, peut perturber la géométrie de certains composants (les 
marbres par exemple). Un air trop sec (gênant pour l’opérateur) peut 
également générer des dysfonctionnements électroniques. 

Les champs électriques et magnétiques : les champs magnétiques 
peuvent être provoqués par l’environnement du processus de pro- 
duction ou plus simplement par le processus lui-même (ablocage, 
frottement de la meule en rectification). La pièce peut se déformer 
sous sa propre aimantation résiduelle. Le CETIM 1 montre qu’une 
opération supplémentaire de désaimantation dans la gamme permet 
de minimiser le défaut de circularité. On peut noter également 
comme influence électrique, les variations de tension du réseau d’ali- 
mentation. 

La latitude et l’altitude : le mesurage des masses est souvent une 
mesure indirecte, où l’on détermine le poids et ensuite la masse en 
introduisant la gravité g. La gravité 2 n’est pas une constante mathé- 
matique, elle dépend de la situation géographique et de la nature des 
sols. 

La pesanteur : la gravité a également une influence sur la déforma- 
tion de la machine et de la pièce (par exemple, une pièce longue 
n’aura pas la même longueur en position horizontale et verticale). 


1 . Grandadam A., La Métrologie dimensionnelle en mécanique industrielle, CETIM, 1 977. 

2. Lorsque l’on veut prendre en compte son influence, il faut tenir compte principalement 
de l’altitude et de la latitude {http://www.ptb.de/en/org/lH H115l_index.htm). 
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Tout objet se déforme sous son propre poids (rappelons que les appa- 
reils de métrologie sont souvent massifs). Il sera important de bien 
situer les appuis pour minimiser les déformations résiduelles (points 
supports 1 pour les marbres : points de Bessel 2 situés à 0,22 de la lon- 
gueur ou de la largeur par rapport aux bords). 

La propreté : à la fois du milieu ambiant (poussière, brouillard 
d’huile) et de la pièce (trace d’huile par exemple), surtout dans le cas 
de mesures optiques. 

Les éclairages latéraux : par exemple dans le cas de mesure par 
vision, une baie vitrée située à proximité peut être la source de per- 
turbations (type : coucher de soleil). 

Le temps : action à long terme, il provoquera des dérives lentes (par 
exemple, fluage de certains composants) et à court terme (décharge 
de condensateurs dans des mesures de rugosités ou de forces, base de 
temps. . .). Il est recommandé un temps de stabilisation avant utilisa- 
tion. 

» Etc. 

Classer les grandeurs d'influence liées au milieu 

Les influences étant nombreuses, on demandera au métrologue d’avoir 
une démarche de physicien. Cela permettra de réellement retenir les 
influences prépondérantes pour pouvoir éventuellement faire des cor- 
rections. Les autres seront négligées, elles augmenteront le bruit de fond 
de mesure (incertitude résiduelle). Pour illustrer cette démarche de clas- 
sement, nous allons l’utiliser lors d’une comparaison de masse entre un 
étalon de 1 kg en platine (étalon de référence) et un étalon supposé de 
1 kg en acier inoxydable (par simplification, les éléments géométriques 
de la balance seront supposés parfaits) (fig. 4.8). 


© 


1. NF E 10-100-1 décembre 2002, Méthodes de mesurage dimensionnel. 

2. A. Friedrich Bessel (1784-1846), astronome et mathématicien allemand. 
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Figure 4.8 : Principe de comparaison de deux étalons de matières différentes 



Dans l’exemple précédent on croit comparer deux masses, et dans ce cas 
avoir : m 1 — m p . 

En réalité on compare des couples, qui peuvent se ramener à la compa- 
raison de forces avec l’hypothèse de bras de levier identiques. Dans le 
vide, on peut garder comme forces uniquement les poids, alors que dans 
l’air il faut en toute rigueur tenir compte de la « poussée d’Archimède ». 
L’équation devient : 

m r g — Vr P Air ' g - m p' g -Vp' P Air ' g 
m,-V r p Atr — m p — V p ■ p Air (indépendant de g) 


m, m p 

m l- — -p A ir= m p- — -P A ir 

Pi Pp 

I P Air 



Pi 


Avec les valeurs approximatives suivantes : 

► Pp ~ 22 000 kg/m 3 (platine à 10 % d’iridium) ; 

* Pl~8 000 kg/m 3 ; 

» PAir ~ 1>2 kg/m 3 (à 20 °C, 50 % d’humidité relative et faible alti- 
tude) 1 . 

k est légèrement inférieur à 1,0001. Ce qui veut dire que l’erreur systé- 
matique commise est d’environ 0,1 gramme. 

1 . Voir le site www. metgen. org! pour calculer p Atr avec plus de précision. 
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Le plus souvent cette erreur sera négligée. Avec des étalons de niasses 
volumiques voisines ; les poussées d’Archimède sont pratiquement 
identiques, elles se compensent. 

Étude de l'influence de la température 

Pour cette grandeur particulière, nous allons reprendre les quatre méthodes 
déjà citées, à savoir : 
g> élimination ; 

► insensibilisation ; 

► compensation ; 

► correction. 

Élimination 

Rappelons que les longueurs d’un objet sont implicitement définies à la 
température de 20 °C 1 . Lors de métrologie fine (recherche du /2 m), il 
est naturel de chercher à faire des mesures à cette température. Cela 
nécessite d’avoir un local climatisé, avec un sas d’entrée pour éviter les 
perturbations brusques, de maîtriser les flux d’entrée des personnes 
(ordre de grandeur : 60 à 80 watts par personne), de choisir un emplace- 
ment à l’abri des changements climatiques (par exemple, éviter une face 
sud avec des baies vitrées), etc. Toutes ces précautions ne suffisent pas. 
Une pièce venant du lieu de production a probablement une tempéra- 
ture interne différente. Elle ne passera pas instantanément à la tempéra- 
ture de 20 °C. Le temps de transition pour atteindre l’équilibre 
thermique est complexe à modéliser. Empiriquement, on peut constater 
qu’il sera proportionnel (pour une même matière) au rapport du 
volume par la surface de l’objet (peau de la pièce). Notons que ce rap- 
port n’est pas sans dimension, il faut travailler dans les mêmes unités 
pour faire des comparaisons (fig. 4.9). 


e> 

© 


1. NF E 10-100-1 décembre 2002, Méthodes de mesurage dimensionnel. 
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Figure 4.9 : Exemple de quelques rapports V/S 



Z 


7 


Sphère : 0 D 


Cube : arête L 


Plaque Lx Lx h 


V_D 
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Pour une même morphologie (par exemple, une pièce cubique), ce rap- 
port est fortement influencé par la taille : 


Le cube de 10 mm d’arête se refroidira environ cent fois plus vite. Le 
CETIM 1 donne pour un tampon cylindrique de diamètre 60 mm, un 
temps de 6 heures pour passer de 25 °C à 20 °C. 

Nous verrons qu’il est possible de faire des corrections pour prendre en 
compte une température différente de 20 °C. Ces corrections ne seront 
efficaces que si toute la matière de la pièce est à la même température. 
En cas de gradient de température à l’intérieur de la pièce, cela va engen- 
drer des déformations, particulièrement dans le cas d’architecture com- 
plexe (par exemple des nervures, rayons, etc.). 

Insensibilisation 

L’idée est de retenir des matériaux peu sensibles aux variations de 
températures. Cette aptitude est quantifiée par un indicateur appelé 
« coefficient de dilatation » exprimé en K - 1 (tableau 4.1). 


1 . Grandadam A., La Métrologie dimensionnelle en mécanique industrielle, CETIM, 1 977. 


cube de 1 m d’arête : — = 1 ; 

5 6 



5 6 


O 
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Tableau 4.1 : Exemples de coefficients de dilatation 


Matériaux 

Coefficient (10 _6 /K) 

Vitro-céramique 

0.02 

Invar 

0,6 à 1 ,2 

Céramiques AI203 

7 

Céramiques SiC 

4 

Céramiques Si3N4 

3.5 

Céramiques Zr02 

2,3 à 3,2 

Carbure de tungstène 

4,5 6 7 

Diabase 

5 

Verres 

669 

Granit 

7,5 6 8 

Alliages titanes 

8,5 6 9,5 

Aciers 

11,5 

Inox ferritiques & martensitiques 

10.5 

Inox austénitiques 

156 16 

Fontes grises 

106 11 

Fontes grises GS 

12,5 

Fontes blanches 

20 

Bronzes 

17,5 6 19 

Laitons 

1 8(CuZn5) à 21 (CuZn40) 

Alliages légers 5356 

22,5 

Alliages légers 2017A 

24 

Alliages de Zinc 

26 

Magnésium 

26 


L’insensibilisation du moyen de mesure à la température ne sera intéres- 
sante que si tous les composants de l’acte de mesurage sont également 
concernés. Par exemple, si l’on compare une pièce en acier (A^ ae , r 
= 11,5 • 10" 6 K -1 ) avec une cale-étalon céramique en Zr02 (A céramique 
= 2,8 • 10“ 6 K" 1 ) à une température différente de 20 °C, la différence 
de dilatation sera perturbatrice. 
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Compensation 

La méthode par compensation ne cherche pas à annuler l’influence de la 
température (dilatation), mais à essayer de contourner ses différents 
effets (au niveau de chaque composant), de manière à ce que le résultat 
global ne soit pas perturbé. Citons, par exemple, le problème historique 
de mesure de la longitude. John Harrison 1 pour améliorer la qualité de 
son horloge marine chercha à compenser l’influence de la température 
sur la période du pendule (maîtrise de la longueur) en combinant judi- 
cieusement des matériaux de coefficients de dilatation très différents 
(1725) (fig. 4.10). 

Figure 4.10 : Principe de la compensation de longueur 
sur un balancier 



Lors de la mesure de cales de haute qualité, on emploie souvent une 
méthode par comparaison en utilisant des cales de références de même 
matière que les cales à vérifier (fig. 4. 1 1). 


1. John Harrison (1693-1776), modeste horloger/menuisier, eut beaucoup de mal à faire 
reconnaître à la communauté scientifique anglaise le bien-fondé de son approche de la 
mesure de la longitude grâce à une horloge de haute précision. Ses différentes inventions 
firent progresser l’horlogerie et eurent un impact financier important sur la fiabilité du 
trafic maritime en permettant d’éviter de nombreux naufrages. 
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Figure 4.1 1 : Comparaison de cales de même matériau 



Correction 

Faire une correction suppose que l’on est capable de prévoir les consé- 
quences de la grandeur d’influence. Dans le cas d’une influence variable, 
non prévisible a priori comme pour la température, seule la correction 
de type modélisation du phénomène est envisageable. 

Pour des pièces métalliques (dilatation thermique isotrope), l’allonge- 
ment peut être modélisé par une fonction linéaire sous réserve qu’il n’y 
ait pas de changement de phase (par exemple, passer d’une structure 
« cubique centrée » à « cubique à face centrée ») : 

Le - L 2 0° ' [l + A Matière 20°)] 

À partir de la mesure faite à la température 0, on peut en déduire la lon- 
gueur de la pièce à 20 °C (la mesure pouvant être faite également en kel- 
vins puisqu’il s’agit d’une différence de température) : 

^20° ~ Lq ~ L 2 o° ‘ ^ Matière 20 °) 

L 2 o° ~ Le - 20° 

I, avec : C e ^ 20 „ = L 20 „ ■ ^ Matière ■ (ô - 20°) ~ L e ■ A Matjire ■ (d - 20°) 

Q_ 

ô 
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Exemple de correction de température 

Nous comparons dans un atelier à la température de 26 °C, une pièce en 
acier avec une cale-étalon en céramique (Zr02) de valeur 75,000 mm. 
L’écart e donné par le comparateur est de + 5 /im. Quelle est la longueur 
de la pièce à 20 °C (fïg. 4.12) ? 

Figure 4.1 2 : Exemple d'illustration de correction de température 



Retenons les valeurs de coefficients de dilatation suivants : 
hpièce= 1 1.5 • 10- 6 K- 1 et X Étalon = 2,8 • 10- 6 K- 1 
Si la pièce et l’étalon sont dans l’atelier depuis suffisamment longtemps, 
on peut considérer qu’ils sont à la même température 6 : 


soit : 

Ap,20° + -(0-20°)] = L E 20 o -[l+ A jj ■(0-2O°)] + £ 

Rappelons que, par définition, la valeur gravée sur la cale correspond 

^ L e,2 o° 

_ L E , 20 ° ' [ l + ^E • { e - 20 °)] + e 

p ’™° = [l + V(0-2O°)] 
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soit la valeur approchée : 

\ + X E -(0 — 20 °)] 

’l + V(0-2O°)] 

Valeur numérique : 

► sans correction : L p 20 o = 75,005mm ; 

* avec correction : L p 20 o = 75,00\mm. 

Remarques : 

• Le lecteur peut remarquer que lorsque l'étalon et la pièce ont le même 
coefficient de dilatation, la correction est inutile. 

• La valeur corrigée n'est qu'une estimation de la valeur vraie, l'intervalle 
de confiance associé à cette valeur (incertitudes) dépend de nombreux 
paramètres (en particulier des coefficients qui ne sont pas connus 
exactement). Nous reviendrons dans le chapitre 8 sur ce point. 



o 
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Chapitre 5 


Maîtriser l'influence 
du posage 



Par idéal, lors de métrologie dimensionnelle, le posage de la pièce ne 
devrait pas avoir d’influence perturbatrice. Pourtant bien souvent, il y 
a un risque important d’imputer à la mesure d’une caractéristique de 
la pièce un complément d’écart provenant en réalité d’un mauvais 
posage. 

Comme pour les autres grandeurs d’influence, nous chercherons à inter- 
venir par : 

élimination : plus modestement, nous chercherons à garantir un 
posage de qualité garantissant si possible l’unicité dans le temps des 
differentes mises en position ; 

insensibilisation : retenir un processus de mesure qui soit insensible 
aux variations de posage ; 

compensation : combiner des signaux de mesure corrélés pour 
s’affranchir des erreurs de mise en position ; 

correction : nous essayerons de modéliser l’erreur de mise en posi- 
tion (à partir de données initiales) de manière à pouvoir faire une 
correction pour éliminer l’influence des variations de posage. 

Ces différentes méthodes seront illustrées à l’aide d’exemples empruntés 
à la production mécanique. La lecture de ce chapitre n’est pas indispen- 
sable pour la suite de l’ouvrage. 


Élimination : améliorer l'isostatisme 
du montage de contrôle 

Si l’on prend une pièce et qu’on la monte plusieurs fois de suite sur le 
montage de contrôle, il faut que la pièce soit toujours positionnée de la 
même manière (en translation et en orientation). Si ce n’est pas le cas, 
cela peut provenir principalement : 

► de la déformation de la pièce ; 

► du frottement qui interrompt le glissement pour venir en appui ; 
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c du système d’ablocage engendrant des efforts qui vont faire bouger la 
pièce lors du serrage ; 

► et du choix « judicieux » des surfaces (ou points) d’appui réalisant 
l’interface « pièce-montage ». Plus particulièrement, le nombre de 
points d’appui ne doit pas être surabondant. 

Notions de degrés de liberté 

Un solide dans l’espace (non contraint) peut se déplacer en translation 
(trois composantes de déplacement) et en rotation (trois composantes 
angulaires). On dira qu’il a six degrés de liberté. 

Dans notre cas particulier, il s’agit de caractériser un ensemble de posi- 
tions dans un voisinage d’une position initiale cible. Comme il s’agit de 
micro-déplacements, la transformation géométrique permettant de 
passer de la position cible à la position réelle (ou réciproquement) peut 
être modélisée par un torseur (linéarisation) appelé « torseur de petits 
déplacements 1 » (les composantes sont considérées comme des infini- 
ment petits du premier ordre) : 


a 



Mesurer la position, revient à connaître les six grandeurs scalaires corres- 
pondantes. Remarquons que le déplacement (en mm) est modélisé par 
un champ de moments, le lieu où il est exprimé doit être clairement 
indiqué : 


Dp. ~ D 0 + &S/R A DPi 


1. Bourdet P. et Clément A., « Controlling a Complex Surface with a 3 Axis Measuring 
Machine », Annah ofCIRP, vol. 25/1/1976, p. 354-361. 
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ou sous forme matricielle : 


u p. 


u 0 


" 0 

-Y 

P' 


x p 

V P 

= 

v o 

+ 

y 

0 

-a 

• 

y p. 

w p i 


_ W 0_ 


_-P 

a 

0 


z p. 


La qualité d’un posage de contrôle sera caractérisée par la mesure de ces 
six grandeurs scalaires. Cette mesure « multicritères » va être fortement 
influencée par le choix des points d’appui réalisant « l’interface » entre la 
pièce et le montage de contrôle. 


Optimisation de la position des points d'appui 

L’objectif est de créer une liaison mécanique complète entre la pièce et le 
montage de contrôle. Le plus souvent, celle-ci sera réellement réalisée 
par six « petites surfaces de contact » matérialisant les six contacts ponc- 
tuels. Dans le cas de surfaces de contact différentes (plan sur plan, par 
exemple), on pourra toujours imaginer le nombre de liaisons ponc- 
tuelles équivalentes (exemple, trois points pour un appui plan). 

Exemple : mise en position de Lord Kelvin 1 (fig. 5.1) 


Figure 5.1 : Mise en position de Lord Kelvin et points de Bessel 



1. William Thomson (1824-1907), physicien britannique mieux connu en tant que Lord 
Kelvin. 
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Ce modèle de mise en position (connu également sous l’appellation : 
plan-trait-point) est très répandu. Nous l’utiliserons comme exemple 
conducteur. Si le choix paraît relativement évident (bien que l’appui 
plan puisse être construit différemment), il n’en est pas toujours de 
même lors de pièces à la morphologie tourmentée (pièces de forge ou de 
fonderie par exemple). 

Note : nous posons comme hypothèse que la pièce se comporte 
comme un solide indéformable. Dans le cas contraire, la disposition 
des points cherchera à lutter contre l'influence de la gravité (déforma- 
tions). Par exemple on pourra retenir les points de BesseP pour mini- 
miser la déformation due à l'action de la gravité (cas des marbres 
notamment). 

Coordonnées plückériennes d'une droite 

Partant d’un point de contact entre la pièce et le montage et après avoir 
défini le plan tangent nous pouvons construire la droite normale à ce 
plan et passant par le point de contact (fig. 5.2). 

Figure 5.2 : Droite normale (D ; ) au point de contact M, 



I 

e 


1. Les points de Bessel sont retenus pour minimiser les déformations de la longueur à g" 

mesurer en positionnant les appuis par rapport aux bords à respectivement 0,22.L et o 

0,22.1 (NFE 10-100-1). ® 
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Si nous associons un repère R(0,x,y,z) à l’ensemble « machine-mon- 
tage de contrôle », nous pouvons définir chaque droite normale par ses 
coordonnées plückériennes 1 . 

Figure 5.3 : Position spatiale de la droite (D ; ) dans le repère R 



n . est un vecteur unitaire de (Dj ) . Partant des coordonnées de : 
n xl Xj 

n. n yi et OM i yi 

n zi z i 

on calcule le produit vectoriel : g Q = OMi a n i 

Les coordonnées plückériennes sont définies par les six grandeurs (com- 
posantes) scalaires suivantes : 




x i- n yi-y r n > 


Julius Plücker (1801-1868), mathématicien et physicien allemand. 
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Ces six composantes ne sont pas indépendantes, elles sont liées par : 

2 , 2,2 -, 

n xi + n yi + 71 zi = 1 


{*** ~ Z i- n p) + n yi i Z i- n xi ~ x r»*)+- { X i- n yi-yi- n J) = 0 

Note : le point Mi peut être pris n'importe où sur la droite (propriété 
d'équiprojectivité). 

Rang d'un système de droite 

Les six droites normales (D ^ ) à (Dg) constituent un système de droites. 
On appelle rang du système, l’ordre le plus grand des déterminants non 
nuis que l’on peut extraire de la matrice construite à partir des coordon- 
nées plückériennes des droites appartenant au système. 



U x2 

n x3 

n xA 

n x5 

n x6 

"yl 

n y2 

n y3 

U yA 

n y5 

n y6 

U z\ 

n z2 

n z3 

n zA 

n z3 

n z6 

y v n zl -z v n yl 

So.x 2 

£o,x 3 

êo,xA 

£o,x 5 

y6- n z6~ Z 6- n y6 

Z V n x\ ~ X V n z\ 

&0,y2 

&0,y3 

&0,yA 

&0,y5 

z 6- n x6~ x 6’ n z6 

X V n y\~ y Y n x\ 

So,z2 

So,z3 

So,zA 

£o,z5 

X 6' n y6-y6- n x6 


Remarque : le rang d'un système est au maximum de 6. 

Théorème de Hunt 

Théorème : soit r le rang du système S des droites de contact (nor- 
males), les degrés de liberté restant entre les deux solides sont au nombre 
de d - 6-r 

Application : dans ce cas, nous cherchons à immobiliser (sans sur- 
contraindre) la pièce par rapport au montage. L’objectif est donc 
d’obtenir r= 6. 


92 


Maîtriser l'influence du posage 


Optimisation de la localisation des points d'appui 

Nous ne rappellerons pas les règles de bon sens pour positionner judi- 
cieusement les six points d’appui (par exemple, éviter trois points alignés). 
L’optimisation est multi-objectifs. Il faut minimiser simultanément 
six sorties (trois translations et trois rotations) en choisissant six points 
d’appui (les normales associées dépendent des interfaces « pièce- 
montage »). Nous proposons de remplacer l’optimisation multi-objec- 
tifs par un indicateur mono-objectif simple à calculer. Nous cherche- 
rons à maximiser le déterminant extrait de la matrice des coordonnées 
plückériennes. 

Exemple 

Reprenons le cas précédent en observant l'influence du coefficient k 
(pour les points 1 à 5 de la figure 5.4) sur le déterminant de la matrice 
construite avec les coordonnées plückériennes (rappel : points de 
Bessel pour k= 0,22). 

Figure 5.4 : Influence du coefficient k 
sur le déterminant des six normales 



J Les points d’appuis latéraux étant à mi-hauteur, cela donne en fonction 
de k pour L = 100 mm ; /= 50 mm ; h = 20 mm : 


© 
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Figure 5.5 : Le déterminant comme indicateur de qualité du montage 



On remarque que pour k= 0,5 nous n’avons que trois appuis efficaces, 
ce qui revient à une liaison rotule (trois degrés de liberté). 

Cette technique 1 permet de comparer différentes solutions, sous réserve 
d’un encombrement voisin (voir figure 5.6 pour la liaison de Boys qui 
donne de meilleurs résultats). 

Figure 5.6 : Exemple d'une liaison de Boys 



1 . La mesure globale de position fournie est fortement bruitée par le couple matière « pièce- g" 

montage » (déformations locales, frottement) et par le système d’ablocage. Cette propo- o 

sition de mesure globale est en cours de test (laboratoire SYMME). 
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Il est également possible de définir une mesure de la qualité de la liaison 
en utilisant une norme de la matrice des coordonnées plückériennes 
(P. Bourdet et A. Clément). Les normes matricielles associées à la 
matrice n’étant pas toujours aisées à calculer, on retiendra la pseudo- 
norme euclidienne : 

l k ™ 

Les éléments de cette matrice n’étant pas homogènes, on ne peut 
qu’optimiser un montage dont le type est déjà retenu. 


Insensibilisation 

Le processus de contrôle est tel que pour une position quelconque (en 
restant dans un voisinage acceptable), il donne la même valeur de sortie. 
Pour illustrer ce concept, nous retiendrons la mesure d’une pièce type 
« cale » par mesure différentielle avec une cale-étalon. 

Figure 5.7 : Exemple de mesure sans dégauchissage 



Si l’on ne prend pas garde au dégauchissage de la pièce (orientation de la 
pièce sur la machine suivant deux axes), le résultat sera erroné (fig. 5.7). 
Par exemple pour une pièce de 75,000 mm, inclinée de a = 2°, la 
mesure donnera : 


L - —^— = 75,048 mm 
cos et 
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Pour annuler cet effet, on envisage une mesure par comparaison, en uti- 
lisant la procédure suivante (fig. 5.8) : 

► mettre la pièce et la cale-étalon côte à côte (placage) ; 

► poser l’ensemble sur la machine ; 

► faire le zéro sur la cale-étalon ; 

» déplacer l’ensemble à l’aide du chariot transversal ; 

► mesurer l’écart entre la pièce et la cale-étalon. 

Figure 5.8 : Exemple de posage avec défaut d'orientation 
(vue de dessus) 



La longueur observée de l’étalon ne correspond pas à la valeur certifiée : 

LpûceObs ~ L EtdonObs + 6 
LpifreVraie _ ^EtalonVraie j , 

co s a co s a 

LpïèceVmie ~ ^ EtalonVraie + e ' 

Si CL est petit, le processus de mesure est pratiquement insensible à la 
perturbation d’orientation. 



I 

e 


o 
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Compensation 

L’objectif est de trouver une technologie de montage où les sorties 
locales (mesures) vont pouvoir se compenser. Intuitivement, c’est ce que 
nous pouvons vérifier sur la figure 5.9 pour la mesure d’un diamètre 
avec deux jauges inductives A et B (type stop cote). 

Figure 5.9 : Somme algébrique des sorties 
des capteurs inductifs 



Il est possible d’obtenir le même résultat avec un seul comparateur 
monté sur un cé flottant. 

Étude de l'incidence du déplacement 
de l'étalon dans la mesure d'un rayon 
(exemple avec trois capteurs) 

Présentation 

On veut mesurer le rayon d’un cercle à l’aide de trois capteurs linéaires 
inductifs placés à 120 degrés (mouvement plan sur plan, points de pal- 
1 page P-), avec ^ la sortie du capteur n° i. 

Q_ 

o 
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Figure 5.10 : Position des capteurs 



Imaginons l’étalon bien « centré », nous avons alors ^ = 0. Si 

nous déplaçons l’étalon, pour ne pas perturber le résultat, nous devrions 
avoir : + £ 3 = 0 


Position d'un capteur 

Chaque capteur est positionné par un point de palpage P i : 
et orienté par une normale n i n yi 


0 


0 

La droite support du capteur peut être représentée par ses coordonnées 
plückériennes : 




0 

0 

0 

X i' n yi~yi M xi 


! 

e 


o 
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Ces coordonnées sont liées par : ||«. || = 1 , soit : n xi + — 1. 

Caractérisation du déplacement de solide 

Le déplacement du solide étant petit, il sera défini par un torseur de 
petits déplacements (linéarisation) : 



Mesurer le déplacement revient à connaître les trois grandeurs scalaires 

»0« v 0 et Y 

Mesure locale au niveau d'un capteur 

Le capteur mesure la projection (£ z ) du déplacement du point Pj sur la 
normale n . (au premier ordre) : = D p •n i 

soit : 

*=/ ~{^ > o + ^ sir A ) * n i 

^ = (^.^) + (Ô^.Q^a^) 

Le produit scalaire D pih • n i peut se représenter comme le comoment de 
deux torseurs, le torseur de petits déplacements déjà étudié {S} et le tor- 
seur représentant les coordonnées plückériennes {9?} (ce point sera géné- 
ralisé au paragraphe traitant l’approche avec correction p. 97 et suiv.). 

£. = {2)}*{g>} 

= { u o- n x i + v o- n yi + Y-( x r n yi -Ji-nj) 
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Incidence du déplacement {2)} sur les capteurs 

On obtient un lien entre les trois composantes uq, Vq et y et les sorties 
capteurs et £3 en résolvant le système : 

H x\ n y\ X \ n y\ - y \ n x\ U 0 

n x 2 n y2 X 2 n y2-J2 n x 2 * v 0 
n x3 n y3 X 3 n y3-y3 n x3\ l Y . 


k' 
= ^ 
U> 


Exemple 

En reprenant la configuration de la figure 5.8, posons : 



Nous obtenons : 



Résolution avec choix de l'origine au centre du cercle 


2 


r 2 


I 

e 


o 
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À partir des différentes valeurs Ç., nous obtenons le système suivant : 


0*u 0 -l*v o + 0*y = Ç 1 



Nous pouvons calculer ^ + <ij 2 + en fonction de uq, vq et y. Nous 
avons bien : 

11 + £ 2 + ^3 “ 0 

Conclusion : le système est insensible à la position de l’étalon (sous 
l’hypothèse de petits déplacements). 


Correction 

Nous avons montré que pour la température, il était possible de faire 
une mesure à une température quelconque et ensuite prédire, grâce à 
une modélisation adéquate, ce qu’il en était réellement à 20 °C. Pour 
cela, il était nécessaire de connaître l’écart de température. Nous retrou- 
verons ici la même démarche. L’objectif est d’identifier l’erreur de 
posage et ensuite d’être capable, par logiciel (déplacement virtuel) ou 
par déplacement physique (table motorisée par exemple), de corriger 
cette influence perturbatrice. 

Il existe plusieurs approches pour proposer des corrections. Ce type 
d’étude pouvant donner lieu à un ouvrage complet, nous nous 
limiterons : 

► à l’association d’un élément géométrique à un nuage de points 
(méthode d’identification d’un élément optimisé) ; 
æ * à l’analyse modale d’une surface (mise en évidence des défauts de 
1 position et de forme). 

Q_ 

o 

© 
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Association d'un élément géométrique 
à un ensemble de points 

Identification de l'élément réalisé (fig. 5.11) 

Figure 5.1 1 : Palpage des points dans le repère machine 



L'élément est identifié 
par palpage. Sa connaissance 
se résume en un ensemble 
de coordonnées de points 
de l'élément réel dans un repère 
mesure entièrement localisé 
par rapport à la machine 


Création d'un repère local de dégauchissage (balançage) 

Ce repère peut être construit soit physiquement en étant lié à la pièce, à 
l’aide d’un nombre de points minimum (choisis en fonction du type de 
surface), soit géométriquement en s’appuyant sur la définition nominale 
de la pièce (modèle CAO). Par exemple, pour la première solution, dans 
le cas du plan, le choix de trois points, judicieusement choisis, sera suffi- 
sant pour définir ce repère local de calcul Rj^Op, xp, y p, zp). 

Figure 5.12 : Définition d'un repère brut de dégauchissage 



Les points P ; sont tous 
sensiblement à la 
même altitude z P = 0 

Notation : 

lO M P] = [X,Y,Am 

[ OpP J = [x,y,z] RP 


I 

e 


o 
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Remarquons que : 

► le choix des axes va permettre de définir un sens d’accostage lors du 
palpage d’un point. Par exemple si la surface réalisée est du type plan, 
l’axe zp correspond à la direction de palpage ; 

► pour définir le repère Rp nous avons retenu trois points assez éloignés. 
Ces trois points sont par définition d’altitude zp = 0. D’autres tech- 
niques pourraient être retenues. 

Mesure de l'élément réel 

En chaque point, nous allons mesurer (par calcul) la différence constatée 
entre la surface réelle (contact physique) et la surface nominale. Cette 
surface nominale est supposée reconstruite géométriquement dans le 
repère Rp. L’écart mesuré au point P; sera nommé (fig. 5.13). 

Figure 5.13 : Position de l'élément réel par rapport 
à l'élément nominal 


zp 



yp 



Nous avons ainsi un ensemble de mesures {£ z } qui caractérise la surface 
réelle. 


Exemple 

„ Dans le cas d'une planéité, le plan reconstruit correspond au plan (Op, 

=| x P , y P ). Le défaut de planéité brut est égal à : 

^ A Planéité = Sup(§ ) - lnf ( ) 

ô 
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Optimisation du dégauchissage 

Pour la définition du repère Æples trois points choisis ne conduisent pas 
forcément à la meilleure association, il reste une part d’arbitraire. Un 
autre choix aurait conduit à un repère Rp> voisin de Rp 

Figure 5.1 4 : Choix du repère de dégauchissage 



En utilisant la même technique de mesure, nous obtiendrions l’ensemble 
de valeurs (mesures) {<f 2 j. La technique proposée par Bourdet et Clé- 
ment 1 consiste à représenter la transformation géométrique de passage 
du repère Rp à Rpàe manière simple à l’aide d’un torseur (linéarisation) 
appelé torseur de petits déplacements. (Les composantes sont considé- 
rées comme des infiniment petits du premier ordre.) 


a 



1. Bourdet P. et Clément A., « Controlling a Complex Surface with a 3 Axis Measuring 
Machine », Annah ofCIRP, vol. 25/1/1976, p. 354-361. 
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Imaginons qu’après avoir fait la mesure, où pour le point Pj nous avons 
mesuré la valeur nous recommencions la mesure avec un nouveau 
repère de dégauchissage (fig. 5.15). 

Figure 5.15 : Déplacement de l'élément nominal 



Le plan contenant s’est déplacé. Le déplacement local est représenté 
par le vecteur D Pith> la nouvelle mesure vaut Cj. (En toute rigueur, la 
mesure ne s’est plus faite suivant mais suivant «/, l’erreur commise 
est du second ordre.) 

L’ensemble de mesures {^j est remplacé, suite au léger déplacement, par 
l’ensemble {e t } où chaque valeur e - est estimée par : 

e i=Zi~ D M %n i 


Par commodité de calcul, on peut prendre le point ou le point 
réel palpé P ( -, cela ne modifie pas le produit scalaire : 

J VÜ)'", 

I D pith %n i= D pi %n i 

© avec : p.P ith —k» n i 
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Le produit scalaire D Pith * n i peut se représenter comme le comoment de 
deux torseurs : 

► le torseur de petits déplacements {S} ; 

» et un torseur représentant le point de palpage et la direction de pal- 
page associée. Cela revient à situer une droite dans le repère associé. 
Ce torseur {SP fi appelé torseur des coordonnées plückériennes, est défini 
comme suit : 


M- 


°p 

Vérifions que : {2)}*{9L} = D Pith • n { 

Nous obtenons en développant le comoment (attention de bien 
exprimer les torseurs au même point) : 

Le champ de moment permet d’écrire : 

%, ~ D Pith + ^ A P lth O P 

Ce qui donne : 

On remarquera que les deux produits mixtes s’annulent. La nouvelle 
expression de l’écart e i sera donnée par : 

e i =<*.-( n xi u + n^ v + njju + La + M fi + N y) | 

En résumé, à un déplacement exprimé par le torseur de petits déplace- §- 
ments {2)} correspond un ensemble de mesures {e t ] déduit à partir de ô 
l’ensemble initial de mesures i& 


L i= y,- n zl - z r n yi 
M i ^z i .n xt -x r n zi 
N i - x i .n.-y j .n xt 
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Notre objectif sera de déterminer le déplacement optimal {S} permet- 
tant de minimiser l’ensemble { e ,■}. 

Note : en plus des paramètres géométriques de position, il est possible 
d'ajouter des paramètres de forme Sr, intrinsèques à la surface (comme 
une variation de rayon pour un cylindre). L'expression précédente 
devient 1 ; 


Dans la suite, nous ne chercherons pas à optimiser en intégrant ces cri- 
tères de forme (la démarche reste identique, il suffit de prendre en 
compte ces nouvelles variables complémentaires). 

Critères d'optimisation 

Minimiser l’ensemble {e t ) n’a de sens que si l’on se donne un critère de 
minimisation. 

Exemples de critères 

• Critère des moindres carrés : rendre minimum : Z 2 = 

• Critère du défaut de forme minimum (critère de Chebyshev) : 

rendre minimum : Z =Max|je.|J. 


• Critère de la somme minimale des valeurs absolues : rendre 
minimum : Z 1 = ^|e ; |. 

Application du critère des moindres carrés 
à un élément du type plan 

Rappelons que chaque valeur est donnée par l’expression : 


(§ - Sr ) - 1 n xi u + n yj v + n zj w + La + Mfl + Nj j 




soit : 



© 


Bourdet P., Schneider F., Spécification géométrique des produits. Cotation et tolérance- 
mentlSO, Paris, Dunod, 2007. 
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Z 2 sera minimum lorsque l’on aura simultanément : 


dZ 2 

da 


= 0 


dZ 2 


= 0 


dZ 2 

dw 


= 0 


Exemple didactique 

L'erreur suivant l'axe z a été fortement amplifiée (fig. 5.16). 


Figure 5.16 : Exemple de choix de critères d'optimisation 


Déformée de la surface avant optimisation Déformée de la surface avec norme Z2 






Calcul du défaut de forme (exemple de la planéité) 

Maintenant que l’on a une solution optimisée (a 0 > A)> ®o)> on P eut ca ^" $ 
culer chaque e l0 correspondant à ce petit déplacement. i, 

Le défaut de forme (planéité) sera donné par : §■ 

A Planéité Qptimisé = Sup(e i0 )-Inf(e i0 ) 
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Insistons sur le fait qu’il est étroitement lié à la méthode d’optimisation 
retenue (on peut, en particulier, rechercher un défaut de planéité 
minimum en utilisant un solveur). 


Calcul du défaut de position 

Soit la pièce réelle de la figure 5.17. 

Figure 5.17 : « Position relative » des surfaces réelles St et S 2 



À chaque surface réelle S;, nous sommes capables, suivant son type, 
d’associer un élément géométrique type (par exemple deux plans dans 
l’exemple). 

Pour la mesure de défaut de forme, nous pouvions nous contenter de 
raisonner avec les centres de sphères et non sur les points réels de contact. 
Par contre pour des mesures dimensionnelles de positions relatives, cette 
ambiguïté devra être levée. 

Par exemple, pour le plan associé, il faudra le décaler suivant la direction 
des Zp de manière à ce qu’il soit tangent à la pièce (matière). Compte 
tenu des remarques ci-dessus, on associe respectivement le plan Pj à la 
surface S ^ et le plan P 2 à la surface S 2 . 
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Souvent ces plans sont stockés dans la base de données uniquement par 
la position d’un point et d’une normale. Si la normale est parfaitement 
définie, le choix du point le plus représentatif de la surface est moins évi- 
dent (effets de bords dus aux surfaces latérales). Soient Oj et O2 ces 
deux points (fig. 5 . 18 ). 

Figure 5.18 : Définition de la position relative de deux plans 




Les logiciels de machines à mesurer tridimensionnels ne mesurent pas 
tous la position relative entre les deux éléments plans de la même façon. 
Mathématiquement, cela n’a pas de sens puisque les deux plans, non 
parallèles, se coupent (formation d’un dièdre). Il faut être extrêmement 
prudent lors de l’établissement de la gamme de contrôle (par exemple, 
en utilisant la mesure de l’erreur de position, mais alors l’ordre des sur- 
faces n’est plus indifférent). 

Mesure par calibres virtuels 

Il est possible d’imaginer la pièce réelle emprisonnée par deux plans 
parallèles, libres en rotation. Ces deux plans peuvent se rapprocher pro- 
gressivement jusqu’à une position limite autorisée par la forme de la 
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pièce. Cette technique * 1 permet d’englober les défauts de forme et de 
position. Elle présente un grand avenir pour vérifier la conformité des 
spécifications type GPS (fig. 5.19). 

Figure 5.19 : Principe du calibre virtuel 



Analyse modale d'une surface ou d'une courbe 

L’analyse modale 2 est une méthode générique d’analyse de forme, elle 
s’applique à tout type de géométrie pour obtenir automatiquement une 
base de défaut. À partir d’un nuage de points représentant la surface 
réelle dans sa position, la méthode « modale » propose, à la manière des 
séries de Fourier, de décomposer ce nuage en une somme incluant la 
valeur moyenne et les défauts associés, triés selon leur degré de com- 
plexité (position, forme, ondulation, etc.). Méthode à fort potentiel, sa 
présentation sort du cadre de cet ouvrage. Par contre, son application est 
facile à interpréter par le métrologue. Par exemple, cette approche est 
intuitive lors de mesure de circularité. 


1. Pairel É., Logiciel GeoVerif, laboratoire SYMME, www.polytech.univ-savoie.fr. 

1. Adragna P. -A., Samper S., Pillet M., Favreliere H., « Analysis of Shape Déviations of 

Measured Geometries with a Modal Basis », Journal of Machine Engineering, vol. 6, 
2006, p. 134-143. 
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Mesure de la circularité par « analyse de Fourier 1 » 

Lorsque l’on mesure une circularité, celle-ci est fortement perturbée par 
la non-coïncidence entre « l’axe » de la pièce et l’axe de rotation du pla- 
teau. Une table à chariots croisés permet de s’affranchir progressivement 
de ce problème (fig. 5.20). 

Figure 5.20 : Position du profil par rapport 
au plateau de la machine 



Plusieurs solutions permettent un dégauchissage assisté. L’analyse par 
série de Fourier en est une. 

En effet, dans la formule suivante : 

f(0) — a 0 + a x • cos(0) + b x • sin (6) + ... + a n »cos (nd) + b n *sin (nd) + ... 


1. Joseph Fourier (1768-1830), mathématicien français connu pour la décomposition des 
fonctions périodiques en séries trigonométriques et de ses études sur la propagation de 
la chaleur. 


112 


Maîtriser l'influence du posage 


le terme en æ 0 peut se corriger en déplaçant plus ou moins le compara- 
teur, et les termes a,\ et b\ correspondent respectivement aux petits 
déplacements u et v à exécuter avec la table à chariots croisés. 

Si l’on soustrait ces termes, il reste les défauts intrinsèques de la pièce 

(fig. 5.21). 


Figure 5.21 : Exemple de dégauchissage par analyse de Fourier 
Rayon = f(0) 



Cette démarche peut se généraliser à tout type de surface, au prix d’une 
complexité accrue. 


© 
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Chapitre 6 


Éléments de statistiques appliqués 
à la mesure 



Estimation d'une valeur vraie 


Le mesurande est quantifié par une valeur vraie x 0 qui ne sera jamais 
exactement connue, mais simplement estimée soit par une valeur 
unique, soit par une moyenne de valeurs. En revanche, la vraie valeur Xq 
peut être présentée comme appartenant probablement à un intervalle 
(nommé intervalle de confiance) fonction d’un seuil de probabilité 
d’appartenance fixé. 


Estimation par une moyenne 

Si l’on était capable de faire un nombre de mesures tendant vers l’infini, 
nous pourrions calculer la moyenne fl de cet ensemble. La valeur vraie 
peut être assimilée à fl, sous réserve que l’appareil ne présente pas de 
biais. Maîtriser l’erreur de justesse, revient à estimer ce biais pour pou- 
voir faire éventuellement des corrections : 

x 0 = fl + b 

Pour des raisons économiques, techniques et de temps, nous limiterons 
le nombre de mesures. Cela revient à dire que fl sera simplement 
estimée par les informations que l’on pourra retirer du ou des échan- 
tillon^) en notre possession. À partir des mesures x z - de l’échantillon, 
nous calculerons une moyenne arithmétique x, qui est une valeur parti- 
culière de X 1 . X est un estimateur de la moyenne fl de la population : 



© 


Rappel : X est une variable aléatoire. Pour chaque échantillon, on 
obtiendra une valeur particulière : 



1. NF ISO 3534-1, Statistique - Vocabulaire et symboles - Partie 1 : Termes statistiques 
généraux et termes utilisés en calcul des probabilités. 
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Dispersion de X 

Les n mesures de l’échantillon sont censées appartenir à une population 
de taille TV supposée infinie (TV» ri). La dispersion de la population de 
mesures est caractérisée par sa variance CT 2 . Cette variance sera estimée à 
partir des données de l’échantillon avec l’emploi de l'estimateur S 2 : 

s 2 = -^-X(x.-x) 2 

n 1 i = 1 

S 2 est une variable aléatoire qui dépend de l’échantillon retenu. Pour un 
échantillon particulier, nous obtenons une variance d’échantillon observée 
notée s 2 : 

v = n - 1 représente le nombre de degrés de liberté (n variables indépen- 
dantes, liées par l’équation donnant la moyenne, qui est nécessaire pour 
le calcul de s 2 ). 

Nous sentons intuitivement que la performance de l’estimateur X va 
dépendre de la taille de l’échantillon n mais également de la dispersion 
de la population mère. Cette dernière est soit connue (du fait d’un 
grand nombre de mesures préalables), soit quantifiée à l’aide de s 2 à 
partir des valeurs de l’échantillon. 

Cas où <7 est connu 

En toute rigueur, ce cas n’existe jamais. En revanche, si l’on connaît la 
reproductibilité du moyen de mesure depuis un certain temps, on peut 
admettre que ce cas est applicable. 

Partant de la loi des mesures supposée normale : 

Loi(X i ) — X(ju,(T 2 ), 
on en déduit la loi pour X : 

Loi(X) = X(jU,—) 

qui s’écrit lorsque l’on passe en variable centrée réduite : 

Loi(^éHM( 0 , 1 ) 

oyv» 
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La formule donnant la variance de X, soit 

V(X) = —, 


n’est valable que si n est petit devant N (le tirage sans remise n’affectant 
pas alors les caractéristiques successives de la population mère). 


Cela donne l’intervalle de confiance suivant ( u : variable de la loi nor- 
male centrée réduite ou loi de Gauss) : 


b \x-u a/ »-^= 
L ^ -fn 


—j=< il<X + u /•—=={ = 

x ~ a A 


OC représente le risque que l’hypothèse d’appartenance de la moyenne à 
cet intervalle ne soit pas vérifiée. Par exemple, si l’on retient un risque de 
5 %, nous obtenons : 


Prob - 1, 96 • -j= < < X + 1, 96 - -j= 




Cas où cest inconnu 

Faute de mieux, nous allons travailler avec l’estimateur S. S donne une 
estimation assez médiocre de <7, particulièrement si n est petit. L’inter- 
valle de confiance pour <7 est donné à l’aide de la loi duX 2 par : 



© 


Exemple 

Si S donne une estimation à partir de cinq valeurs, nous obtenons : 



Prob [0,60 • S < a < 2,87 • S] = 0,95 


La méthode pour obtenir l’intervalle de confiance sur la moyenne est 
dans le même esprit que dans le cas où (7 est connu, mais cela revient à 
utiliser une sorte de coefficient de sécurité en remplaçant la variable u 
par la variable t (loi de Student). 
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Prob \X-t LL<X + t „/*—r= =l-a 

L n ~ lü -Â 4n * 

Avec l’exemple précédent, cela donne : 

Prob 1-2,78*^= <u<X + 2, 78*^= =0,95 

L S S\ 

Comparaison de deux estimations 
de moyennes 

Ce cas est assez fréquent lorsque l’on veut comparer les résultats : 

► de deux machines ; 

► de deux opérateurs ; 

► de deux fournisseurs ; 

» etc. 

Paramétrage du problème 

Chaque contexte de mesure d’indice 1 ou 2 est caractérisé par un échan- 
tillon propre (tableau 6.1). 

Tableau 6.1 : Paramètres des deux sources de mesures 


Taille 

Moyenne 

Écart-type 

Moyenne : estimation 

- 

Ecart-type : 

échantillon 

population 

population 

échantillon 

estimation échantillon 

n i 

m 

o\ 


si 

n 2 

P-2 

02 

x 2 

s 2 


Les variables aléatoires d’échantillon et X 2 seront supposées être dis- 
tribuées suivant une loi normale. Les échantillons de mesures ne sont 
pas appariés (c’est-à-dire qu’il n’y a pas de lien par exemple ente la f™ 
valeur de l’échantillon 1 et lay'™ f valeur de l’échantillon 2). 

Notre objectif est de savoir si la différence observée des moyennes notée 
D est significative ou non. 
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Nous retenons a priori l’hypothèse nulle H 0 (/ij = /l 2 ). D est une 
variable aléatoire : 

D-X l -X 2 avec : E(D) — 0 (Hypothèse H 0 ) 



Test de comparaison 

D étant une variable aléatoire, nous devons prendre une décision à partir 
de sa valeur observée d. Nous considérerons que nous sommes dans le 
cas ou CT j et(J 2 ne sont pas connus (donc a priori différents). Pour pou- 
voir faire une comparaison avec la variable t de Student, nous utiliserons 
la variable suivante : 


Z sera estimée par 2 — — . 

S D 

La distribution de Z reste très proche de la distribution de Student, on 
considérera que l’hypothèse H 0 sera acceptée si 



s D est calculé à partir des estimations 5] et r 2 ■ 



Le nombre de degrés de liberté V est le plus souvent approché par 
V = n x + n 2 - 2. 

Pour plus de rigueur, on peut utiliser le test de Welch 1 : 



1. Pagès J., Statistiques générales pour utilisateurs, PUR, 2005. 
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Comparaisons d'échantillons appariés 

Ce cas se produit par exemple lorsque qu’un échantillon de pièces est 
mesuré par deux personnes différentes ou sur deux moyens de mesures 
différents. De plus, chaque pièce est identifiée, si bien qu’il est possible 
de comparer les deux mesures obtenues pour chaque pièce. C’est un cas 
particulier de l’analyse de variance à deux facteurs sans répétition (le fac- 
teur pièce peut présenter plusieurs modalités, par contre l’autre facteur 
n’a que deux modalités). 

Échantillon des différences 

Pour chaque pièce i, on calcule la différence 

d i — x xi — x 2 j- 

De cet échantillon de différences, on calcule les estimations suivantes : 


Recherche de l'influence d'un facteur 
(analyse de la variance à un facteur) 

Influence du facteur identifié 

L’objectif est de savoir si les différences de moyennes lorsque le facteur 
varie (modalités), sont imputables au hasard (hypothèse nulle : HO) ou 


d = ~%d t 



2 



L’hypothèse nulle H 0 sera acceptée (risque a) si 


d 
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si au moins une des moyennes diffère des autres (hypothèse : Hl). Par 
exemple, le facteur peut être un changement : 

► d’opérateurs ; 

► de pièces ; 

► de machines ; 

* de posages ; 

► d’atelier (environnement) ; 

* etc. 


Notations et présentation 

Nous avons : 

* n mesures (échantillon total) ; 

► L niveaux du facteur d’influence étudié (modalités) ; 

► r répétitions par modalité (r sera constant pour chaque modalité, 
donc n = L - r). 


Tableau 6.2 : Présentation des données 
Répétitions 

1 ... j ... r Moyennes 



yij est une valeur particulière de la variable Y-. Cette variable suit la loi : 
Loi{Yfl = J{(jUj,<J 2 ) 

Nous supposerons que le facteur étudié agit seulement sur la moyenne. La 
non-influence sur la variance G 2 doit être vérifiée par un test ( Cochran, 
Hartley, Bartlett. . .). 
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Chaque variable Y- peut être représentée par une valeur moyenne, 
fonction de l’influence du facteur et un écart aléatoire (dispersion de 
mesure) : 

Y ÿ ~ fli + e ÿ 
Loi(e..) = >f(0, CT 2 ); Vz, y 

Cm{e-,e ü ) — 0;V (*',_/) -£ (£,/) les mesures sont indépendantes. 

Pour mettre en évidence l’influence (a ? ) explicite de la modalité i du 
facteur, on retiendra l’écriture suivante {fl représentant la moyenne 
globale) : 

Y t] =ii + a t +£ t] 

avec : = 0. 


Estimation des paramètres 

fl est estimée par J. Avec la notation précédente, il vient : 

a i est estimé par à i = J ■— J 
Estimation de la variance des résidus : 

II (>,-*? ny 2 

â 2 = s 2 =^ — 

n — L n — L 

ou encore, comme ici n — r-L: 

Hfef I 2 ? 

2 _ > j J_ J 

Z(r-l) L 1 

Q- 

ô 
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Étude de la variabilité 
(sommes des carrés des écarts) 

À chaque résultat de mesure y^, nous pouvons associer deux écarts : 

► un écart à la mesure moyenne du facteur (influence pour la modalité î) : 

un écart à la moyenne générale : 

yÿ-ÿ 

qui peut se mettre sous la forme 

(y 9 +(-ÿi + ÿi)~ÿ)- 

L’équation d’analyse de variance donne : 

i j i t j 

ou encore : 

lb,-ÿh^*y, 

Chaque terme représente une somme de carrés d’écarts, notée SCE. 
Soit : 

SCE j ( variabilité totale) = SCEp ( variabilité due au facteur) 

+ SCEpfaduelU (variabilité résiduelle). 

À chaque somme correspond un nombre de degrés de liberté ( ddt) 
(tableau 6.3). 


Tableau 6.3 : Degrés de liberté des sommes d'écarts 


SCE 

Écart 

Somme 

Relations 

ddl ( v) 

sce t 

h-y) 

lï{y«-ÿf 

XX(y,y-ÿ)=o 

vr=n- 1 

sce f 

(y,-ÿ)=« ; 

'■ïfo-yf 

X«;=° 

v F = L - 1 

sce r 



V/ — > X e ;/ = 0 

v R = n - L 
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Coefficient de détermination R 2 

La variabilité se décompose en deux parties : 

» une partie non maîtrisée, variabilité résiduelle, correspondant à 
SCE r ; 

» l’autre partie subissant l’influence du facteur étudié : SCEp. 

Dans le cas particulier, de l’analyse de la variance à un facteur, cela 
donne : 

r2 _ SCE t -SCE r _ SCE F 
SCE T " SCE T 

Si ce rapport est égal à un, cela revient à dire que 1 00 % de la variabilité 
sont expliqués par le facteur d’influence. 

Pour savoir, si la part de variabilité attribuée à la grandeur d’influence 
est significative, nous allons utiliser le test suivant. 


Test sur l'effet supposé du facteur 

Carrés moyens 

En divisant chaque somme de carrés des écarts par le nombre de degrés 
de liberté associés, on obtient un carré moyen d’écarts (CME) représen- 
tant physiquement une variance. 

Pour estimer la variance inter-lots (variance résiduelle), nous avons 
écrit : 


XXi.v, XXi-,i 


Les lots étant de taille identique, nous pouvons utiliser également : 


CME R correspond à l’estimation de la variance résiduelle. 


I 

e 


o 
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De même, posons : 



r-m-ÿf m-ÿf 


CME F — 


L - 1 


L - 1 


Si l’on retient l’hypothèse H 0 , cela revient à dire que le changement de 
modalités n’a pas d’influence. Avec cette hypothèse, les différentes 
modalités peuvent être considérées comme des échantillons appartenant 
à la même population. La variance sur les moyennes d’échantillons est 
estimée par : 


Donc, sous l’hypothèse H 0 , l’espérance de CMEp devrait être : 

O 2 = CME r . 

Le test consiste à comparer CMEp à CMEp. Pour cela, on calcule le 
rapport : 



CME F 


L - 1 


L’espérance de r^est égale à — . 



Si H 0 est vraie, la variable FQp ervc ; suit la loi : 


w(w=0' a) - 


Loi ^ F Observé'. 


L’hypothèse H 0 sera retenue si : 



!(l-a) 


(test de Snedecor). 
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Exemple : mesure d'une cale sur machine à mesurer 

Suite à l'acquisition de mesures d'une cale-étalon sur machine à 
mesurer tridimensionnelle, on a des doutes sur l'incidence de la position 
de la cale par rapport à la machine. En effet, ce n'est pas les mêmes 
masses de chariots mises en mouvement suivants les axes utilisés 
(tableau 6.4). 


Tableau 6.4 : Mesure d'une cale en faisant varier sa position 


Mesure d'une cale de 75 mm suivant différentes orientations sur MMT 


n° 1 

n° 2 

n ° 3 1 n ° 4 

P * « 

Plan YZ, mesure axe Z 

75,002 

75,002 

75,001 

75,001 

75,000 

75,0012 

0,0008 

Plan YZ, mesure axe Y 

74,996 

74,996 

74,996 

74,995 

74,996 

74,9958 

0,0005 

Plan XY, mesure axe X 

74,994 

74,996 

74,994 

74,996 

74,996 

74,9952 

0,001 1 

Plan XY, mesure à 45° 

74,999 

74,998 

74,998 

74,997 

74,999 

74,9982 

0,0008 


Résultats de l'analyse de variance à un facteur (la position de la cale) 
(tableau 6.5). 


Tableau 6.5 : Résultats de l'analyse de variance 


Source 

ddl 

SCE 

CME 

MMT (résiduel) 

16 

0,00001120 

0,00000070 

Position 

3 

0,00011160 

0,00003720 

Total 

19 

0,00012280 

0,00000646 


R 2 = 

90,88 % 

F observé = 

53,143 

F calculé = 

3,24 


Remarque : l'utilitaire d'analyse Excel™ permet d'obtenir directement 
cette étude (tableau 6.6). 


Tableau 6.6 : Analyse directe avec Excel™ 


ANALYSE DE VARIANCE 

Source des 
variations 

Somme 

carré 

Degré 

liberté 

Moyenne 

carrés 

F 

Probabilité 

Valeurcritique 
pour F 

Entre groupes 

0,0001116 

3 

3,72E-05 

53,143 

l,53E-08 

3,23886695 

À l'intérieur 
des groupes 

1,12 E-05 

16 

7E-07 
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Conclusion : l'hypothèse Hq est refusée, la position de la cale semble 
avoir de l'influence sur la mesure. 

Test non significatif (H 0 est acceptée) 

Lorsque H 0 est acceptée, cela revient à dire qu’il n’y a pas de détection 
d’influence du facteur. On peut considérer que l’on a un seul gros 
échantillon de taille n = L ■ r : 


2 _ SCE r 

Sr ~~^Ï 

soit comme intervalle de confiance pour la variance : 


V R 2 < ^2 < _R_ 2 

2 S R-°R~ 2 S R 

X a Xa 


Z 

La moyenne générale /i est estimée par y avec une variance —p 
Soit comme intervalle de confiance sur la moyenne : 

où t est la variable de Student à n — 1 degrés de liberté. 


Test significatif (H 0 est refusée) 

Nous avions modélisé la mesure par : Y- — jl + ct l + £-■ 

Si Hq est vraie, c’est que l’on pense avoir : a — 0; Vf. 

Dans le cas contraire, c’est que manifestement, il existe au moins un 
«.très différent de 0 (au moins une modalité du facteur F a une grosse 
influence). 
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Il va falloir distinguer une variance résiduelle al et une variance due au 
facteur a p . 

La variance résiduelle est estimée par : 


avec :n — L degrés de liberté. 

Soit comme intervalle de confiance pour la variance résiduelle : 


Prob 




1 -a 


1 -a 


Pour l’influence du facteur, nous ferons la distinction entre deux cas : 

* soit la modalité du facteur est aléatoire (par exemple, position de la 
cale au hasard, choix des opérateurs. . .) ; 

* soit il s’agit de modalités fixes (par exemple, comparaison inter-labo- 
ratoires). 


Modalités aléatoires 

CME p a. comme espérance : E(CME P ) -a R + r-o p . 
La variance o 2 P peut être estimée par : 

2 _ CME p -CME r _ CME F — s R 

s F ~ - 

r r 

soit comme intervalle de confiance pour la variance : 



r \ 


1 

CME P 2 


r 

M SR 

l^-f) J 


I 

e 


o 
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Application à la mesure de la cale : nous considérerons que la position 
de la cale sur la machine est aléatoire et dépend du choix arbitraire de 
l’opérateur. 

S Position ~ 0 . 0027 t ? 27?2 

Prob[0.0013 < G Position < 0.0055] = 0.95 


Modalités fixes 

Dans ce cas, le choix de la modalité n’est pas dû au hasard, il relève 
d’une volonté de comparaison (par exemple, comparaison inter-labora- 
toires). 

Les effets de la modalité n° i sont notées CL i et estimées par ôt i =J.—J 
On peut rechercher si une modalité a un effet très fort par rapport aux 
autres ou si l’on peut classer les modalités en groupes d’influence com- 
parable (dans ce dernier cas on peut utiliser la méthode de Tukey- 
Kramer). 

Nous allons envisager deux hypothèses pour la modalité OC l : 

► hypothèse H 0 : la modalité n’a pas d’influence (a z = 0) ; 

) hypothèse Hj : la modalité est nettement différente de la valeur 
moyenne (a z - <> 0). 

Rappelons que oc { est estimé par d t — ~ÿ i — J, qui peut se mettre sous la 
forme : 



Elle est plus faible que la variance sur la moyenne des modalités (parti- 
culièrement quand L est petit). 

Comme nous travaillons avec un estimateur, nous utiliserons la loi de 
Student : 
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L’acceptation de H 0 se fera à l’aide du test suivant (avec le risque (X) : 


Ct; 

s» 


<t 



Remarque : si l'on a une valeur vraie du mesurande (notée yg), l'effet de 
la modalité (erreur systématique) est estimée par â. = y. - y Q . 


Recherche de l'influence de deux facteurs 
(analyse de la variance à deux facteurs) 

Nous devons traiter des données de mesures, obtenues pour différentes 
modalités de deux facteurs. Ce cas est souvent rencontré lorsque l’on 
cherche des influences potentielles sur la reproductibilité. Par exemple, 
dans les études « GRR » (chapitre 7) on cherche à connaître l’influence 
de la pièce et de l’opérateur sur la dispersion des mesures. 
Expérimentalement (plan factoriel), on se fixera un nombre de moda- 
lités pour chaque facteur (par exemple trois opérateurs mesurant dix 
pièces). À chaque combinaison de modalités (un opérateur/une pièce) 
correspond une ou plusieurs mesures (r répétitions) (tableau 6.7). 
Chaque mesure y^ r est une occurrence de la variable Y- repérée par trois 
indices : 

► i : indice de ligne (par exemple, opérateur A, B ou C) ; 

* j : indice de colonne (par exemple, numéro de la pièce mesurée) ; 

► r : indice de répétition. 
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Tableau 6.7 : Présentation des données avec deux facteurs 
Modalités facteur 2 














(y/zv-ov-y^Hy;/. 














Convention d'écriture des différentes moyennes 

» Moyenne des R répétitions pour la combinaison modalité i du fac- 
teur 1 et modalité y du facteur 2 : 

y g» ~ ^ 2 * 

► Moyenne (QR) de toutes les mesures correspondant à la modalité i 
du facteur 1 : 

* Moyenne (PR) de toutes les mesures correspondant à la modalité^ du 
facteur 2 : 

ÿ.j. = -^ly yr 

» Moyenne générale (PQR) de toutes les mesures : 

j7 *** ” PQR^/ ,jr 
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Modélisation 

Pour mettre en évidence les influences explicites (a,) de la modalité i du 
facteur 1 et (jGy) de la modalité^ du facteur 2, on retiendra l’écriture sui- 
vante (/ 1 représentant la moyenne globale) : 

Y ijr = n + a i +P j + aP ÿ + e ii r 
Loi(e. v ) = Jvf(0,C7 2 );Vi,j,r 

Notion d’interaction : sia/? est different de zéro, on parlera d’interac- 
tion entre les modalités i et j. Par exemple, pour la mesure de la pièce i 
par l’opérateur y, cela peut représenter la recherche systématique du dia- 
mètre maximum pour une pièce présentant un défaut d’ovalisation. 

Analyse de la variabilité 

Dans ce chapitre, nous nous contenterons d’étudier sommairement la 
décomposition de la variabilité. Pour une étude plus approfondie, nous 
renvoyons le lecteur aux ouvrages spécialisés. 

La variabilité totale ( SCEj ) se décompose dans ce cas en quatre parties 
(tableau 6.8). 


Tableau 6.8 : Sommes des carrés d'écarts et carrés moyens d'écarts 


Source 

SCE 

Ddl(v) 

CME 

Totale 

III ( y ,-,- r - ÿ ...) 2 

vr=PQR- 1 


Facteur 1 

Ig%„. -y...) 2 

Vfi =P- 1 

— & 

Facteur 2 

I pr[y.,. - y ...) 2 

Vr2 = Q - 1 


Interaction 

«II(y, r -y,., -y.,. +y...) 2 

V F,F 2 = (P- b * (Q~ b 

sœ f 

CME f i f2 -( p _i)( Q _i) 

Résidus 

III(y i/r -y ; ,) 2 

vr = (R - 1 ) * (PQ) 

CMEr = PQ(R - 1) 


o 
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On peut vérifier que : 


et : 


SCE t — SCE Fl + S CE p 2 + S CE p- lF2 + SCE R 


Test sur la validité de l'influence des facteurs 


À partir des carrés moyens d’écarts calculés précédemment, on va com- 
parer une valeur F 0 y servé à un F seui i (test de Snedecor) : 

F observé 

Nous allons vérifier : 


CME H 2 Q_i 
cme r 

CMEp XF2 (jP _ 1 )( q_ i; 

CME r ~ mR ~ x) 


Si Foy servé est plus grand que la valeur seuil, l’hypothèse de non- 
influence sera rejetée. 


Exemple : analyse type G RR 

Nous allons analyser une série de dix pièces, mesurées par trois opéra- 
teurs, chaque mesure étant répétée trois fois (tableau 6.9). 

Tableau 6.9 : Exemple d'étude GRR 



Pièce 1 

Pièce 2 

Pièce 3 

Pièce 4 

Pièce 5 

Pièce 6 

Pièce 7 

Pièce 8 

Pièce 9 

Pièce 

10 

Opérât. 

87,830 

89,770 

90,140 

89,750 

89,440 

87,550 

88,470 

89,570 

88,160 

87,350 

87,420 

89,680 

90,160 

89,450 

89,380 

87,520 

88,030 

89,450 

88,700 

86,970 


87,170 

89,540 

90,180 

89,110 

89,590 

87,490 

88,150 

89,590 

88,600 

87,430 


88,120 

89,870 

90,670 

89,820 

89,060 

87,630 

88,370 

89,300 

88,850 

87,730 

88,080 

89,970 

90,430 

89,580 

89,040 

87,560 

88,410 

89,430 

88,890 

87,460 


87,955 

90,070 

90,440 

89,530 

89,190 

87,640 

88,540 

89,430 

89,150 

87,610 


87,920 

89,970 

91,040 

89,510 

89,360 

88,250 

88,700 

89,750 

89,190 

87,570 

88,050 

90,030 

90,810 

90,140 

89,610 

87,980 

88,940 

89,760 

89,010 

87,790 


87,580 

90,160 

90,950 

89,840 

89,580 

88,120 

88,620 

89,770 

88,860 

87,700 
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En reprenant le tableau 6.8, nous obtenons les résultats suivants 
(tableau 6.10). 


Tableau 6.10 : Exemple d'étude GRR 


Source 

SCE 

ddl(v) 

CME F observé F seyii 

Totale 

91,6855 

89 

1,030 



Facteur « Opérateurs » 

2,7838 

2 

1,392 

48,983 

3,150 

Facteur « Pièces » 

85,6567 

9 

9,517 

334,933 

2,040 

Interaction 

1,5400 

18 

0,086 

3,011 

1,778 

Résidus 

1,7050 

60 

0,028 




Les deux facteurs sont bien influents et l'on peut noter une interaction 
non nulle entre ces deux facteurs. Cela se manifeste en figure 6.1, par 
des courbes qui ne sont pas simplement décalées. 


Figure 6.1 : Étude de l'interaction opérateurs/pièces 



Cette exemple sera repris au chapitre 7. Notons que nous pouvions 
avoir directement ces résultats en utilisant l’utilitaire d’analyse d’Excel™ 
en choisissant « Analyse de variance : deux facteurs avec répétition 
d’expérience » (tableau 6.1 1). 

I 
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Tableau 6.1 1 : Exemple direct avec Excel™ 


ANALYSE DE VARIANCE 

Source 

des variations 

Somme carrés 

Degré 

liberté 

Moyenne 

carrés 

F 

Probabilité 

Valeur critique 
pour F 

Echantillon 

2,783787222 

2 

1,391893611 

48,983 

2,441 32E-1 3 

3,5150411487 

Colonnes 

85,65668472 

9 

9,517409414 

334,933 

7,87534E-48 

2,040096092 

Interaction 

1,540029444 

18 

0,085557191 

3,0109 

0,000711265 

1,778445835 

À l'intérieur 
du groupe 

1,70495 

60 

0,028415833 




Total 

91,68545139 

89 






Régression linéaire 

Cet outil statistique sera intéressant pour comparer deux séries de 
mesures d un même échantillon de pièces , obtenues dans des contextes de 
mesure différents (changement de machine, de lieu, etc.). Par exemple, 
il permet de calculer les corrections de justesse en comparant les résul- 
tats d’une machine à mesurer étalon, en salle de métrologie, avec ceux 
donnés par une machine à mesurer en production (nous reviendrons 
plus en détail sur ce point aux chapitres 7 et 8). 

Pour une même pièce n° i, nous reporterons la valeur x i de la machine 
d’atelier en abscisse et la valeur y z de la machine à mesurer étalon en 
ordonnée, comme le proposent C. Perruchet et M. Priel 1 . Cette disposi- 
tion n’est pas naturelle. Pourtant, elle est justifiée car la machine d’ate- 
lier devrait fournir un résultat correspondant à la machine étalon. x z - 
n’est pas une valeur aléatoire, c’est la valeur lue. Par contre, pour une 
même valeur Xj, on peut s’attendre à différentes valeurs y; correspon- 
dantes. Cette dernière est une valeur aléatoire. y t est une valeur particu- 
lière de la variable Yi (fig. 6.2). 


1 . Perruchet C., Priel M., Estimer l’incertitude, AFNOR, avril 2000. 
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Figure 6.2 : Régression linéaire 



Nous retiendrons comme modèle : y i — (X ■ x i + /3 + £. 

£ représente les résidus. C’est une variable aléatoire supposée distribuée 
par une loi normale ^V~(0, Cf). 

Malheureusement, ce modèle ne pourra qu’être estimé. En effet, nous 
n’avons qu’un échantillon de points (couples Xj, y-). À l’aide de cet 
échantillon, nous proposons un modèle estimé: y. - a-x +b + e-, en 
cherchant à minimiser ^ ej (méthode des moindres carrés). 

En résumé, a et b ne sont que des estimations de a et /3. La « qualité » de 
cette estimation sera fortement influencée par le nombre de points et 
l’étendue de mesure (valeur maximum — valeur minimun ) . 


Estimation des paramètres 

À partir des moyennes x et ÿ, on en déduit : 


Cov(x,y) 
f Var(x) 
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Intuitivement, des machines identiques dans leur comportement 
devraient donner a = 1 et b = 0. 

* Si b est différent de zéro, cela correspond à une erreur systématique 
correspondant à un décalage d’origine. 

► Si a est différent de 1, cela correspond à un changement d’échelle 
(métriques différentes suivant les moyens de contrôle). 


Incertitudes sur les paramètres 

Il est évident que l’étendue de mesure, ou plage des « x », va améliorer la 
qualité des coefficients a et b mais également la distribution des points. 
Cela se traduira par l’emploi de l’indicateur « variance sur les x ». 
D’autre part a et b ne sont pas indépendants (covariance non nulle). En 
utilisant l’écart-type <7 des résidus, cela donne : 



O n’est en général pas connu. On devra travailler avec son estimation s. 
Ysbi-a-Xi-b) 2 

s 2 = — (il y a deux paramètres estimés). 

n - 2 

Intervalle de confiance de la correction 

L’utilité de l’étude précédente est de prévoir une correction sur la 
mesure brute x* sans avoir recours systématiquement à une comparaison 
à la machine étalon. En résumé, quelle est la confiance que l’on peut 
avoir dans le résultat y* calculé en fonction de la mesure brute x*. 

Dans le résultat donné par y * = a ■ x* + b. 
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rappelons que nous ne sommes pas sûrs des coefficients a et b, soit : 
1 


+ r - 


i+— 


+ 2-(*.)-W-T 


Le dernier terme correspond à l’influence de la covariance entre a et b 
(pour plus de détails voir le chapitre 8 sur le calcul d’incertitudes). 


:^-\xl+sl + x 2 

n ■ s 2 x L J 


Nous venons de quantifier l’incertitude propre à la méconnaissance des 
facteurs Cf et /?. Cela permet de calculer un intervalle de confiance sur la 
moyenne des y* pouvant correspondre à la valeur x*. Pour une valeur 
individuelle, il faut rajouter la variance propre et utiliser : 




ou encore : 


1 (x t -x 


L’intervalle de confiance (fig. 6.3) (pour une probabilité de refus de Cf) 
est donné par : 


>-*) 2 


I 

e 
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Figure 6.3 : Exemple d'intervalle de confiance pour la correction 



Invervalle de confiance pour E(y*) 
Invervalle de prévision pour y* 
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Chapitre 7 


Valider un processus 
de contrôle 



Pour avoir confiance dans une mesure, il faut que l’ensemble des res- 
sources humaines et matérielles du processus de contrôle conduisent à : 

► une mesure fidèle (le processus réitère des résultats voisins) ; 

► une mesure juste (la moyenne de ces résultats s’approche d’une 
valeur vraie). 

Dans ce chapitre, nous chercherons à valider un processus de contrôle 
dédié à une tâche bien définie. Il s’agit d’appareils utilisés en production 
de série. Pour ceux d’emploi plus universel (un pied à coulisse par 
exemple), la justesse globale (sur toute l’étendue de mesure) sera estimée 
en utilisant l’outil « régression linéaire » vu aux chapitres 6 et 8. 


Concept de fidélité 

Cela sera perçu comme l’aptitude à donner des résultats de mesure voi- 
sins pour un même mesurande. Ce degré de voisinage est quantifié par 
une statistique du type écart-type ou étendue de mesure. Pour évaluer la 
robustesse de cette fidélité, il est judicieux d’envisager l’influence de 
diverses grandeurs pouvant perturber la qualité de cette statistique. Ces 
grandeurs constituent l’ensemble des 5M 1 du processus de mesure : 

► moyen de contrôle : l’appareil (avec ses étalons) ; 

& main-d’œuvre : l’opérateur (dextérité physique et compétence) ; 

► matière : la pièce (problèmes de la géométrie, de la forme, de l’état de 
surface, etc.) ; 

» méthode de mesure (protocole, traitement, etc.) ; 

► milieu (température, vibrations, hygrométrie, etc.). 


1 . Duret D., Pillet M., Qualité en production : de l’ISO 9000 à Six Sigma, Paris, Éditions 
d’Organisation, 2005. 
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Il est pratiquement impossible de prendre en compte toutes les 
influences, soit simplement par ignorance, soit parce qu’on a estimé que 
cette influence particulière engendrait une variabilité très faible. 

Lorsque les influences sont figées et que l’on ne prend en compte que 
l’influence de l’équipement de mesure sur une courte période de mesure 
(variabilité résiduelle), on définira la fidélité par la répétabilité 1 . Si au 
moins une des influences varie au cours des mesures, on parlera alors de 
reproductivité. 

Répétabilité des résultats de mesurage ( repeatability ) : étroitesse de 
l’accord entre les résultats des mesurages successifs du même mesu- 
rande, mesurages effectués dans la totalité des mêmes conditions de 
mesure : 

- même mode opératoire ; 

- même observateur ; 

- même instrument de mesure utilisé dans les mêmes conditions ; 

- même lieu ; 

- répétition durant une courte période de temps. 

Reproductibilité des résultats de mesurage ( reproducibility ) : étroi- 
tesse de l’accord entre les résultats des mesurages du même mesu- 
rande, mesurages effectués en faisant varier les conditions de mesure 
décrites ci-dessus. 

Attention, car certaines entreprises limitent la reproductibilité ( repro- 
ducibility ) uniquement au complément dû à l’opérateur {appraiser 
variation). La prise en compte du temps est alors nommée stabilité 
C stability/drift ) (fig. 7.1). 


1 . Les termes spécifiques à la métrologie (comme grandeur ou unité) sont définis dans le 
VTM ( Vocabulaire international de la métrologie), PR NF ISO 99999 octobre 2006 
(VIM), 3 e éd. 
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Figure 7.1 : Deux perceptions du concept de reproductibilité 



Méthodes Milieu 


Définition VIM 
Main-d'œuvre 



Industriellement, deux approches sont souvent retenues 1 , à savoir la 
méthode du GUM 2 et la méthode GRR du MSA 3 . 


Méthode du GUM 


Très pratique lorsque le résultat final Y dépend de mesures de bases A z - 
(entrées), elle nécessite de connaître la relation fonctionnelle : 


Y=AX X ,X 2 x n ) 


1. Le lecteur pourra également consulter les normes CNOMO : E41.36.110.N et 
GE41-013N, téléchargeables sur www.cnomo.com. 

2. Guide to the Expression of Uncertainty in Measurement (GUM), NF ENV 13005 août 
1999. L’AFNOR publie ce document sous la référence XP X 07-020. 

3. Automotive Industry Action Group (AIAG) (2002). Measurement Systems Analysis 
Reference Manual, 3 e éd., Chrysler, Ford, General Motors Supplier Quality Require- 
ments Task Force ( www.carwin.co.uk/qs ). 
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En particulier, une mesure de base brute peut être affinée en appliquant 
des corrections c z - (vues comme des entrées). À la mesure corrigée, on 
peut associer une incertitude obtenue à partir du développement de f 
autour des moyennes des c ,■ en série de Taylor du premier ordre (étude 
reprise en détail au chapitre 8). 

M cor = M bn* + 'L C i 

rt M cor)=°M bm '+'L° 2 ( C i) 

L’écriture précédente suppose que les corrections ne sont pas corrélées. 
Lorsqu’une correction n’est pas prise en compte, cela revient à faire une 
correction nulle. Ce qui ne veut absolument pas dire qu’il faut négliger 
l’incertitude associée à cette correction. 

Les mesures corrigées sont plus du domaine des laboratoires. En pro- 
duction, pour un résultat immédiat, on demandera à l’équipement de 
donner directement le résultat. En cas de corrections, celles-ci devront 
être intégrées à ce dernier. 

Méthodes GRR 

Très employées dans l’industrie automobile, elles sont basées sur l’ana- 
lyse de variance à deux facteurs avec répétition d’expériences. On retiendra 
les facteurs d’influence : 


I 

e 


o 
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► opérateurs ; 

» pièces. 

En réalité, l’influence de la méthode est incluse dans l’influence de 
l’opérateur au cas où la procédure de mesure ne serait pas parfaitement 
explicite. Les autres influences seront masquées dans la variance rési- 
duelle et s’ajouteront à la variance intrinsèque de l’appareil, c’est pour- 
quoi l’on parlera d’incertitude de l’équipement de mesure. 

Il existe principalement deux approches : 

» GRR type « average and range method ». Cette approche suppose 
qu’il n’y a pas d’interactions entre les pièces et les opérateurs. Chaque 
mesure peut être modélisée par : T y = /i + a • + / 3j + e- r , c’est-à-dire 
que la mesure de la f eme répétition est influencée par l’effet de la 
pièce (0Cj) et par l’effet de l’opérateur (fij) . Il reste une part aléatoire 
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£ÿ r non identifiée explicitement qui vient s’ajouter au résultat 
moyen / 1 ; 

► GRR type « ANOVA ». L’interaction « pièce-opérateur » peut se ren- 
contrer du fait de méthodes de mesures propres à chaque opérateur. 
Par exemple, en cas d’ovalisation d’un alésage, un des opérateurs 
recherchera systématiquement le diamètre maximum, alors que, 
pour un autre, c’est la valeur minimum qu’il faut retenir. La modéli- 
sation sera légèrement différente, elle prendra en compte cette 
interaction : Y~ r - pi + a i + pj + ap- + £- r - 

Pour les deux types d’études, nous associerons au processus de mesure 
un écart-type global appelé <JqrR qualifiant l’incertitude de reproducti- 
bilité. 

Notations des différents écarts- types recensés : 

* ®GRR '• variabilité de reproductibilité de mesure ; 

* Opy : variabilité de répétabilité (équipement de mesure) ; 

► (?AV : variabilité due à l’opérateur (au niveau de la mesure) ; 

» <7py: variabilité du moyen de production (si non connue par MSP 1 , 
elle sera estimée à l’aide de l’échantillon de dix pièces par exemple) ; 

* O j y: variabilité totale (production & mesure). C’est celle qui est 
estimée en MSP. 

GRR type « average and range method » 

Cette méthode se veut simple et pragmatique. Nous la présenterons à 
l’aide de l’exemple donné en figure 7.2 sous forme de feuille de calcul. 

Mode opératoire 

On retient un nombre de n pièces différentes sélectionnées au hasard (ce 
qui laisse entendre qu’elles ne sont pas prélevées sur un intervalle trop 
court de production). Elles doivent donner, si possible, une bonne 
image de la dispersion globale du moyen de production. Ces pièces 


1. Pillet M., Appliquer la maîtrise statistique des procédés, Paris, Éditions d’Organisation, 
4 e éd., 2005. 
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seront référencées et mesurées par k opérateurs, sans que ceux-ci voient le 
numéro d’identification. Cela permettra de faire mesurer les pièces plu- 
sieurs fois par chaque opérateur (r répétitions ) en évitant des mesures 
successives sur une même pièce (il faut éviter la perturbation due à un 
« effet mémoire » de l’opérateur). 

Pour pouvoir utiliser cette méthode sans correction particulière, il est 
demandé d’avoir un nombre d’opérateurs k et un nombre de pièces n 
vérifiant la relation suivante : 


\k»n 


> 16 1 


Formules statistiques de base 

La dispersion de mesure est le plus souvent estimée à 6. (T (ce qui corres- 
pond, pour une distribution gaussienne centrée, à une probabilité de 
99,73 %). Il est bien sûr possible de retenir un autre coefficient si on le 
désire. 

L’écart-type est calculé directement à partir des moyennes d’étendues 
par une formule du type : 



où R~ est calculée à l’aide de N sous-groupes (N = k opérateurs X n 
pièces). Chaque sous-groupe (pièce i mesurée par l’opérateury) a comme 
étendue R~ (obtenue à partir des r répétitions de mesure Xÿ^) . 


Tableau 7.1 : Tableau coefficient 62 


Nombre 
de répétitions 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Coef c/j 

1,128 

1,693 

2,059 

2,326 

2,534 

2,704 

2,847 

2,970 

3,078 


Si N <15, on retiendra un coefficient légèrement modifié appelé 

d\\ 


1. Duncan A.J., Quality Control and Industrial Statistics,4 c éd., Homewood, IL, Richard o 

D.Irvin Inc., 1974. ® 
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Tableau 7.2: Tableau coefficient d| 



Nombre N de sous-groupes 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

2 répétitions 

1,414 

1,2/9 

1,231 

1,206 

1,191 

1,181 

1,173 

1,168 

1,163 

1,160 

1,157 

1,154 

1,153 

1,151 

1,149 

3 répétitions 

1,912 

1,806 

1,769 

1,750 

1,739 

1,731 

1,726 

1,722 

1,719 

1,716 

1,714 

1,712 

1,711 

1,710 

1,708 


Calcul des moyennes et étendues 

Rappelons que pour chaque pièce i mesurée (en général) deux ou trois 
fois par l’opérateur y, on calcule l’étendue Rÿ. 

Pour chaque opérateur y, on obtient : 

* la moyenne de ses mesures : 



la moyenne des différentes étendues : 



Avec la même démarche que celle utilisée en analyse de la variance 
(chapitre 6), en comparant les carrés moyens CM Résiduel et CM Qp’ rateur , 
nous allons utiliser : 

la moyenne des moyennes d’étendues : 



l’étendue des moyennes des opérateurs : 

- X Maxi ~ X mmï 


O 
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Figure 7.2 : Exemple d'étude G RR (méthode moyenne et étendue) 
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Validation des mesures 

On peut utiliser les méthodes de la maîtrise statistique des produits pour 
vérifier que la différence de mesures obtenue par opérateur pour une 
même pièce reste dans des limites raisonnables. 

Une différence trop importante entre les mesures d’une même pièce 
pour un opérateur donné (étendue forte) peut détecter un problème 
avec l’appareil ou un problème d’utilisation par l’opérateur. Récipro- 
quement, des différences trop voisines de zéro sont douteuses (appareil 
bloqué ou en butée). 

La limite supérieure de contrôle est donnée par la formule : 

LSC r 


avec : 


Nb mesures/pièce/opérateur 

2 

3 

d a 

3,267 

2,574 


Figure 7.3 : Exemple de surveillance de la dispersion par opérateur 
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Répétabilité (Equipment Variation) 

La répétabilité de l’équipement de mesure est calculée à partir de la 
moyenne des moyennes d’étendues : 

R 

°EV ~ , 

^2 

Nous avons k ■ n échantillons de taille 2 ou 3 (répétitions). 
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Incidence de l'opérateur sur la reproductibilité 
( Appraiser Variation ) 

Nous avons calculé l’étendue sur les moyennes des opérateurs : 


Nous pouvons en déduire l’écart-type sur les moyennes d’échantillons, à 
l’aide du coefficient : 


Note : nous retenons le coefficient 62 car nous avons un seul sous- 
groupe de deux ou trois moyennes (dans l'exemple traité, un sous- 
groupe de trois moyennes, trois opérateurs, soit d 2 = 1,912). 

Rappelons que les moyennes d’échantillons extraits d’une même popu- 
lation mère de variance G 2 , sont distribuées avec une variance égale à 
°l = c 2 /(r-ra) (r ■ n correspondant ici à la taille des échantillons de 
mesure par opérateur). 

Si l’influence des opérateurs est réellement nulle, nous devrions avoir 


Deux cas peuvent se présenter : 

► si le premier membre de l’équation est inférieur ou égal au second, 
nous considérerons que l’influence des opérateurs est nulle : 

<*AV = 0 

» sinon, l’influence de l’opérateur sera caractérisée par la variance 
<J 2 AV-> calculée à partir de la formule suivante : 


R. 


°MoyOp ~ 


R- 



R x ; °ev 

^2 ylr-n 



2 



o 
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Dispersions du processus de contrôle 

La dispersion globale de l’équipement de mesure est donnée par : 


OqRR permet d’obtenir une estimation de la reproductibilité du pro- 
cessus de contrôle. 

G RR type « ANOVA » 

Ce cas correspond à une analyse de variances à deux facteurs avec répéti- 
tions d’expériences (traitée au chapitre 6). Elle permet de détecter des 
éventuelles interactions « pièces-opérateurs » et dans ce cas l’écart-type 
du processus <7 GRR de contrôle est légèrement modifié : 


Ce type d’étude donne une analyse plus fine et est préférable (exemple 
d’étude complète en figure 7.5). 

Avant de faire une analyse de variance, il est souhaitable de vérifier si les 
dispersions par opérateurs sont homogènes ou si, dans le cas contraire, 
elles mettent l’accent sur un problème d’utilisation de l’instrument (ins- 
trument bloqué par exemple). 




Figure 7.4 : Exemple d'interactions entre pièces et opérateurs 
(méthodes de mesure différentes) 



Accentuation du défaut de forme 
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On peut utiliser pour cela les tests de Hartley, de Cochran ou de Bartlett. 
Exemple : test de Cochran 

(sous réserve du même nombre de mesures par opérateur) 

Soit sf l'estimation de la variance pour l'opérateur i. On calcule le rap- 
port g = que l'on compare à la valeur tabulée gi_ a (a représentant 
le risque encouru). 

Cette valeur g Ua est fonction du nombre d'opérateurs et du nombre de 
degrés de liberté pour chaque estimation (v= n mesures - 7). 


I 

e 


o 
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Figure 7.5 : Exemple de fiche GRR de type « ANOVA w 1 
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Capabilité d'un processus de contrôle 

Les indicateurs sont calculés en comparant les dispersions de mesures 
soit avec la dispersion de production, soit avec la tolérance de la caracté- 
ristique étudiée. 


Conformité de l'appareil en contrôle de réception 

Il s’agit dans ce cas de vérifier la conformité du produit destiné au client 
(externe ou interne). Cet indicateur de capabilité maximum potentielle 
du moyen de mesure revient à comparer la dispersion de mesure ( k • <7rr) 
avec l’intervalle de tolérance IT. 

Soit par exemple, les règles d’usage qui suivent. 


Règle de l'édition 2 du MSA 

Nous pouvons calculer les quatre indicateurs suivants (PV sera définie 
au chapitre suivant) : 


%£K = ioo - 6 Gev 

IT 

%AV = 100- 6 ° AV 
IT 

%R&cR = 100- 6 ° GRR 
IT 

6 * o 

%PV = 100 ^ 


Tableau 7.3 : Règles « Édition 2 du MSA » 


%R&R < 10 % 

Équipement conforme 

1 0 % < %R&R < 30 % 

Acceptable (dérogation) si contraintes fortes (économique, 
disponibilité, limites techniques, etc.). Une démarche 
d'amélioration doit être entreprise. 

30 % < %R&R 

Problème de mesure (instruction de mesure à revoir, 
environnement, formation...). Sans amélioration potentielle, 
penser à changer de moyen. 


o 
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Règle type « Maîtrise statistique des procédés » 


C mc 


IT 


>4 


Cette règle est un compromis, elle conduit à retenir o m 


< IT/ 


Zone de conformité en contrôle de réception 

La zone de spécification doit être réduite si l’on veut garantir la qualité 
fournie au client (ISO 14253 - Partie l 1 ). Si cette réduction n’est pas 
faite, il y a partage des risques entre le client et le fournisseur. 

Figure 7.6 : Restriction de la zone de spécification 
pour garantir la qualité 



1. NF EN ISO 14253-1 : 1998, mars 1999, Spécification géométrique des produits 
(GPS) - Vérification par la mesure des pièces et des équipements de mesure. Partie 1 : 
règles de décision pour prouver la conformité ou la non-conformité à la spécification. 
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Conformité de l'appareil en validation de processus 

Industriellement, cette approche est moins utilisée. Dans la MSA 1 2 , il est 
demandé de comparer la dispersion de l’appareil avec la tolérance (%/7) 
et également avec la dispersion du moyen de production ( %TV, TV : 
dispersion totale ). L’édition 3 insiste sur ce point en introduisant un 
nouvel indicateur le « ndc 1 » ( Number of Distinct Categories ) défini par : 

ndc = J 2-^- 
°GRR 


Origine du « ndc » 

Partant de l’idée que la variation du processus de contrôle intervient 
dans la variabilité totale, on définit un facteur de discrimination 
DR : 


DR - 


DR - 




— %/ ndc 1 + 1 



2 , 2 
PV + °GRR 
2 

°GRR 


ndc = > !i.^L 

°GRR 

Suivant l’importance de la part de variabilité due au processus de 
contrôle, le ratio DR peut varier de un à l’infini. 


1. Automotive Industry Action Group (ALAG) (2002), Measurement Systems Analysis 
Reference Manual, 3 e éd., Chrysler, Ford, General Motors Supplier Quality Require- 
ments Task Force (www. cartvin. co. uk/qs) . 

2. Wheeler D. J., Lyday R.W., Evaluating the Measurement Process, 2 e éd., SPC Press Inc, 
1989. 


I 

e 


o 

© 


160 


Valider un processus de contrôle 


Figure 7.7 : Intégrer un processus de contrôle 
en fonction du « ndc » 



Choisir une machini 
capable de produire 
en conformité avec 
les exigences du client 
(indicateur Cp) 


Machine de production 
(Dispersion ô.apy) 


Choisir un appareil, 
non pas en fonction 
des données du client (/7) 
mais en fonction 
de la dispersion 
du moyen de productioi 
(indicateur ndc] 



I Appareillage de contrôle 
(Dispersion ô.Ogrr ) 


Si la dispersion de production n’est pas connue, alors par défaut, elle 
sera estimée à partir du lot de pièces mesurées. Chaque pièce étant 
mesurée k ■ r fois ( k opérateurs faisant chacun r mesures), on estime le 
mesurande (dimension de chaque pièce) par la moyenne des mesures. 
Pour cela, il faut s’assurer que l’échantillon de pièces soit représentatif 
(image de la dispersion du moyen de production et 50 % environ des 
pièces en dehors du « bruit de mesure ») (fig. 7.8). L’idéal est d’estimer la 
dispersion de production par une approche de type MSP. 

Figure 7.8 : Vérification de la représentativité 
de l'échantillon de pièces 
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À partir des n moyennes, on détermine l’étendue Rpv À l’aide de ce 
résultat, on calcule un écart-type de production O py {Part to part Varia- 
tion) (dans l’exemple traité, il y a un sous-groupe de dix mesures). 



Partant de l’additivité des variances, on calcule la variation totale appa- 
rente du processus, incluant les dispersions de production et de mesure 
{Total Variation) : 


°TV ~ + {°Pv) 

Attention, si l’on a des données provenant d’une démarche MSP, cela 
correspond probablement & (Jjy, car il y a de grandes chances que les 
mesures aient été faites avec cet appareil. On en déduit Opye t l’on cal- 
cule le nombre de catégories distinctes « ndc » qui, dans l’idéal, doit être 
supérieur ou égal à cinq : 

°GRR J 

Conformité de l'appareil en production dans le cas 
d'une capabilité du processus élevée 

La démarche avec l’indice ndc vu précédemment est bien adaptée 
lorsque le processus de production est capable sans excès. Il peut arriver, 
suite à des problèmes de disponibilité et de charge machine, que l’on 
accepte de faire de la sur qualité. L’indicateur Cp est alors élevé (par 
exemple Cp = 4). Il est tout à fait possible, dans ce cas de figure, de 
piloter une production avec un appareil dispersant dans le même ordre 
de grandeur que le moyen de production, et cela sans risque pour le 
client. 

I 
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Figure 7.9 : Cas d'un processus de production 
correspondant à un indicateur Cp élevé 
(possibilité d'une dispersion de mesure plus importante) 



Cas des contrôles automatiques 

De plus en plus, on voit apparaître en production de grande série des 
postes automatiques de contrôles finaux. C’est souvent une exigence du 
client (donneur d’ordre), rassuré par la mesure systématique des pro- 
duits achetés. La mesure de toutes les caractéristiques de définition du 
produit exige des solutions technologiques d’acquisition très rapide 
(inférieure ou égale aux temps de production). Reste à savoir si la 
cadence élevée de contrôle ne détériore pas la qualité métrologique. 
Tout étant automatisé, l’influence de l’opérateur n’existe plus. Pourtant, 
comme la culture GRR est présente dans l’entreprise, on aimerait bien 
avoir « quelque chose d’équivalent » pouvant servir d’indicateur de qua- 
lité de mesure. Rappelons qu’une étude GRR n’est qu’une analyse de la 
variance à deux facteurs d’influence (pièces et opérateurs) avec répéti- 
tion de la mesure (voir chapitre 6). On peut très bien faire le même type 
d’étude, où le facteur opérateur est remplacé par un facteur présumé 
g influent (prise de pièce, indexation, milieu ambiant. . .). Une enquête de 

1 terrain, auprès des utilisateurs journaliers, est souvent riche d’enseigne- 

ments pour détecter les facteurs d’influence. 

ô 
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Concept de justesse 


L’étude précédente a permis de mettre en évidence l’influence de diffé- 
rents facteurs. Par exemple, si un des opérateurs donne des résultats plus 
forts que ses collègues, cela se remarquera. Par contre, si l’ensemble des 
opérateurs dérivent systématiquement de façon identique, cela ne se 
remarquera pas. En conséquence, il est nécessaire d’avoir une valeur 
vraie matérialisée par un étalon. La valeur vraie est par définition 
inconnue. Une mesure très fine de l’étalon (correction incluse) permet 
d’approcher cette dernière, on parlera d’une valeur conventionnellement 
vraie. C’est à partir de cette donnée que l’on testera l’appareil. 

Une valeur étalon ne correspond pas forcément à une valeur nominale 
(par exemple, 50,000 mm). Physiquement, cela peut être une pièce, 
extraite de la production, qui a été mesurée par des moyens plus perfor- 
mants. L’estimation de l’écart moyen entre un échantillon de mesures 
de l’étalon et une valeur supposée vraie de l’étalon permet d’estimer un 
écart de justesse locale ( biais ). C’est pourquoi, si possible, il faudra 
retenir un étalon voisin des pièces à mesurer car bien évidemment 
l’erreur de justesse peut varier en fonction de l’entrée (défaut de linéarité 
et d’échelle). 

L’erreur de justesse est définie par : 

A J = fi-x Q 

» H correspond à la moyenne d’un nombre infini de mesures ; 

» Xq est la valeur de l’étalon (supposée vraie). 

En réalité, nous ne pouvons qu’estimer l’erreur de justesse par : 


ôj = x-x 0 

» x est l’estimation de [1 à l’aide des n valeurs de l’échantillon (l’incer- 
titude sur cette estimation est de type A, voir chapitre 8) ; 

► Xq est la valeur donnée par le certificat d’étalonnage (estimation de 
*o)- 


La tolérance d’étalonnage permet de calculer l’incertitude associée à la 
valeur étalon (incertitude de type B). 

L’incertitude sur AJ est caractérisée par : 
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Par exemple, si le certificat d’étalonnage propose : X 0 - x 0 ± T 0 , on peut 
retenir avec l’hypothèse d’une distribution uniforme : 


2 Mesure , 0 

S “ J ~ n + 3 

La qualité de l’estimation peut être bornée en utilisant la loi de Student 
(oc représentant le risque, voir exemple figure 7.10) : 


Probabilité [<5/ - f„_ 1;1 _ a/2 • s AJ < AJ < 8 J + t„_ hâ _ al2 • s AJ ] = 1 - a 


Figure 7.10 : Illustration de l'intervalle de confiance 
de l'erreur de justesse 



Note : tant que la valeur « x 0 » est incluse dans cet intervalle de con- 
fiance, il est possible que le biais mesuré ne soit que l'expression de la 
variabilité due à l’échantillonnage. 


Bornes d'acceptation pour l'erreur de justesse 

Rappelons que l’on espère avoir un appareil avec une erreur de justesse 
tendant vers 0 (espérance nulle). Plusieurs règles peuvent être envisa- 


» si possible la valeur « 0 » doit rester fortement probable, c’est-à-dire 
qu’elle doit appartenir à l’intervalle de confiance (exemple figure 7. 1 1). 
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» si ce n’est pas le cas, cela peut éventuellement provenir également de 
l’estimation de l’écart-type sur AJ (valeur très faible due à des mesures 
pratiquement identiques), qui conduit à un fort resserrement des 
bornes de l’intervalle de confiance. Si cette hypothèse se confirme, on 
peut comparer AJ avec la tolérance et exiger par exemple (norme 
interne d’usage) : 


Dans le cas contraire, on peut conclure à une erreur systématique de jus- 
tesse effective et envisager une retouche ou un rebut de l’appareil. Il est 
possible de faire de la bonne mesure en corrigeant l’erreur de justesse, 
mais cela est dangereux en production. Cette approche est plus réservée 
à une mesure de laboratoire. 

Rappelons qu’il s’agit ici d’une mesure de justesse locale. Dans le cas 
d’appareil de type universel à forte étendue de mesure, il faudra vérifier 
l’erreur de justesse en plusieurs points de cette étendue (chapitre 6 et 8). 



2 


O 

© 


166 


Valider un processus de contrôle 


Figure 7.11 : Exemple d'estimation d'erreur de justesse 


ETUDE DE LA JUSTESSE LOCALE A J 

Référence de l’appareil : | 

CIMG365 

lOpérateur : |Jean 

Référence pièce :| 

MC014 

23/05/2008 19:31 


Valeur de l’étalon de référence | 

3,000 iRéférence étalon 

E/M3CSL 

Tolérance d’étalonnage à ± 

0,005 Uniforme 

▼ Estimation o Et : 0,003 


</> 

LU 

CL 

=> 

CO 

UJ 

2 

2.980 

2.990 

2,990 

2,990 

3,000 

3,010 

3,000 

3,020 

3,030 

3,030 

3,000 

3,010 

3,020 

2,980 

3,010 










































| Moyenne des mesures X : 


Ecart-type estimé (mesures) : 


| |Ecart-type estimé : 


||\J estimée Moyenne des mesures - Valeur de l’étalon = 6 | 

0,0040 lie 

||Ecart-type estimé sur Aj : 

0,0051 H 


ENCADREMENT DE AJ 
(Loi de Student) 
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Capabilité du processus de contrôle : 
approche inertielle 

Dans la validation des processus de contrôle, l’approche classique 
consiste, soit à valider de façon séparée la justesse et la fidélité, soit 
simplement à additionner les variances de ces deux éléments. La propo- 
sition du ndci est une approche permettant de donner plus de cohérence 
aux deux études de fidélité et de justesse en se fixant une limite qui 
dépend de la combinaison de ces deux éléments (illustration en 
figure 7.12). 


Rappel sur les contraintes relatives aux dispersions 


Si pour le choix du moyen de production, on retient un indicateur Cp 
qui doit être supérieur à 4/3, cela donne : 


soit : 


Cp = 


fJçUeM_ > ^ 

6 • (J py 3 


Si l’on retient un indicateur « ndc », supérieur ou égal à 5, cela donne : 

Æi 


En combinant ces inéquations, nous obtenons une valeur seuil pouvant 
être utilisée directement en contrôle de réception : 

0 <dL.J T ITclient 

U CRR ~ 11 Client 
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Influence du décentrage 

La qualité livrée au client est caractérisée par sa dispersion Opy et par un 
décentrage 8p. En réalité, l’estimation de la qualité intrinsèque est 
perturbée par le processus de contrôle entraînant une dispersion et un 
décalage parasites. Cette identification globale de la qualité est quanti- 
fiée par : 


Nous ne chercherons pas à limiter ces deux caractéristiques séparément 
(dispersion et décentrage) mais globalement. Le concept de tolérance 
inertielle nous permettra une analyse permettant des compensations 
entre ces deux composantes. 

Tolérancement inertiel 

Le tolérancement inertiel 1 a comme principe de vérifier qu’une sortie Y 
(pour un lot de k produits) présente une dispersion acceptable en regard 
de l’écart systématique constaté par rapport à la cible. Ces deux carac- 
tères (dispersion résiduelle et écart par rapport à la cible) peuvent être 
représentés comme les composantes d’un pseudo-vecteur : 



Le seuil de conformité sur la sortie Y permet des compensations entre 
ces deux caractères (sous forme de lien). Par exemple, une faible disper- 
sion permet de compenser un écart systématique par rapport à la cible 


Le lien entre ces deux caractères correspond à la norme | || 2 appliquée 
à/. 


1. Pillet M., « Inertial Tolerancing », The Total Quality Management Magazine, vol. 16, 
n° 3, mai 2004, p. 202-209. 
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Nous appellerons cette norme « inertie » et la noterons /. La décision de 
conformité pour la sortie Y sera donnée si : 


Sauf mention particulière du client, l’inertie maximum sera déduite de 
la tolérance classique par : 


Dans le cas contraire, cela suppose que le client a attribué à chaque 
composante de définition une inertie spécifique que devra respecter le 
fabricant. 

Indice de capabilité Cpi 

Par analogie avec l’indicateur Cp, il est possible de définir un nouvel 
indicateur : 



Cpi ~ — 



En réalité, l’indicateur évalué est égal à : Cpi ohs = 
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Figure 7.12 : Indicateurs de capabilité inertielle 


Capabilité inertielle Cpci 


Tolérance : 0,400 

Inertie l M „ : 0,067 


Justesse 5 0 : 

0,016 1 |lnertie l Q : ] 0,016 T 

0,004 | 

Indice de capabilité Cpci mini : 

5 1 

Indice de capabilité Cpcin M (l„„/l G ) 

Seuil d'l 0 pour garantir Cpci mini : 

0,013 

lertie l 0 : 


Inertie l Q / limite client & limite production 



:ie processus de contrôle IG 


Inertie l P : | 0,033 | 

Indice de capabilité prod. Cpi (WI P ) : 

2,01 

? 

1 

f 

1 

Indice de capabilité CpcipRQQ (lp/l G ): 


Seuil d'l G pour un ndc,„, M limite | 0,009 

Inertie l 0 : 
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Couplage inertie du processus de contrôle 
et inertie du processus de production 

L’inertie résultante vectorielle I T est la combinaison de l’inertie vecto- 
rielle de production I p et de l’inertie vectorielle de contrôle I G . 

Loi inter ne 

L’ensemble des inerties vectorielles est muni d’une loi de composition 
interne, qui appliquée aux inerties vectorielles précédentes, donne : 


L’inertie résultante, combinaison des deux processus (production et 
mesure), est égale à : 


Influence des décalages 

L’équation précédente ne peut pas se réduire à la loi classique : 


car les sources de variabilité n’étant qu’au nombre de deux, il n’y a pas 
d’effet de compensation a priori. 

* 5q sera estimé par une analyse de la justesse du processus de contrôle. 
Cela nécessite l’emploi d’un étalon voisin des pièces de la produc- 
tion. 

* ôp sera estimé indirectement à partir des observations d’un journal 
de bord d’un pilotage par SPC nous donnant le décalage résultant ôp 
instantané. 

Nous obtenons la relation : 



soit au niveau des composantes : 
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« ndc » inertiel 

Suivant l’importance de la part de variabilité due au processus de 
contrôle par rapport à la variabilité totale, nous avions défini un indica- 
teur « ndc » à partir d’un ratio DR. Par analogie, ce calcul de ratio peut 
s’appliquer aux inerties de ces deux processus. Il peut varier, de même, 
de un à l’infini : 




ndc ! - yfl ■ 


f 2 • 8 p ■ 8 G 


dl r + 2 -8 p S r 


p ' Ip + Iç + 2 • S p ■ Ô G 
I 2 g + 2S p S g 


^ndc 2 j 


+ 1 


Comme Sp n’est pas directement mesuré, nous retiendrons la formule 
pratique suivante : 


ndc r = s[2- 


\JIg + 2-(S t <5 g )-<5 g 


Remarquons que : 


ndc ; - 


1 c; 


Partant du principe qu’un pilotage permet de faire tendre vers zéro Sp 
nous poserons que : 


sf2-—> ndc t limite 
Ir. 


En résumé, avec l’hypothèse d’un Sp tendant vers 0 (c’est-à-dire ôj 
= 8q}, il faut choisir un processus de contrôle dont la valeur Iç, relative- 
ment à l’inertie de production Ip, respecte la règle suivante : 


i e s ^ ■/, 

ndcjkmite 
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» Iç est déterminé suite à l’analyse de la « répétabilité-reproductibilité » 
du processus de contrôle qui nous permet de calculer Gçpp [MSA 
2002\ et d’une étude de justesse donnant l’estimation ôç. 

» Ip est obtenu par comparaison de l’inertie globale Ip, déterminée à 
partir de l’analyse MSP et de l’inertie du processus de contrôle Iç. 

Si l’on veut comparer directement Iç à Ip, nous obtenons : 




^[ndcj limite) 2 + 2 


soit par exemple pour un ndc\ limite de 5, cela donne sensiblement la 
contrainte suivante : 


I r <-I T 

g 4 r 


Figure 7.1 3 : Illustration graphique avec exemples 
des inerties l P et l c 
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Chapitre 8 


Estimer l'incertitude 
de mesure 



Rappelons qu’une mesure sert à prendre une décision. Cela peut servir : 

► à décréter la conformité d’un produit (par rapport aux exigences du 
client) ; 

► à valider un processus (capabilité, performance, pilotage, sécurité, 
développement durable, etc.). 

Cette mesure a comme objectif d’être le reflet de la qualité réalisée. Mal- 
heureusement, le processus de mesure va donner une image perturbée 
(ou bruitée) de la qualité réellement produite. Pour prendre une déci- 
sion à bon escient, il faut que l’on soit conscient de cette perturbation. 
Pour cela, on cherchera à définir une zone liée à la mesure réalisée, où la 
valeur vraie (jamais connue par définition) à de grandes chances d’être 
située. C’est le concept même d’incertitude. 

Le métrologue doit gérer son processus de contrôle comme un fournis- 
seur (de prestation). Il doit situer sa performance par rapport aux exi- 
gences de mesure de son client. On retrouve la notion classique 
d’indicateur de capabilité (voir chapitre 7). 

Nous nous contenterons dans ce chapitre d’étudier des incertitudes sur 
des mesures simples nécessitant des corrections d’influence (influence 
de la température par exemple) ou sur des mesures indirectes à partir de 
grandeur de base (comme la mesure d’une masse volumique). 

Pour estimer cette incertitude, nous nous appuierons sur une approche 
statistique de la mesure. La mesure résultante sera modélisée suivant la 
démarche proposée par le GUM ( Guide to the Expression of Uncertainty 
in Measuremeni) . Cette démarche, initialisée par le Bureau international 
des poids et mesures et pouvant être appliquée à de multiples secteurs 
d’activités, a été retenue par de nombreuses instances nationales de 
métrologie. Ce guide est proposé comme norme française sous la 
référence : NF ENV 13005 août 1999, Guide pour l’expression de 
l’incertitude de mesure - Indice de classement : X07-020. 

Pour approfondir l’approche statistique de la mesure, on pourra se 
reporter entre autres, aux normes suivantes : 
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NF ISO 3534-1 Janvier 2007 - Statistique - Vocabulaire et sym- 
boles - Partie 1 : Termes statistiques généraux et termes utilisés en 
calcul des probabilités. 

ISO/TS 21748 Mars 2004 - Lignes directrices relatives à l’utilisa- 
tion d’estimations de la répétabilité, de la reproductibilité et de 
justesse dans l’évaluation de l’incertitude de mesure. 

NF ISO 5725-1 Décembre 1994 - Application de la statistique - 
Exactitude (justesse et fidélité) des résultats et méthodes de 
mesure - Partie 1 : Principes généraux et définitions. 

NF ISO 5725-2 Décembre 1994 - Application de la statistique - 
Exactitude (justesse et fidélité) des résultats et méthodes de mesure 
- Partie 2 : Méthode de base pour la détermination de la répétabilité 
et de la reproductibilité d’une méthode de mesure normalisée. 

NF X07-001 Décembre 1994 - Normes fondamentales - Vocabu- 
laire international des termes fondamentaux et généraux de métrologie. 
FD X07-021 Octobre 1999 - Normes fondamentales - Métro- 
logie et applications de la statistique - Aide à la démarche pour 
l’estimation et l’utilisation de l’incertitude des mesures et des résul- 
tats d’essais. 

FD X07-022 Décembre 2004 - Métrologie et applications de la 
statistique - Utilisation des incertitudes de mesures : présentation 
de quelques cas et pratiques usuelles. 


Présentation de la méthode du GUM 

Exemple conducteur 

Pour donner un aperçu de la méthode, nous retiendrons comme 
exemple conducteur la mesure indirecte de la masse volumique d’un 
cylindre. Partant de grandeurs de base (longueurs et masse), avec des 
hypothèses simplificatrices (en considérant le volume comme un 
cylindre de révolution parfait), nous allons calculer la masse volumique 
p (modélisation de la sortie en fonction des entrées). Nous rappellerons 
les formules génériques et leurs applications à cet exemple (fig. 8.1). 
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Figure 8.1 : Exemple conducteur : mesure d'une masse volumique 



Modélisation du processus de mesure 


C’est un des points les plus délicats. Il faut recenser les variables de base 
( X ù permettant de calculer la variable finale (Y) et proposer un modèle 
(en essayant d’être le plus exhaustif possible). On pourra s’aider d’une 
démarche 5M ou de la liste publiée dans le GUM. 

» Écriture générique : Y - f{X v X 2 , X n ) 


M M 



p-—-M- D~ 2 ■ h~ l 

K 


Variables de bases cachées 

Dans l’exemple précédent, avec la formule retenue, nous n’avons pas pris 
en compte l’influence de la température. Si son influence n’est pas direc- 
tement modélisée dans la formule donnant p, de fait elle sera implicite- 
ment prise en compte dans les paramètres dimensionnels {D et h) sous 
réserve d’être en conditions de reproductibilité. Ce type de variable 
engendre de la covariance entre les variables identifiées. 

Note sur les constantes mathématiques : n est une constante mathé- 
matique (a priori invariable). Pour calculer les incertitudes, nous ferons 
de l'analyse numérique où K ne sera pas exactement connue (simplement 
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estimée par une valeur décimale). En toute rigueur, cette définition 
incomplète engendre de l'incertitude. Nous négligerons, dans les exem- 
ples suivants, cette part d'incertitude (ordre considéré comme supérieur). 

Variables de bases ignorées 

Il s’agit de variables ayant une influence sur le résultat final de mesure, 
mais qui n’interviennent pas implicitement dans les autres variables de 
bases recensées. Bien qu’influentes, elles peuvent être complètement 
ignorées du métrologue et fausser son estimation. Nous en donnerons 
un exemple simple en mesurant la surface d’une tôle rectangulaire à 
l’aide des dimensions de base ( a et b), en ignorant le défaut de perpendi- 
cularité a (fig. 8.2). 

Figure 8.2 : Exemple de variable ignorée : a 


Le modèle simple donne : S — a- b. 

alors que nous avons en réalité (a est infiniment petit du premier 
ordre) : 



cos a 



La mesure du même mesurande, par des méthodes différentes, peut per- 
mettre de « démasquer » ce genre de variables. 
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Loi de propagation de l'incertitude 

Le modèle Y — f{X v X 2 , X n ) permet d’obtenir une estimation y 

de Tà partir des estimations de base (entrées) x\, x 2 . ...x n : 


y = f(x v x 2 , x n ) 

La qualité finale de l’estimation y va dépendre de la variabilité rencon- 
trée sur les estimations Xj. 

Le paramètre caractérisant la variabilité dans l’estimation de X^ est 
appelé incertitude type, il sera noté u(x t ). Mathématiquement, il repré- 
sente l’écart-type <J 2 - (la lettre u, comme « uncertainty », permet de se 
démarquer de la dispersion de production). 

Le GUM retient deux approches pour quantifier les incertitudes types 
de base (évaluation de type A ou de type B). Nous reviendrons plus en 
détail sur ces deux types d’évaluation. 

L’incertitude composée est donnée par l’équation suivante (obtenue à 
partir du développement de / autour des moyennes des x z - en série de 
Taylor du premier ordre) : 


u c oo=X 

i=\ 

’dffi 

dx i 

2 

■u\x t ) 

+2 I É 

Ésà j=i+V 



coefficient de sensibilité pis rapport à la variable x t , se lit « d rond/ 
dx i 

sur d rond x i ». Il correspond à la dérivée partielle de /par rapport à la 
variable Xj au voisinage de la valeur Xj. Cette dérivée se calcule comme 
pour une dérivation classique, il suffit de considérer les autres varia- 
bles comme des constantes pendant l’opération de dérivation ; 
u(x t ) est l’incertitude type de la variable Xj (écart-type), u 2 {x-) corres- 
pond à la variance ; 


u(xpxj) = u(xj,x ■) est la covariance des variables x;,xj (attention aux 
unités). 
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Il est souvent plus pratique de présenter l’équation précédente sous la 
forme suivante : 


f 

*£00= É 

y=i l 

'd£ 

dx i 

W 

il 

*=1 j=i+V 



Rappels 

u(xi,Xj) = Cov {xj,xj) = E {xjyj) - E(^).E (xj) = x { x- — x i • x- 
avec E(x z ) : espérance mathématique de la variable aléatoire Xj. 
u 2 (xj) = u (Xp xj) = Var (x) = E (x 2 ) - (E(x)) 2 = xj ~[x i j 


Coefficient de corrélation : r(xj,Xj) = u(x;,x-)lu(x-) u(xj = 


Cov(x i ,x •) 


Lorsque les variables ne sont pas corrélées (r = 0), l’incertitude com- 
posée se calcule comme suit : 


?(r) = Z 


dX: 


' (X:) 


Cette formulation est beaucoup plus simple. Il ne faut l’utiliser que 
lorsque l’on a vérifié que les covariances sont nulles (ou au moins que 
l’on a de fortes présomptions). 

Application au calcul de la masse volumique p : 

+2 ■ ^r~ ~ • r{M, D) ■ u(M) ■ u(D ) 
oM oD 

+2 ■ • r(M,h) ■ u(M) ■ u(h) 

dM dh 

+2-^ ~-r(h,D)-u(h)u(D) 

dh dD 
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Si l’on fait l’hypothèse que la masse (choix de la matière) est complète- 
ment indépendante du volume, on peut admettre que r(M,D ) et r(M,h) 
sont nuis. L’expression précédente se simplifie fortement : 




+2-^--^--r{h,D)-u{h)-u{D) 
dh ôD 


Exemple 

Pour un cylindre en alliage d'aluminium, on a fait cinquante mesures de 
M, D et h. Les triplets de mesures (M, D et h) sont mesurés successive- 
ment sur un temps relativement court (répétition) à 20 °C. On a estimé 
ensuite, avec une variation de température de l'ordre de 18 °C à 26 °C, 
les mesures corrigées que l'on aurait eues dans des conditions de repro- 
ductibilité (tableau 8.1). 


Tableau 8.1 : Mesure des paramètres avec et sans influence 
de la température 



Masse M 

D 

(0 = 20 °C) 

h 

(0= 20 °C) 

P 

(0=20 °C) 

D 

(0* 20 °C) 

h 

[6*20 °C) 

P 

(0*20 °C) 

(g/mm 3 ) 

992.1656 

80.0083 

70.0088 

0.0028 

80.0117 

70.0118 

0.0028 

S 

(g/mm 3 ) 

0.0975 

0.0049 

0.0051 

4.81 E-07 

0.0065 

0.0048 

5.84E-07 


Si l'on mesure à température constante, l'influence de la température 
sera négligée. Elle n'introduira pas de covariance entre D et h (fig. 8.3). 
Nous retiendrons la formule simplifiée suivante : 



À l'aide des valeurs moyennes des variables, on obtient : 

u 2 c (p) = 8.072E-1 2 • 9.51 0E-03 + 4.965E-09 • 2.390E-05 + 1 .621 E-09 • 2.61 7E-05 

u 2 (p) = 2.378E-13 

u c (p) = 4.88E-07g/mm 3 


183 


Qualité de la mesure en production 


Figure 8.3 : Distribution des couples (D et h) en fonction du temps 
à température constante de 20 °C 


Diamètre/Hauteur 



Si l’on fait le même type d’étude, avec variation de la température, sans 
tenir compte de la covariance entre D et h, nous obtenons : 

u c {p) = 8.070E-12 • 9.5 10E-03 + 4.963E-09 • 4. 195E-05 + 1 .621E-09 • 2.276E-05 
« 2 (p) = 3.218E-13 
u c (p) = 5.67E-07g/mm 3 

Or, à l’examen de la répartition des couples ( D et h) mesurés à des tem- 
pératures différentes, on subodore une légère corrélation (fig. 8.4). 
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Figure 8.4 : Distribution des couples (D et h) en fonction du temps 


avec variation de la température 
Diamètre/Hauteur (influence de la température) 



En tenant compte de la covariance, nous obtenons : 



soit pour cet exemple : 

u c (p) = 8.072E-12* 9.510E-03 + 4.965E-09 • 2.390E-05 + 1.621E-09 • 2.617E-05 


+5.672E-09 • 4.866E-06 
U'(p) = 3.494E-13 
u c (p) - 5.91E-07g/mm 3 

Note : il est parfois plus simple de traiter la covariance que de définir 
une relation fonctionnelle plus complexe (par exemple, retenir comme 
variables supplémentaires la température et le coefficient de dilatation 
du matériau). 
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Risque dû au développement limité au premier ordre 

On peut remarquer de légères différences entre les valeurs calculées et les 
valeurs estimées par échantillonnage. Ceci est particulièrement vrai 
lorsque la réponse Y n’est pas une fonction linéaire des entrées X i 
(fig. 8.5). Dans un cas de forte non-linéarité, il faut envisager un déve- 
loppement de Taylor au deuxième ordre. De toute évidence, si la varia- 
bilité sur les entrées est forte, il faut rester prudent. 


Figure 8.5 : Influence de la non-linéarité sur l'estimation 
de l'incertitude composée 



Quantifier les incertitudes types u(xj) : 
évaluation de type A ou de type B 

Le GUM retient deux modes d’évaluation des incertitudes types u(x-), à 
savoir le type A et le type B. 

Évaluation de type A 

Ce type est basé sur l’utilisation des techniques statistiques (variance expé- 
rimentale, variance de distribution des moyennes, etc.) (voir chapitre 6). 
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Pour un échantillon de n résultats de mesure, la moyenne arithmétique 
est le meilleur estimateur d’un nombre infini de mesures. 


L’estimation de la moyenne \l l du paramètre X i est donnée par : 

x- - X - -Y y . C’est la valeur à prendre en compte dans la formule 

y — f(x 1 ,x 2 , x n ) (voir exemple précédent sur la mesure de masse 

volumique). 

L’estimation de l’écart-type expérimental est donnée par : 



Si l’on dispose de plusieurs séries de mesures pour un même paramètre, 
avec pour chaque série une estimation s lk , une meilleure estimation de la 
variance - sous réserve d’appartenance à une même population - sera 
donnée par : 

2 _ ( n i ~ 1)4 + ( n 2 ~ 1)4 +—+(n k ~ 1)4 
(«i-l) + («2 -!) + ••• + («*-!) 


L’incertitude type sur la moyenne X de l’échantillon est obtenue par : 

2 

u 2 (x ) = — (sous réserve d’indépendance des mesures x-j) . 

Ce type d’évaluation est souvent préféré au type B (parfois à tort). Il faut 
rappeler que répéter la mesure permet de faire chuter l’incertitude type 
(au détriment du coût et du temps d’acquisition), que cette diminution 
varie en inverse de V«, c’est-à-dire que la chute est rapide uniquement 
au début (fig. 8.6). 

s i est une pauvre estimation de il faut être assez prudent sur la perti- 
nence de la valeur obtenue (chapitre 6). 


soit : 



Prob [M?r”>* s / ~ a t - V/2- s), ' s .] = 1 -° 
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Pour une probabilité de 95 %, les bornes k^e . t k$ u p sont rej 
sur la figure 8.7 (en fonction de la taille de l’échantillon). 

Figure 8.6 : Chute de l'incertitude type en fonction de n 
Évolution de 1 /Racine (n) 


^présentées 


Figure 8.7 : Fourchette (coefficients) de l'estimation de a,- 
Fourchette d'estimation à 95 % 


♦ klnf 
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Nb mesure 


Évaluation de type B 

Ces méthodes sont utilisées, soit lorsque l’outil statistique est d’utilisa- 
tion difficile (trop onéreux ou destructif par exemple), soit parce que 
l’on pense bénéficier de meilleures estimations, provenant d’expertises 
internes ou externes. Ainsi, si l’on a une garantie sur la variable X- (par 
exemple, certificat d’étalonnage donnant l’incertitude élargie), on peut 
estimer l’incertitude type (sous hypothèses de distribution) comme on le 
ferait à partir d’un intervalle de confiance. 

Cette démarche nécessite une analyse de praticien. Elle peut faire appel 
à une expertise construite sur des observations (mémoire). Parfois, c’est 
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une information donnée par un tiers qui sera utilisée (certificat d’étalon- 
nage, données constructeurs, valeurs encyclopédiques, etc.). Souvent, 
on se contente des lois de probabilité suivantes, représentées sur l’inter- 
valle ± a (figure 8.8 graduée en a! 3) : 

► loi uniforme (ou loi rectangulaire) : 

f( x ) = J- u (x) = ^=. 

2-a V 3 

Cette loi est retenue lorsque l’on a peu d’information. Par exemple 
un certificat d’étalonnage nous garantit que la valeur X ) appartient à 
± a, mais ne nous donne aucune information complémentaire sur la 
densité de probabilité f{x). Par défaut, toutes les valeurs internes 
seront considérées comme équiprobables ; 

► loi normale (ou loi de Gauss) : 


i -fc# 

/(*) = — f=e 2(7 u{x) = -, 

Oylln k 

Toutes les valeurs n’ont pas une probabilité de 100 % d’être com- 
prises dans ± a, mais une probabilité qui dépend du choix de niveau 
de confiance (1 — a) retenu (par exemple si k = 3, cela veut dire que 
la probabilité d’appartenance à ± a est de 99,73 %. Si k = 2, elle est 
voisine de 95 %). Cette loi est d’un usage fréquent, elle est bien 
appropriée lorsque les causes de variabilité sont nombreuses et du 
même ordre de grandeur (théorème central limite). 

Il existe d’autres lois pouvant être retenues. Ce choix peut être guidé par 
l’analyse d’histogrammes avec l’appui de tests statistiques. Nous citerons 
simplement : 


* loi triangulaire avec u(x) — 




comparée à la loi uniforme, on estime 


que la position moyenne est la plus probable. C’est un compromis 
des deux précédentes lois ; 

* loi en U (ou dérivée d'arc sinus 1 ) : 



«(*)=- 7 *. 

V2 


'j. metgen.org. 
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Cette loi est bien adaptée à une grandeur centrée sur b et pouvant 
évoluer entre ± a de façon sinusoïdale. On rencontre ce problème 
avec des variables asservies et présentant un problème de pompage 
(climatisation par exemple). 

Figure 8.8 : Comparaison des quatre lois usuelles 



Incertitude élargie 

En statistique, on donne pour une grandeur étudiée, un intervalle de 
confiance ayant une certaine probabilité de contenir la valeur vraie 
(chapitre 6). Pour construire cet intervalle, on utilise généralement : 

► une estimation donnée par la moyenne de n observations indépen- 
dantes ; 

un paramètre de variabilité (par exemple l’estimation de l’écart-type) ; 

► une fonction de densité de probabilité f{x) (par exemple, la loi nor- 
male). 

Nous n’avons pas cette dernière information. Cependant, il est pratique 
de donner une mesure sous la forme : « 

1 

| Mesure = Valeur estimée ± U | & 

o 

U est appelée « incertitude élargie » et calculée par : U - k -u c {y). 
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Le passage de u c à U ne donne pas de nouvelle information, car le coeffi- 
cient k, appelé « facteur d’élargissement » est choisi arbitrairement 1 . Le 
plus souvent, par analogie avec le calcul de probabilité à l’aide de la loi 
normale, on retient : 2<k <3- 

Les indicateurs de capabilité d’un processus de production sont calculés 
à l’aide d’une dispersion définie à 6.(7. Par souci de normalisation, sauf 
indications contraires, il paraît normal d’étendre cela à la mesure et de 
retenir k = 3. 

Expression de l’incertitude : on peut soit donner l’estimation du 
mesurande et l’incertitude composée associée, soit donner l’estimation 
et l’incertitude élargie. Nous retiendrons la seconde solution, qui est 
plus dans l’esprit de la production. Il faut prendre garde de retenir la 
même résolution pour ces deux composantes. 

Exemple 

M = 99,166 g ±0,095 g (k= 3) 


Incertitude sur les corrections 

Nous distinguerons deux types de correction : 

* une correction locale pour un processus de mesure utilisé en série avec 
des mesurandes voisins (par exemple un lot de pièces usinées sur la 
même machine) ; 

* une correction globale, applicable sur toute l’étendue de mesure de 
l’appareil (appareil d’emploi plus universel, comme un pied à cou- 
lisse). 


Pour un calcul plus fin de k, le lecteur pourra consulter l’annexe G, « Degrés de liberté 
et niveau de confiance » de la norme NF ENV 13005 août 1999, Guide pour l’expres- 
sion de l’incertitude de mesure. 
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Correction locale 

À partir d’une estimation (mesure brute), nous appliquerons diverses 
corrections pour obtenir une mesure corrigée (meilleure image du 
mesurande) : 


Mesure corrigée = mesure brute + correction de justesse + correction 
de dérive + correction des facteurs d’influence + autres corrections. . . 


Le plus souvent, ces corrections (supposées non corrélées) ne sont pas 
faites pour différentes raisons (ignorance des influences, complexité de 
modélisation, coût, qualification, etc.). Ne pas faire ces corrections ne 
veut pas dire qu’elles ne sont pas prises en compte, mais qu’elles ont été 
estimées dans un premier temps à la valeur « 0 ». 

Mathématiquement, sous réserve d’indépendance des différentes cor- 
rections, on peut écrire : 

M„ r = M^'Lc t 

»^>=i'»L_ + X< 1 '“ 2 w> 


Correction locale de justesse 

Pour faire une correction de justesse, il faut avoir une estimation de la 
différence entre la mesure et une valeur vraie. 

Rappels (chapitre 7) : 

L'erreur de justesse est définie par : 

AJ = p-X 0 

• n correspond à la moyenne d’un nombre infini de mesures ; 

• X 0 est la valeur de l'étalon (supposée vraie). 

En réalité, nous ne pouvons qu'estimer l'erreur de justesse par : 

I 

8J = x-x 0 | 

• x est l'estimation de p à l'aide des n valeurs de l'échantillon (l'incerti- §- 

tude associée est de type A) ; J 

© 

• x 0 est la valeur donnée par le certificat d'étalonnage. 
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La tolérance d'étalonnage permet de calculer l'incertitude associée à la 
valeur étalon (incertitude de type B). 


L'incertitude sur AJ est caractérisée par : 



T L/on 


et estimée par : 


„2 _ ^ MesureEta/on 

AJ ~ n 



Par exemple si le certificat d’étalonnage propose : X Q — X Q ± T Q , on 
peut retenir avec l'hypothèse d'une distribution uniforme : 



3 


L'erreur de justesse est estimée à SJ 

Cela ne nécessite pas une courbe complète d’étalonnage, mais simple- 
ment une connaissance locale pour le voisinage de la valeur nominale du 
mesurande. 


L'erreur de justesse n'est pas connue 

Dans ce cas, on fera de fait une correction nulle. Ce choix peut aussi être 
dicté par des considérations économiques et de rapidité d’acquisition. Si 
l’appareil a été vérifié (chapitre 3), cela veut dire que son défaut de jus- 
tesse (inconnu de l’utilisateur courant) ne doit pas dépasser une valeur 
seuil maximum (a = EMT : écart maximum toléré). 

-a<AJ < +a 

La méconnaissance de l’erreur de justesse invite à la prudence. On 
J retiendra comme fonction de densité de probabilité, une loi uniforme : 



M = M, 




© 


u c (M c J - u 





a ‘ 


3 
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Remarque : il y a une grande similitude avec le concept d'inertie du pro- 
cessus de contrôle Iq (proposé au chapitre 7) : 



l G nécessite un étalonnage pour être quantifiée. Pour améliorer l'inertie 
Iq du processus de contrôle, on peut chercher à améliorer les deux com- 
posantes, en acceptant des compensations. Pour améliorer u c , il faut en 
priorité intervenir sur l'incertitude de la méthode ou, en dernier recours, 
changer de classe d'appareil (réduire la valeur de a). 

Correction locale de dérive 

Ce cas est très proche du précédent, sauf qu’en plus, l’erreur de justesse est 
une fonction du temps. La solution la plus utilisée est d’essayer d’éliminer 
cette influence (voir chapitre 9). L’incertitude associée est contenue dans 
l’incertitude de justesse. 

Correction locale de grandeurs d'influence 

La méthode de traitement correspond à la méthode classique, où l’on 
cherchera autant que possible à soustraire les facteurs d’influence (sous 
forme de corrections algébriques c t ) dans la relation : 

= ~5 J+ 1c, 

u L S:yJ lrii ) = +s aj + X ^ 2 

À titre d’exemple, nous présentons la correction à faire lors d’une 
mesure directe avec un appareil qui n’est pas à la même température que 
la pièce fabriquée (fig. 8.9). En particulier, si la pièce vient d’être usinée, 
sa température interne peut-être nettement plus élevée que celle de 
l’appareil. 


I 

e 


o 

© 
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Figure 8.9 : Exemple de correction de température 


Etude du milieu 


I 23/05/2008 19:19 




àL rM 
àh AM 
dL P , 2 p. 
dÀ A 
ÔL h _ 2 


dO A 

SL P 

dA P 

SL „, 2 q . < 
dÛ P 


0,0012 
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Correction globale 


L’appareil retenu couvre une étendue de mesure importante : 



Par commodité de présentation, nous retiendrons Xi n f= 0 (mais ce n’est 
pas toujours le cas, par exemple pour un micromètre : 25-50 mm). 

Prenons le cas d’un micromètre 0-25 mm. À l’aide de cales-étalons de 
classe 1, on compare la valeur fournie par l’appareil (xÿ) et la valeur cor- 
respondante lue sur la cale-étalon y 0 i (l’exactitude de cette valeur sera 
supposée grande par rapport aux indications fournies par l’appareil). 
Pour chaque groupe de k mesures Xÿ, on calcule une estimation d’écart : 


Note : nous utiliserons la présentation de C. Perruchet et M. Pr/e/ 1 qui 
conseillent de retenir la donnée de l'instrument comme grandeur 
d'entrée et comme grandeur de sortie la valeur de l'étalon correspon- 
dante. Cette présentation permet d'obtenir des régressions linéaires 
par la méthode des moindres carrés, l'entrée (valeur lue) n'étant pas une 
variable aléatoire (exemple figure 8.10). 


Figure 8.10 : Exemple d'erreurs de justesse par rapport à l'étendue 



Cas A E Cas B E Cas C 


► Cas A : l’estimation de l’erreur de justesse est constante sur toute 
l’étendue (cas typique d’une lame d’air entre un marbre et une cale- 
étalon). 

► Cas B : au défaut précédent, s’ajoute un défaut d’échelle (par exemple, 
le pas du micromètre est un peu fort). 
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* Cas C : cas le plus général, un défaut de linéarité vient se superposer 
aux deux précédents (variation non constante du pas de vis du micro- 
mètre). 

Remarque : dans l'exemple précédent, nous avons considéré que la 
courbe S = f(E) est la même en charge ou en décharge. Il n’y a pas d'effet 
d'hystérésis. 

La représentation de la figure 8.10 est peu exploitable si l’on a une 
grande étendue relativement aux écarts estimés. Il est plus pratique de 
« redresser ce graphique » en représentant les écarts par rapport à la pre- 
mière bissectrice (appareil sans défaut). De plus, rappelons qu’il ne s’agit 
que d’estimations en nombre limité (fig. 8.1 1). 

Figure 8.11 : Exemple d'erreur de justesse pour un micromètre 
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Incertitude associée à une correction nulle 

Nous sommes dans le même cas que pour une justesse locale, son défaut 
de justesse (inconnu de l’utilisateur occasionnel) ne doit pas dépasser 
une valeur seuil maximum (a = EMT : écart maximum toléré), ceci 
quelle que soit la valeur d’entrée. 

-a<àJ{X.)<.+a VX 

On retiendra comme fonction de densité de probabilité, une loi uniforme : 
a 2 

u c ( - M cJ^ u M irMi (x,) + Y 
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Parfois, le bornage ± a se fait plage par plage. On a des limites en 
« marches d’escalier » (par exemple, instrument de pesage 1 en figure 8.12). 

Figure 8.12 : Exemple d'erreur tolérée en fonctionnement 
pour une balance de classe II d'étendue 30 kg 


Erreur maximale tolérée en service 



Incertitude associée à une correction spécifique 
à une valeur locale 

Cela nécessite d’avoir la courbe d’étalonnage avec l’incertitude locale 
associée (suivant la position dans l’étendue, nous n’avons pas les mêmes 
valeurs) : 

M mr = M bme -8 J {X i ) 

Vu la faible taille des échantillons, on peut éventuellement retenir une 
incertitude constante pour toute l’étendue de mesure, calculée avec 
l’ensemble des valeurs. 


NF EN 45501 novembre 1993, Aspects métrologiques des instruments de pesage à 
fonctionnement non automatique - Indice de classement : E12-000. 
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Bien que ce genre de courbe soit donnée avec chaque appareil, ce n’est 
pas d’un emploi pratique en production. Cette approche est plutôt 
réservée à la mesure en laboratoire. 

Figure 8.13 : Courbe d'étalonnage brute (non modélisée) 
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Incertitude associée à une correction modélisée 

C’est un compromis entre les deux méthodes précédentes. La correction 
est modélisée, cela évite de mémoriser un nombre important de valeurs. 
Par exemple, si l’on retient un modèle linéaire, il suffit de mémoriser 
deux paramètres pour faire les corrections (fig. 8.14). 


Figure 8.14 : Modélisation de la courbe d'étalonnage 


Erreur de justesse 


Cette approche est d’une mise en œuvre pratique et facile (tableur). 
Néanmoins, il faut rester très prudent. Il ne faut pas oublier que ce n’est 
qu’un modèle et qu’il ne colle pas toujours parfaitement à la réalité. On 
peut choisir un modèle plus complexe, mais il demandera la mise en 
mémoire de plus de paramètres. 
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Dans l’exemple de la figure 8.14, la correction proposée est du type : 
8 J =a-X i + b 


a et b ne sont que des estimations des vrais paramètres a et j3. L’incerti- 
tude va devenir importante dès que l’on s’éloigne de la valeur moyenne 
des X i (à peu près le milieu de l’étendue). Cette incertitude nécessite 
d’estimer un écart-type des résidus s R( 5 s idus (voir chapitre 6). Nous pou- 
vons proposer : 


M cor ~ M bme — a- X { — b 



2 


O 

© 
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Maintenir la qualité 
métrologique 



L’entreprise doit vérifier au cours du temps que la qualité de la mesure 
ne se dégrade pas. En cas de constat trop tardif de dérive excessive de 
l’appareil, il faut reprendre les mesures ayant permis de décréter à tort la 
conformité. Pour limiter ce risque, l’entreprise devra mettre en place et 
piloter un processus de surveillance pour certifier cette garantie de qualité. 
Il est tentant d’exiger une fréquence de surveillance élevée pour minimiser 
le risque client. Ce processus de surveillance nécessite des ressources maté- 
rielles et humaines et immobilise les appareils. Cela engendre des coûts 
supplémentaires qu’il faut mettre en regard avec les risques de non- 
détection (reprise des contrôles, rappels, etc.). Le choix de la fréquence 
de surveillance devra essayer de répondre à ce souci d’optimisation 
financière et de qualité. On commence souvent par une vérification à 
périodicité calendaire fixe, associée à un réajustage si nécessaire. Cela 
permet de voir comment se comportent les appareils. Avec plus d’expé- 
riences, on peut appliquer les méthodes employées en maîtrise statis- 
tique des procédés (MSP). En effet, il serait inutile et coûteux de réajuster 
systématiquement un appareil qui ne s’est pas déréglé (en fait où il n’y a 
pas confirmation de déréglage). Quelle que soit la stratégie retenue, elle 
doit être décrite dans une procédure précisant les méthodes utilisées pour 
déterminer ou modifier les intervalles de confirmation métrologique 1 . 

Il va de soi qu’un appareil qui aurait subi des problèmes de fonctionne- 
ment (mauvaise manipulation, chocs, surcharge, etc.) ne doit pas 
attendre la date prévue avant d’être vérifié à nouveau (certains appareils 
sont très sensibles aux déménagements et doivent être vérifiés systémati- 
quement). 

Rappelons qu’à l’issue de chaque vérification de conformité métrolo- 
gique, l’appareil peut être : 

» maintenu pour l’utilisation prévue (conforme) ; 

* ajusté ou réparé et remis en service ; 


1. NF EN ISO 10012 septembre 2003, Systèmes de management de la mesure - Exigences 
pour les processus et les équipements de mesure et FD X07-014 novembre 2006, 
Métrologie - Optimisation des intervalles de confirmation métrologique des équipe- 
ments de mesure. 
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» déclassé (utilisation nécessitant une précision plus faible) ; 

* mis au rebut. 

Les deux derniers cas relèvent de la maîtrise de produits non conformes, 
ils doivent être clairement identifiés. 


Correction de la dérive 
à périodicité calendaire fixe 

Cette tradition de périodicité calendaire est héritée de la métrologie légale 
(vérification des balances commerciales, des pompes volumétriques à 
essence, etc.). Cet état d’esprit est passé naturellement de la métrologie 
légale à la métrologie industrielle. L’OIML (Organisation internationale 
de la métrologie légale) publie le document n° 10 : ILAC- G24/OIML 
D10 : Conseils pour la détermination des intervalles de réétalonnage des 
équipements de mesure utilisés dans les laboratoires d’essais l . Ce document 
retient cinq méthodes de révision des intervalles de vérification : 

méthode 1 : ajustement automatique en escalier (inspirée des cartes 
CUSUM, on reporte sur un graphique une augmentation si l’on 
dépasse la valeur tolérée ou une diminution dans le cas contraire) ; 
méthode 2 : utilisation des cartes de contrôles de la MSP (cette 
méthode permet d’appréhender les dérives, lentes ou rapides). Elle 
est bien adaptée pour les entreprises maîtrisant ces outils de contrôle 
de la qualité par les statistiques (nous reviendrons sur cette méthode 
au paragraphe 2) ; 

méthode 3 : on reprend les principes précédents mais en se basant 
sur une durée réelle d’utilisation. Cela nécessite un système ou une 
technique de mesure du temps de fonctionnement. Cette méthode 
n’est pas forcément judicieuse lorsque les cycles de fonctionnement 
sont très courts ; 

méthode 4 : utilisation systématique d’une pré-vérification (appelée 
test « boîte noire »). Bien adaptée aux instruments complexes, cette 
méthode propose à chaque mise en route une vérification intégrée 


1. Téléchargeable sur le site : www.oiml.orglpublicationslDID010-e07.pdf. 
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simplifiée qui tire la sonnette d’alarme en cas de soupçon de dérive 
importante exigeant un complément de vérification (maintenance 
prédictive). Il faut que la pré-vérification soit potentiellement sen- 
sible à toutes les causes de dérives possibles ; 

méthode 5 : approches statistiques basées sur les observations d’un 
grand nombre d’instruments. 


Choix arbitraire d'initialisation de la périodicité 

L’empirisme est souvent de règle au début pour choisir la période de 
correction (par exemple une période de six mois, recommandation du 
fournisseur d’appareils...). Elle est affinée suite aux constats de dérive 
ou non lors des vérifications. On peut diminuer la fréquence si les enre- 
gistrements indiquent des dérives inexistantes ou lentes. 

La périodicité des ajustages (si nécessaire) est fonction entre autres : 
du degré d’exactitude souhaité : en général les appareils de basse 
résolution (comme un réglet) dérivent peu et nécessitent moins 
d’attention. On peut accepter une incertitude de l’équipement plus 
importante ; 

du taux d’utilisation : certains appareils sont très spécifiques 
(comme des piges de contrôle), ils peuvent rester plus ou moins long- 
temps sans être utilisés. On peut ne faire démarrer la période d’utili- 
sation qu’à partir de la date de mise en service, sous réserve que cette 
dernière soit bien identifiée (pige recouverte de cire par exemple) ; 
de l’environnement d’utilisation : si l’instrument est utilisé dans 
un environnement difficile (chaleur, vibrations, etc.), il est évident 
qu’il faudra réduire la période de vérification, de même si le nombre 
d’opérateurs utilisant l’appareil est important. Lorsqu’il n’y a qu’un 
seul opérateur expérimenté et soigneux, cela crée un sentiment 
d’appartenance bénéfique pour l’appareil ; 

de l’aptitude à dériver des appareils : par exemple des équerres uti- 
lisées en mécano-soudage peuvent être vérifiées seulement tous les 
trois ans (sauf accident d’utilisation) ; a contrario certains appareils 
peuvent présenter des dérives rapides et nécessiter des vérifications 
quotidiennes (processus d’équilibrage par exemple) ; 
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du risque économique et/ou d’insécurité encouru : dans l’indus- 
trie spatiale, le coût d’une vérification sera négligeable en regard du 
coût d’un lanceur. La vérification sera alors systématique. Pour des 
pièces de sécurité, une campagne de rappel deviendrait nécessaire en 
cas de doute ; 

» etc. 


Principes de calcul des périodes de vérification 

La détermination de la période inter- vérifications nécessite (fig. 9.1) : 

» de définir ou calculer un écart maximum toléré (EMT) pour la dérive 
(limite d’exactitude retenue) ; 

► d’avoir un modèle de dérive en fonction du temps : d - f(t). C’est 
souvent un modèle linéaire qui est retenu à l’aide d’une régression à 
partir de points de contrôle. 

Figure 9.1 : Période de dérive limite en fonction 
de l'incertitude élargie U 



Pour calculer la valeur limite, il faut faire une hypothèse de point de 
départ (par exemple, d( 0) = 0). Par prudence, comme la vérification de 
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la dérive n’est pas exempte d’incertitude ( U : incertitude élargie, sup- 
posée a priori constante en fonction du temps, vue au chapitre 8), on 
diminuera par sécurité la période. 

EMT -U = f {Période) 

Période = f~\EMT -U) 

Soit avec un modèle linéaire : d — k-t : 


Période — 


EMT -U 

k 


Par exemple pour les contrôles par calibre, si l’on a constaté que l’usure 
était bien linéaire, cette technique s’applique bien. On peut même 
l’affiner si l’on a estimé la dérive d§ initiale pour l’appareil neuf : 


Période — 


EMT — U — d 0 
~k 


Cette approche n’est applicable que lorsque l’incertitude n’est pas trop 
forte en regard de Y EMT. 


Choix de la périodicité par la méthode Opperet 

Le Collège français de métrologie a proposé une méthode originale 
d’optimisation des périodicités d’étalonnage appelée « Opperet 1 » (pour 
« optimiser les périodicités d’étalonnage »). La méthode 2 est une solu- 
tion ouverte, intégrée dans plusieurs logiciels de gestion des moyens de 
mesure. Cette méthode présente le gros avantage de prendre en compte 
des critères techniques et économiques. Sa mise en œuvre est basée sur la 
démarche suivante. 

Prise en compte des facteurs d'influence 

La méthode propose, comme suggestions non limitatives, neuf thèmes 
de réflexion précisant les influences économiques, d’utilisations, de ris- 
ques. À partir de ces items de réflexions, on définira un certain nombre 


1 . www. cfmetrologie. com. 

2. Cette méthode est également présentée dans l’annexe D du fascicule de documentation 
FD X07-014 novembre 2006. 
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de critères Cj de jugement pour les appareils appartenant à un même 
périmètre (c’est-à-dire concernés par le même choix de critères). Ce 
choix de critères doit être consensuel et faire intervenir différents ser- 
vices (analogie avec l’AMDEC). Chaque instrument i appartenant au 
périmètre sera noté pour chaque critère (note Nÿ pour le critère Cj. Il 
est proposé de retenir une base de notation comprise entre — 2 et + 2. 
Prenons comme exemple didactique un périmètre constitué de cinq 
appareils (Aj à A5) pour lesquels, on a retenu quatre critères (C 1 à C4 ; 
normalement la méthode prévoit au moins cinq critères) (tableau 9. 1). 


Tableau 9.1 : Exemple simplifié de notation 


Périmètre 1 

Critères 

Cl 

C2 

C3 

C4 

* 

Ai 

Nu 

N12 

N, 3 

N,4 

a 2 

N21 

N 2 2 

N23 

N24 

a 3 

N31 

N 32 

N33 

N34 

A4 

N41 

N42 

N43 

N44 

A 5 

N51 

N 52 

n 53 

N 54 


Notation de l'appareil 

Pour éviter un problème de notation, on propose pour chaque critère Cj 
un changement de variable transformant la note Nÿ de l’appareil i en un 
écart normalisé ENÿ (qui est analogue à la variable centrée réduite u de 
la loi normale, ou t de la loi de Student) : 


s.j est la valeur de l’écart-type expérimental pour l’ensemble des notes 
correspondant au critère Cj. 

Chaque appareil i se verra affecté d’une note globale NG- en utilisant un 
système de pondération des critères (classement de type ABC). 


i 


I 

e 


o 

© 
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Tableau 9.2 : Calcul des notes globales (pondération des critères) 


Périmètre 1 

Critères 

Ci 

C 2 

c 3 

c 

■4 


Ai 

NGi 


EN, i 


EN, 2 


EN, 3 


EN, 4 


a 2 

ng 2 


en 21 


en 22 


en 23 


en 24 

1 

a 3 

ng 3 

Pl 

EN 3 i 

p 2 

en 32 

P 3 

en 33 

P4 

en 34 

< 

a 4 

ng 4 


EN 4 i 


en 42 


en 43 


EN 44 


a 5 

ng 5 


EN 5 i 


en 52 


en 53 


EN 54 


Cette note globale est à nouveau transformée en écart normalisé suivant 
la même procédure, soit pour l’appareil i : 


NG i - NG 

S NG 


Tableau 9.3 : Calcul des écarts normalisés globaux 


Périmètre 1 


Al 

EN ng _, 

a 2 

EN ng _ 2 

A3 

EN ng _ 3 

A4 

EN ng4 

A 5 

EN ng _ 5 


Distribution des périodicités 

La périodicité initiale a été choisie au pire des cas pour tous les appareils 
d’un même périmètre (pratiquement valable pour tous les appareils, elle 
devient trop excessive pour les appareils plus fiables). La méthode 
Opperet a comme objectif d’individualiser cette périodicité en fonction 
de l’écart normalisé global EN^g i attribué à chaque appareil. On part 
de l’hypothèse réaliste (théorème central limite) que la distribution des 
périodicités pour l’ensemble des appareils du périmètre est normale. 
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Rappelons que ce modèle dépend uniquement de deux paramètres (un 
de tendance, la moyenne et un de dispersion, l’écart-type). Donc, pour 
pouvoir l’appliquer, il nous faut deux informations indépendantes. Elles 
nous seront fournies indirectement par le choix de probabilités qui nous 
permettront de calculer deux variables réduites u in ft t u§ u p (fig. 9.2). 

Figure 9.2 : Présentation des distributions de périodicités 



* Pour la périodicité actuelle, on pourra vérifier le pourcentage d’appa- 
reils qui étaient non conformes (cela suppose que l’on a des enregis- 
trements). Le plus souvent cette probabilité est estimée entre 2 % et 
3 %. Cela donne : 

p mf = 2% -> u M = -2,054 
P'm£ — 3% — > u- m ç — —1,881 
Nous retiendrons pour u in f la valeur arrondie 2. 

► La périodicité maximum est évaluée sur des critères de bon sens (par 
exemple, si l’on a des observations, on peut la définir à l’aide des 
appareils du périmètre qui n’ont jamais nécessité un réajustage). 
Cette valeur est souvent plafonnée à une limite raisonnable. 

La probabilité d’avoir une périodicité supérieure à cette limite est for- 
tement improbable. Par exemple, si l’on retient : 

Psup ~ 99,8% — » u Sup = +2,878 
Psup ~ 99,9% — » u Sup - +3,090 
La valeur u§ U p sera arrondie à 3. 
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En résumé, on appliquera la méthode suivante (fig. 9.3) : 

1. Saisir la période actuelle. 

2. Définir la période maximum raisonnable. 

3. En déduire « 7 - w yen fonction du risque associé à la période actuelle (en 
général u i n f = ~ 2). 

4. En déduire u^ u p en fonction de la probabilité d’avoir une période 
inférieure à la période maximum (en général u$ u p = +3). 

5. Calculer en variable centrée réduite u Sup -u m f (pris à 5 par défaut). 

6. En déduire s Périodictté : 


7. En variable centrée réduite la moyenne est par définition nulle. 

8. En variable physique (temps) la périodicité moyenne est égale à : 
~ p +2-s 


Figure 9.3 : Calcul de la moyenne et de la dispersion de la périodicité 
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Personnalisation de la périodicité d'étalonnage 

En admettant que la distribution de la variable centrée réduite u et la 
distribution de l’écart normalisé global ENq sont identiques, on peut 
écrire : 


EN N G_i 

EN m _ i 


P— P 

i Moyenne 
S Périodicité 


P - P + EN ■ s 

1 Moyenne ~ J - dy NG_i J Périodicité 

La différence entre la nouvelle périodicité et la périodicité actuelle 
Pl Actuelle s’appelle le « correcteur Opperet ». Le plus souvent, il est 
positif, c’est-à-dire que l’on gagne en périodicité. Par contre, dans le cas 
de notation sévère pour les critères Cj, ce correcteur peut-être négatif. 


Note 1 : s'il existe dans l'entreprise une méthode propre de détermina- 
tion des périodicités, il est toujours intéressant de voir s'il y a une corré- 
lation entre les périodicités retenues et les écarts normalisés calculés 
par la méthode Opperet. 

Note 2 : cette méthode demande une notation des critères attentive, 
non partisane, faite par des experts provenant de services différents. Il 
ne faut pas oublier ou sous-estimer un critère important. 


Correction de la dérive 
en utilisant les méthodes de la MSP 

Rappel sur la MSP 1 

La sortie d’un processus de production est sujette à des variations pou- 
vant provenir de causes communes ou de causes spéciales. 

* les causes communes correspondent à un bruit de fond du processus. 

On n’arrive pas à les maîtriser ni à les identifier de façon précise. Il faut 1 
« faire avec ». Elles engendrent la variabilité intrinsèque du processus ; 

o 

1 . Pillet M ., Appliquer la maîtrise statistique des processus, Paris, Éditions d’Organisation, 2005. ® 
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► les causes spéciales, difficilement prévisibles, vont perturber fortement 
la sortie du processus. Cela se manifestera soit par un déréglage 
(variation de la moyenne), soit par un accroissement de la dispersion 
(augmentation de l’écart-type). Le pilote du processus devra identi- 
fier rapidement la cause de cet aléa. Pour cela, il travaillera sur un 
échantillon de sorties afin de diminuer l’effet de la dispersion rési- 
duelle (fig. 9.4) 

Figure 9.4 : Dispersion et déréglage d'un processus 


Cible 

[ Déréglage 

1 valeur ^ 

déréglage J / 

dispersion [ / 

i • • i 

k Er 

s. de 

\ >'« 


. Dispe 

n* 



À partir de la moyenne d’un échantillon, le pilote du processus doit 
prendre une décision : 

► soit considérer que l’écart entre la moyenne de l’échantillon et la 
cible reste faible et imputable aux causes communes ; 

* soit estimer que cet écart paraît inadmissible, qu’il y a de fortes pré- 
somptions de causes spéciales. 

La moyenne de l’échantillon donne une estimation de la moyenne réelle 
du processus. Il reste un certain flou dans cette estimation (chapitre 6). 
Avec un risque a d’erreur (risque fournisseur), nous pouvons donner un 
intervalle de confiance contenant les moyennes d’échantillons de taille 
n, pour un processus parfaitement centré sur sa cible. 

Prob Cible — u / • < X < Cible + u / • -j= = l-a 

L ~ /2 yjn ~/2 yjn J 

On supposera que la variabilité naturelle (7 (causes communes) est par- 
faitement identifiée (observation de longue durée). 
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Le fait de travailler avec une moyenne permet d’avoir des limites natu- 
relles plus fiables pour détecter une dérive éventuelle (réduction en 1/V ri) 

(fig. 9.5). 


Figure 9.5 : Pilotage par moyenne ou par valeur individuelle 



Shewhart 1 a proposé de représenter (fig. 9.6) la double inéquation ci- 
dessus par une carte graphique de pilotage permettant une visualisation 
aisée des dérives (réglage à l’aide de la moyenne ou de la médiane) et des 
dispersions (indicateur avec l’étendue ou l’écart-type). 


I 

e 


1. Shewhart Walter A. (1 891-1967), Statistical Methodfrom the Viewpoint ofQuality Control. §" 

NF X06-031-1 décembre 1995, Application de la statistique - Cartes de contrôle - o 

Partie 1 : cartes de contrôle de Shewhart aux mesures. 
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Figure 9.6 : Exemple simplifié d'une carte de contrôle 


Date 

23/02 

Heure 

6 h 00 

6 h 30 

7 h 00 

7 h 30 

8 h 00 

8 h 30 

8 h 35 

9 h 00 

9 h 30 

Mesure 1 

1 

- 1 

1 

0 

1 

0 

-2 

- 1 

-2 

Mesure 2 

1 

0 

0 

0 

1 

2 

1 

1 

- 1 

Mesure 3 

-2 

-2 

- 1 

-3 

0 

3 

- 1 

1 

- 1 

Mesure 4 

1 

0 

- 1 

- 1 

1 

2 

- 1 

- 1 

0 

Mesure 5 

0 

0 

0 

0 

0 

3 

2 

0 

0 

Total 

1 

-3 

- 1 

- 4 

3 

10 

- 1 

0 

- 4 

Moyenne 

0,2 

-0,6 

-0,2 

-0,8 

0,6 

2,0 

-0,2 

0 

-0,8 

Etendue 

3 

2 

2 

3 

1 

3 

4 

2 

2 



Application de la MSP au suivi de vérification 

J Le processus de contrôle peut être piloté comme un processus de pro- 
I. duction. À la consigne d’entrée (mesurande), donnée par l’étalon, 

§- correspond une « sortie mesure » qui va présenter des dispersions, 

ô Soit ces dispersions sont dues uniquement à la variabilité résiduelle 
intrinsèque, soit il y a une cause spéciale qui provoque des variations 
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inhabituelles. La connaissance de la variabilité résiduelle demande 
une longue observation du processus qui fait défaut ici. Par défaut à 
l’initialisation, on peut la quantifier à l’aide de ou de Ggy si 

c’est toujours le même opérateur qui fait les étalonnages (voir 
chapitre 7 pour les études GRR de capabilité des processus de 
contrôle). 

La fréquence de contrôle étant faible comparée à celles utilisées en pro- 
duction, nous utiliserons une carte « moyennes et écarts-types » plutôt 
qu’une carte « médianes et étendues ». 

Les limites de contrôle, permettant de savoir si un ajustage est néces- 
saire, sont basées sur les inéquations suivantes : 

Prob \cible - <X< Cible + = 1 -a 




1 -a 


Par exemple avec n = 5 et un risque a de 2 %o (c’est-à-dire d’ajuster un 
appareil qui n’en avait pas besoin), nous obtenons : 

Prob Cible - 3.09 • <X< Cible + 3.09 • = 99.8% 

L S S J 



soit : 


Pïob[OMr-l,382 • g grr <X< Cible + 1,382 • a GRR ] = 99.8% 
Prob [0, 1 5 1 • o GRR < S < 2, 149 • o GRR ] = 99.8% 


Pour l’ écart-type S, on surveille surtout la borne supérieure, étant donné J 
le faible effectif de l’échantillon (fig. 9.7). à 

CL 

Ô 

© 
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Figure 9.7 : Exemples de cartes avec <r GRR de 0,01 mm 
et un échantillon de cinq mesures de l'étalon 


0,020 ■ 

0,015 . 

0,010 ■ 

0,005 ■ 

0,000 ■ 

- 0,005 ■ 

-0,010- 
-0,015- 

- 0,020 ■ 

0,025 

0,020 ■ 

0,015 ■ 

0,010 ■ 

0,005 ■ 

0,000 ■ 

01 23456789 10 11 12 

Les règles de base d’interprétation des cartes se feront comme pour le 
pilotage d’un processus classique de production (par exemple, « point 
hors limite, sept points consécutifs d’un même côté, sept points consé- 
cutifs tous en augmentation ou en diminution... », on prend alors la 
décision de réajuster le processus de contrôle). 

La carte de Shewhart permet de voir l’évolution de la justesse en fonc- 
tion du temps. Suivant le VIM 1 , c’est le concept de constance (nommée 
stability pour le MSA 2 ). Elle est particulièrement performante dans le cas 
de dérive rapide (comme un choc par exemple), par contre pour les 
dérives lentes, la carte EWMA est mieux adaptée. 


Carte des moyennes 


Ô 1 2 3 I 5 ] 7 8 9 10 ÏT 

? 6 


Cartes des écarts-types 


1 . Les termes spécifiques à la métrologie (comme grandeur ou unité) sont définis dans le 
VIM ( Vocabulaire international de la métrologie) — Projet de norme X07-001PR (PR 
NF ISO 99999 octobre 2006) . 

2. Automotive Industry Action Group (ALAG) (2002), Measurement Systems Analysis 
Reference Manual, 3 e éd., Chrysler, Ford, General Motors Supplier Quality Require- 
ments Task Force (www.carwin.co.uk/qs). 
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Carte EWMA 1 

Nous décrirons simplement l’esprit de la méthode, pour une application 
concrète le lecteur est renvoyé aux ouvrages spécialisés. 

Mesurons à l’aide d’un appareil un étalon de longueur x 0 . Nous allons 
obtenir une mesure légèrement différente X;. À l’aide de ces deux valeurs, 
nous déterminons un écart : 

e i = x i ~x o 

C’est cet écart que nous allons surveiller. L’écart réel peut être estimé à 
partir d’une seule mesure ou à partir d’un échantillon de n valeurs e^. La 
seconde solution est plus onéreuse mais meilleure quant à la confiance 
de la décision à prendre. 

Les cartes EWMA ( Exponentially Weighted Moving Averages) semblent 
particulièrement bien adaptées à la gestion des moyens de mesure lors 
de dérives lentes (la moyenne est pondérée par les valeurs antérieures) : 


Mq = valeur cible (valeur initiale prise égale à 0). 

A dépend des risques encourus (risques fournisseur a et client /3). Il est 
souvent pris égal à 0,2. Une valeur plus faible (A = 0,1) prend plus en 
compte les valeurs précédentes, elle permet de mieux détecter de petits 
écarts mais est moins performante pour détecter un décalage instantané 
important (A= 0,4). 

On détectera une cause spéciale, lorsque l’indicateur M- franchit une 
limite (inférieure LlC^i ou supérieure LSCj^-). Lorsque i croît, ces 
limites convergent assez vite vers les valeurs suivantes (en retenant un 
facteur d’élargissement égal à trois) : 


1. NF X06-031-3 décembre 1995, Application de la statistique - Cartes de contrôle - 
Partie 3 : cartes de contrôle à moyennes mobiles avec pondération exponentielle 
(EWMA). 


Af . = A • i?) + (1 - A) • M { _ 
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x 

\2-X) 


,{2-X) 

(Pour les échantillons de « une mesure » prendre n = 1 .) 


Carte CUSUM 1 

Il existe plusieurs types de cartes CUSUM ( Cumulative SUMs ). Le 
principe est de faire la somme des écarts précédents. Si cette somme 
dépasse une limite H, cela détectera une dérive probable. Comme la 
dérive peut être positive ou négative, on calculera deux sommes (tou- 
jours positives) pouvant être comparées à la même limite H : 

SH j = Afof[o,(S// y _ 1 + ^. 

SL = Max [0, (SZ._j - ë i - k ■ cr_ )] 

H -h-G 7 ( h est choisi entre 4 et 5, k = 0,5 pour détecter un décentrage 
d’un écart-type.) 


1. NF X06-031-4 décembre 1995, Application de la statistique - Cartes de contrôle - 
Partie 4 : Cartes de contrôle des sommes cumulées {CUSUM). 
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La démarche Six Sigma a permis de passer d’une production de masse 
avec contrôle final à une production optimisée avec contrôle intégré. 
L’objectif principal de cette démarche est de faire la chasse à la disper- 
sion. L’approche processus, fortement préconisée par la norme ISO 
9001 (v 2000), permet de bien identifier où mesurer la dispersion. Les 
sorties du processus demandent bien sûr d’être vérifiées (conformité 
produit), mais on met de plus en plus l’accent sur la maîtrise et le 
pilotage. 

Une des solutions pour mettre en œuvre la démarche Six Sigma, consiste 
à appliquer la méthode DMAICS, présentée à la figure 10.1. 

L’étape « Mesurer » arrive en deuxième position. En effet, si la qualité 
de l’information saisie est douteuse, ce sont toutes les étapes suivantes 
qui seront perturbées. Cette étape doit être robuste. Le processus de 
mesurage nécessite un investissement en ressources humaines et maté- 
rielles. En particulier, la dispersion inhérente au processus de contrôle 
doit être faible vis-à-vis de celle du processus à maîtriser et à surveiller. 
Par exemple, il est utopique de vouloir prétendre faire du SPC ( Statis - 
tical Process ControJ) si le processus de mesure n’est guère plus perfor- 
mant que le processus de production à piloter (en supposant que ce 
dernier est capable vis-à-vis des exigences du client mais sans excès). Le 
remède risque d’être pire que le mal. 

Pour cela, nous avons vu qu’il faut vérifier : 

s l’adéquation processus de contrôle-processus de production (à l’aide 
d’une étude de répétabilité et de reproductibilité) ; 

* la « justesse » de l’information en raccordant les appareils aux étalons 
nationaux, 

► la pérennité de la performance des appareils en utilisation (gestion et 
suivi). 
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Figure 10.1 : Les étapes et revues de projet de la méthode DMAICS 1 



R5 


Ces différentes étapes doivent être compatibles avec un objectif de coût. 
Faire des mesures de qualité, suppose que l’on obtient une information 
pertinente pour les décisions futures à prendre, mais que le degré d’exac- 
titude recherché reste raisonnable en rapport aux exigences du client et 
des dispersions des autres composantes de la chaîne de production. 


I 

e 


1 . In Pillet M., Six Sigma - Comment l’appliquer, Paris, Éditions d’Organisation, 2003. 


o 
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